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INTRODUCTION 


This  report  describes  the  results  of  a  research  program  extending  from 
1  December  1977  to  30  September  1979*  The  objectives  of  the  program  have 
been  to  study  the  relationships  between  geomagnetic  pulsation  activity  in 
the  period  range  10  to  300  seconds  (Pc  3,4,  and  5  by  traditional  classifi- 
cation)  and  solar  wind  parameters.  The  study  has  centered  on  two  particular 
parameters,  solar  wind  speed  and  interplanetary  magnetic  field  (IMF)  cone 
angle  6^  (the  angle  between  the  IMF  and  the  sun-earth  line)  as  representing 
the  critical  quantities  in  two  models  of  pulsation  origin. 

According  to  the  first  model  the  properties  of  the  bow  shock,  and  in 
particular,  turbulence  in  the  downstream  magnetosheath,  are  controlled  by 
the  solar  wind  velocity  and  field  direction.  As  these  waves  are  swept  around 
the  magnetosphere  they  penetrate  or  stimulate  the  boundary  and  excite  internal 
resonant  field  line  oscillations.  These  are  observed  at  the  ends  of  the 
resonant  field  lines  as  ground  magnetic  pulsations. 

In  the  second  proposed  mechanism,  magnetosheath  flow  along  the  magneto¬ 
pause  excites  the  Kelvin-Helmhol tz  instability  producing  surface  waves.  These 
evanescent  waves  couple  internally  to  resonant  field  tines  and  produce  ground 
pulsations.  In  this  mechanism  the  instability  criterion  is  a  function  of 
solar  wind  velocity  and  field  orientation. 

In  this  program  statistical  correlations  between  solar  wind  parameters 
at  Explorer  35  and  two  ground  observatories  reached  fruition.  The  amplitude 
of  ground  pulsation  increased  as  the  velocity  increased,  and  was  larger  when 


the  field  was  aligned  with  the  earth-sun  line.  These  results  have  been 
confirmed  in  several  ways  by  researchers  at  other  institutions  as  well. 

The  ultimate  purpose  of  the  investigation  is  to  attain  the  capacity 
to  use  micropulsation  records  acquired  from  surface  magnetometers  to  infer 
certain  key  parameters  of  the  solar  wind,  such  as  field  magnitude,  field 
direction,  wind  velocity,  and  of  the  magnetosphere,  such  as  magnetopause 
distance  and  subauroral  plasma  density.  The  present  investigation  has  sug¬ 
gested  this  is  a  realizable  goal. 

The  following  pages  discuss  the  results  in  a  general  way  in  terms  of 
activities  carried  out  under  the  initial  work  statement  and  specific  scien¬ 
tific  findings.  Technical  descriptions  of  results  appear  as  Appendices 
consisting  of  copies  of  published  or  submitted  papers.  Additional  results, 
pending,  but  incomplete  or  undocumented,  are  described  briefly  and  the 
status  of  the  project  is  summarized  and  recommendations  offered. 


RESULTS 


Statement  of  Work.  The  following  statements  summarize  the  activities  of  this 
program  In  terms  of  the  work  statement  In  the  proposal  of  August  1977 - 

a.  Study  the  relationship  of  solar  wind  parameters  and  wave  propa¬ 
gation  and  fluid  flow  In  the  magnetosheath  to  micropulsation 
exci tat  ion. 

Theoretical  computations  of  plasma  parameters  in  the  magnetosheath  were 
examined  as  a  guide  to  the  locations  and  local  times  the  magnetopause  is  most 
likely  to  be  affected  by  convection  and  propagation  of  large  amplitude  oscil¬ 
lations  from  the  shock.  Also,  the  literature  on  Kelvin-Helmhol tz  excitation 
of  the  magnetopause  was  reviewed  to  clarify  whether  that  phenomenon  consti¬ 
tutes  a  competing  or  a  complementary  source  of  geomagnetic  pulsations,  and 
to  determine  the  extent  to  which  the  K-H  mechanism  has  been  tested  by  obser¬ 
vation.  Factors  governing  the  resonance  of  waves  in  the  magnetosphere  were 
reviewed  to  assess  the  possible  effects  on  ground  station  data.  Finally, 
upstream  wave  and  energetic  particle  components  of  the  quasi-paral lei  (pul¬ 
sation)  shock  structure  were  studied  in  relationship  to  the  INF. 

b.  Improve  the  methodology  employed  In  handling  surface  and  satellite 
data  appropriate  to  the  Study. 

The  use  of  hourly  distributions  of  interplanetary  field  data  was  developed  and 
applied;  joint  distributions  of  pulsation  amplitude,  solar  wind  velocity,  and 
IMF  cone  angle  were  studied;  the  representation  of  the  cone  angle  by  its 
cosine  was  explored  and  found  to  improve  the  organization  of  the  data;  an 
approach  to  normalizing  ground  station  data  to  produce  a  global  measure  of 
pulsation  activity  independent  of  local  time,  latitude,  and  spectral  distortions 
was  designed  and  proposed  for  future  study. 
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c.  Exploit  the  schedules  and  formats  of  the  available  surface  and 
satellite  data  for  further  correlational  analysis. 

Correlations  between  solar  wind  speed  were  expanded  to  include  a  second  data 
interval,  a  second  ground  station,  and  a  range  of  longer  period  pulsations 
(Pc  S)  than  had  previously  been  investigated. 

d.  Assess  the  prospects  for  initiating  an  operational  system  using 
ground  station  micropulsation  records  to  estimate  selected  solar 
wind  parameters. 

Evaluations  were  made  of  the  methods  necessary  to  turn  data  into  an  operational 
index  (see  b.  above)  of  the  AFGL  network's  potential  contributions  to  an  opera¬ 
tional  system  and  its  relationship  to  other  networks,  of  the  state  of  knowledge 
of  the  sources  of  daytime  pulsations,  and  of  the  prospects  for  employing  pul¬ 
sations  as  event  predictors. 

e.  Publish  significant  results  in  the  scientific  literature  as  they 
are  developed. 

Six  reports  were  published,  completed,  or  submitted  for  publication  during 
the  performance  interval.  Three  international  meetings  and  workshops  were 
attended  and  results  presented  at  each. 

Specific  Accomplishments.  It  was  established  that  a  definite,  strong  correla¬ 
tion  prevails  between  pulsation  amplitude  and  solar  wind  speed  for  Pc  3 
and  4  (Appendix  A). 

The  velocity  correlation  also  applies  to  Pc  5  (Appendix  E) . 

The  correlation  between  IMF  cone  angle  6^^  and  pulsation  amplitude  is 
strongly  dependent  on  transient  excitation  of  pulsation  signals,  so  that  in 
dealing  with  data  in,  say,  hourly  units  of  time,  the  distribution  of  0„R  within 


the  hour,  rather  than  Its  average,  should  be  taken  into  account  for  best 
results  (Appendix  B) . 

Pulsation  amplitude  Is  jointly  dependent  on  both  0Xfi  and  V^;  maximal 
activity  occurs  when  both  0XB  is  small  and  V^y  is  large  (Appendix  C) . 

The  observations  supporting  either  the  direct  excitation  of  pulsations 
by  waves  in  the  magnetosheath  or  the  stimulation  of  pulsations  by  surface 
wave  instability  (Kel vin-Helmhol tz)  are  incomplete  at  best,  leaving  appreciable 
gaps  in  models  attempting  to  connect  signals  at  the  magnetopause  with  signals 
at  the  ground  (Appendix  D) . 

The  quasi-parallel  portion  of  the  bow  shock  is  intimately  associated 
with  hydromagnetic  waves  and  suprathermal  particles  in  the  solar  wind  outside 
as  well  as  inside  the  shock.  This  association  will  have  to  be  taken  into 
account  in  application  of  pulsation  measurements  as  diagnostics  of  plasma 
properties  outside  the  magnetopause  (Appendices  F,G). 

The  AFGL  magnetic  observatory  system  may,  by  virtue  of  site  locations 
and  data  handling  capability,  have  the  virtue  of  recording  measurements  ideally 
suited  to  investigating  and  monitoring  a  part  of  the  pulsation  band  diagnostic 
of  magnetospheric  and  solar  wind  conditions,  but  considerable  improvement  in  the 
data  processing  and  distribution  will  be  required  (Appendix  H) . 

PENDING  RESULTS 

Several  incipient  or  incomplete  investigations  bear  on  the  objectives 
of  this  study.  Two  results  drawn  from  these  investigations  are  sufficiently 
advanced  to  be  described  here. 
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Magnetosheath  Wave  Patterns.  Figure  1  Illustrates  the  relative  influence  of 
two  principal  factors  governing  the  delivery  of  large  amplitude,  quasi -paral lei 
shock-structured  waves  to  the  magnetopause.  These  factors  are:  solar  wind 
velocity  (speed  and  direction)  and  Alfven  wave  velocity.  As  the  streamlines 
in  the  figure  show,  only  solar  wind  passing  through  the  subsolar  region  of 
the  shock  impinges  directly  on,  or  comes  very  close  to,  the  magnetopause. 

Waves  convected  with  the  wind  on  the  flank  of  the  shock  can  reach 
the  dayside  magnetopause  only  if  they  cross  the  sheath  in  less  time  than  the 
plasma  carrying  them  takes  to  sweep  them  downstream.  The  small  pairs  of 
arrows  indicate,  at  the  intersection  of  each  pair,  the  relative  velocities 
of  the  solar  wind  in  its  local  direction  of  flow,  and  the  Alfven  wave  velocity 
toward  the  magnetopause.  The  flow  and  wave  velocities  have  been  estimated  from 
the  computations  of  Spreiter  and  Alksne  (1969)  for  the  IMF  in  the  ecliptic 
(the  plane  of  the  figure)  at  45°  to  the  sun-earth  line.  The  small,  circled 
numbers  in  the  magnetosheath  in  the  figure  indicate  approximately  the  local 
ratio  CA/VSW  (i.e.,  M^  *)  where  is  the  Alfven  speed  (M^  is  the  Alfven  Mach 
number) . 

It  is  clear  in  Figure  I  that  not  only  does  plasma  entering  the  flank 
of  the  magnetosheath  avoid  the  magnetopause,  but  waves  originating  in  the  flank 
plasma  are  so  slow  relative  to  Vgy  that  they  are  unlikely  to  propagate  to  the 
magnetopause  until  far  downstream.  In  contrast,  subsolar  plasma  not  only 
streams  toward  the  magnetopause,  but  its  parametric  properties  support  faster 
wave  propagation  than  at  the  flank,  making  delivery  of  Waves  to  the  magneto¬ 
pause  highly  probable.  The  figure  suggests  that  up  to  some  angle  from  the 
subsolar  point,  quasi-parallel  waves  could  cross  layers  of  solar  wind  and 


reach  the  dayside  magnetopause,  while  beyond  that  angle  they  could  not.  The 

numbers  on  the  streamlines  just  inside  the  nominal  shock  outline  are  rough 

estimates  of  the  length  of  the  corresponding  streamline  segment  up  to  the 

terminator  (the  -Y  axis)  divided  by  the  shortest  distance  from  the  streamline 

to  the  magnetopause.  These  numbers  are  a  crude  guide  to  the  ratio  of  V-.,  to 

5W 

CA  needed  to  pass  shock  oscillations  to  the  magnetopause.  If  the  distance 
ratio  is  greater  than  the  velocity  ratio  (the  inverse  of  the  circled  numbers), 
the  waves  may  reach  the  magnetopause;  otherwise  they  will  not.  The  figure 
suggests,  then,  that  waves  reaching  the  magnetopause  must  originate  in  the 
shock  no  more  than  about  kS°  or  so  from  the  subsolar  point,  when  the  IMF  is 
at  A5°.  Such  a  requirement  is  compatible  with  the  model  of  shock  origin, 
since  a  A5°  IMF  will  cause  quasi-parallel  structure  to  prevail  in  the  sub¬ 
solar  region. 

Further  evaluation  of  Figure  1  and  its  implications  for  other  dir¬ 
ections  of  the  IMF  are  in  progress  as  this  report  is  being  prepared  and  will 
be  released  and  discussed  in  coming  months. 

Correlations  at  Synchronous  Elevations.  All  the  published  correlations  of 
this  study  so  far  have  compared  interplanetary  with  ground  station  data.  The 
phenomena  of  interest,  however,  are  magnetospher ic  waves,  of  which  traditional 
micropulsations  at  the  earth's  surface  are  only  one  manifestation.  An  ad¬ 
ditional  approach  to  tying  the  magnetosphere's  oscillations  to  the  solar  wind 
is  through  satellite  measurements,  of  which  those  made  in  geosynchronous  orbit 
constitute  a  readily  available  data  set.  Two  preliminary  results  follow. 

Figure  2  shows  clearly  that  the  signal  power  in  the  Pc  3  frequency 
range  at  ATS-6  increased  with  Vgy,  just  as  the  hourly  maximal  amplitude  does 
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on  the  ground.  The  data  base  for  Figure  2,  and  Figure  3  below,  is  some  3000 
two-hour  intervals  of  ATS-6  magnetometer  recordings  between  June  and  Sep¬ 
tember  I971*. 

The  inserts  in  Figure  3  reproduce  the  three-dimensional  summary 
diagrams  from  Appendix  C,  showing  the  general  trend  of  the  joint  distribu¬ 
tion  of  hourly  pulsation  maximum  with  hourly  average  V  and  the  cosine  0 

5W  XB 

when  20-percent  of  the  hourly  cone  angles  of  the  IMF  were  below  0VO.  The 

XB 

large  central  diagram  in  Figure  3  summarizes  the  data  from  ATS  6  in  the  same 
format  as  the  inserts.  In  this  case,  two-hour  averages  of  Vcu  and  0WD  were 
used  and  the  vertical  scale  represents  the  probability  of  occurrence  of  Pc  3> 
based  on  the  appearance  of  definable  power  peaks  in  the  magnetic  field 
spectra. 

We  see  that  one  result  of  the  ATS  study  has  been  a  remarkable  veri¬ 
fication  of  the  joint  dependence  of  pulsation  activity  on  VCIJ  and  0  The 

same  pattern  of  pulsation  activity  rising  with  solar  wind  speed  and  cosine  of 
the  cone  angle  appears  at  synchronous  orbit  as  on  the  ground,  despite  the 
somev  1  different  means  of  processing  the  data.  The  dependence  of  6B  on  V_u 
is  much  steeper,  and  its  occurrence  much  higher  at  low  angle  than  high,  and 
the  dependence  of  <5B  on  cos  0^  is  more  pronounced,  and  the  occurrence  much 
higher  at  high  velocity  then  at  low.  There  is  a  hint  in  the  ATS  diagram  that 
Pc  3  occurrence  peaks  or  levels  off  at  c  600-700  Km/sec,  but  a  data  set 
containing  many  more  of  the  rather  unusual  speeds  above  700  Km/sec  will  be 
necessary  to  sustain  such  an  inference. 

The  entire  ATS  6  study  will  be  the  subject  of  a  report  in  prepa¬ 
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ration. 
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In  addition  to  the  studies  of  the  foregoing  paragraphs,  two  other 
investigations  have  been  initiated,  although  neither  has  advanced  enough  to 
provide  results  for  this  report.  The  first  will  attempt  to  find  correlations 
of  pulsations  at  the  ground  with  identified  surface  wave  events  recorded  at 
the  magnetopause  by  ISEE  1.  The  second  will  attempt  to  define,  develop  and 
test  a  composite  index  of  worldwide  pulsation  activity. 

SUMMARY  AND  RECOMMENDATIONS 

This  project  has,  within  its  own  activities,  demonstrated  the  control 
of  daytime,  midperiod  pulsation  activity  by  the  joint  influence  of  the  solar 
wind  velocity  and  the  IMF  direction.  In  addition,  the  project  has,  through 
the  effect  of  its  publications,  stimulated  parallel  results  by  other  investi¬ 
gators  which  have,  except  for  some  details  still  to  be  reviewed,  corroborated 
our  findings. 

Much  of  what  has  been  done  has  established  qualitative  relationships 
between  V,.^  and  cone  angle  0^B,  and  pulsation  signals  represented  in  several 
guises:  hourly  average,  hourly  maximum,  two-hourly  power  peak,  or  probability 
of  occurrence  above  background.  The  major  thrust  of  these  qualitative  results, 
aside  from  mere  proof  of  their  existence,  has  been  to  focus  attention  on 
general  signal  amplitude  rather  than  sporadic  appearance  of  visually  defined 
"wavetrain"  events.  Thus,  the  magnetosphere  must  be  thought  of  as  a  system 
responding  to  solar  influences  in  a  wholly  general  way  and  not  just  resonating 
on  occasion  at  certain  periods  or  transferring  wavetrains  to  the  ground  only 
at  special  times.  The  system  contains  resonant  responses,  and  the  propagation 
of  quasi-sinusoidal  waves  to  any  particular  point  at  any  particular  time  is 


enhanced  according  to  prevailing  parameters  of  the  system,  but  it  appears  that 
at  any  time  the  solar  wind  is  delivering  stimulus  to  the  magnetopause,  the 
effect  can  probably  be  detected  in  the  magnetosphere  at  some  level  by  "listening" 
to  signal  strength  over  most  of  the  midperiod  band. 

This  project  has  also  produced  some  quantitative  results.  Indeed, 
one  of  its  contributions  has  been  to  de-emphasize  the  event-approach  to  pul¬ 
sations  that  has  been  customary  by  operating  on  signal  level  rather  than  on 
selected  sets  of  "continuous  pulsations."  We  have  determined  that  the  dis¬ 
tinction  in  period  between  Pc  3  and  Pc  4  is  of  little,  if  any,  value  when 
correlating  signal  level  with  solar  wind  parameters,  and  that  the  distinction 
between  Pc  4  and  Pc  5  is  probably  no  better;  that  pulsation  signals  should 
be  treated  in  terms  of  transient  amplitude  (or  power)  enhancement  and  that 
correlations  should  be  expected  over  intervals  of  ten  minutes  or  longer, 
depending  on  period;  that  thresholds  of  about  250  Km/sec  and  55°  exist,  such 
that  pulsation  activity  is  negligible  when  V^w  and  0XR  are  jointly  below  and 
above  their  respective  thresholds. 

It  is  recommended  as  an  outcome  of  this  program’s  findings  and 
review  of  the  results  of  other  concurrent  investigations  that  future  effort 
move  to: 

a.  Investigate  direct  observations  of  wave  transfer  into  and 
through  the  magnetosphere, 

b.  Develop  an  index  or  measure  of  global  geomagnetic  pulsation 
act ivi ty , 

c.  Quantify  the  relationship  between  global  activity  and  solar 
wind  properties. 
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Figure  1.  Solar  wind  streamlines  in  the  magnetosheath  and 
wind-wave  speed  parameters  for  *»5°  IMF. 
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Abstract.  The  aaplitude  of  Pc  3,4  magnetic 
pulsations  at  Calgary  is  shown  to  increase  as  the 
solar  wind  velocity  increases.  The  pulsation  am¬ 
plitude  dependence  on  solar  wind  velocity  can  be 
understood  in  terms  of  the  Kelvln-Helmholtz  in¬ 
stability  at  the  magnetopause  boundary. 

Introduction 

Any  method  for  monitoring  the  properties  of 
the  solar  wind  and  magnetosphere  using  ground- 
based  rather  than  satellite  records  would  decrease 
che  cost  of  obtaining  this  much  needed  information 
about  the  space  environment.  Consequently,  cor¬ 
relations  between  solar  wind  properties  and  fea¬ 
tures  of  ground  magnetic  records  have  long  been 
sought  (cf.  Jacobs,  1970).  In  particular,  it  has 
been  suggested  that  the  amplitudes  of  Pc  3  (10-45 
sec)  and  Pc  4  (45-150  sec)  magnetic  pulsations  at 
the  earth's  surface  are  highly  correlated  with 
such  solar  wind  parameters  as  the  interplanetary 
magnetic  field  (IMF)  orientation  (Bol'shakova  and 
Troitskaya,  I968j  Greenstadt  and  Olson,  1976, 

1977;  Webb  and  Orr,  1976),  Alfven  waves  in  the 
solar  wind,  and  solar  wind  velocity  (Vinogradov 
and  Parkhomov,  1974).  Furthermore,  the  periods 
of  the  pulsations  have  been  correlated  with  solar 
wind  velocity  (Troitskaya  and  Gul'yel'ml,  1967), 
IMF  magnitude  (Troitskaya  et  al.,  1971;  Gul'yel'ml 
and  Bol'shakova,  1973;  Gul'yel'ml  et  al.,  1973; 
Russell  and  Fleming,  1976)  and  solar  wind  density 
(Grlngauz  et  al.,  1971).  Each  correlation  devel¬ 
oped  by  the  studies  cited  above  has  been  somewhat 
ambiguous  because  of  compounded  physical  effects, 
methodological  imperfections,  or  both.  In  partic¬ 
ular,  the  recent  work  by  Greenstadt  and  Olson 
(1976,  1977  referred  to  as  G  and  0)  identified 
correlations  of  IMF  direction  with  Pc  3,4  noise, 
despite  wide  scatter  of  individual  data  points. 

In  the  work  to  be  described  in  this  paper  we  used 
the  G  and  0  data  set  of  hourly  Pc  3  and  4  ampli¬ 
tudes,  plus  plasma  measurements  from  two  satel¬ 
lites  to  demonstrate  that  solar  wind  velocity  was 
also  correlated  with  the  amplitude  of  Pc  3  and 
Pc  4  activity. 

Pc  3,4  Pulsation  Amplitude  vs.  IMF  Orientation 

Greenstadt  and  Olson  examined  the  relationship 
between  the  amplitude  of  dayslde  magnetic  actlv- 
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lty  in  the  Pc  3  and  4  frequency  range  and  Byg, 
the  angle  that  the  interplanetary  magnetic  field 
makes  with  the  sun-earth  line.  They  found  that 
large  aaplitude  activity  was  associated  with  small 
6g,  l.e.  when  the  solar  wind  and  IMF  are  nearly 
aligned.  Their  results  were  consistent  with 
greenstadt 's  previous  work  (Greenstadt,  1972a 
and  b)  regarding  the  propagation  and  convection 
of  large  aaplitude  quasi-parallel  shock  oscil¬ 
lations  through  the  magnetosheath  to  the  magneto¬ 
pause. 

Using  data  from  Calgary  during  September 
through  November  1969,  G  and  0  displayed  Pc  3 
and  Pc  4  hourly  maximal  noise  amplitude  vs.  hour¬ 
ly  minimal  angle  0xb-  This  display  Illustrated 
that  larger  amplitudes  are  associated  with  smaller 
SrB  min;  however,  the  large  scatter  In  the  data 
which  contained  many  large  amplitude  events  at 
relatively  large  Bxb  min,  suggested  that  other 
parameters  may  have  a  significant  influence  on 
the  pulsation  amplitude. 

Pc  3,4  Pulsation  Amplitude  vs.  Solar  Wind  Velocity 

A  plot  of  hourly  Pc  3  activity  amplitude  vs  .so¬ 
lar  wind  velocity  using  the  G  and  0  data  is  shown 
in  Figure  1.  The  solar  wind  velocities  are  hour¬ 
ly  averages  from  Explorer  33  and  HE0S  1  [normal¬ 
ized  to  Explorer  33  values  (Moreno  and  Slgnorlnl, 
1973) )  from  the  composite  Solar  Wind  Plasma  Data 
Set  compiled  by  J.H.  King  (personal  communication 
1976).  The  median  amplitudes  are  shown  In  50 
km/ sec  bins,  and  the  amplitude  lncreaaes  with  in¬ 
creasing  velocity.  The  amplitude  In  the  500-600 
km/sec  range  is  nearly  twice  (1.8)  the  amplitude 
in  the  300-400  km/sec  range,  and  67Z  of  the  high 
amplitude  points  (>2y)  have  velocities  >.500 
km/sec. 

The  trend  for  larger  amplitude  ground  activity 
to  be  associated  with  large  velocities  Is  also 
demonstrated  in  Figure  2  illustrating  hourly  max¬ 
imum  amplitude  of  Pc  4  activity  as  a  function  of 
solar  wind  velocity.  Again,  the  median  ampli¬ 
tudes  clsarly  demonstrate  that  large  amplitudes 
were  associated  with  larger  velocities,  and  73Z 
of  the  high  amplitude  points  (^  4y)  were  associ¬ 
ated  with  velocities  >500  km/sec. 

There  are  clearly  exceptions  to  the  general 
trends  of  the  dependence  of  pulsation  amplitude 
on  8xb  "in  and  on  solar  wind  velocity.  At  times 
of  small  8xb  or  large  solar  wind  velocity,  the 
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Figure  1.  Hourly  maximum  amplitude  of  Pc  3 
waves  (10-45  sec)  measured  at  Calgary, 
Alberta  for  local  times  from  0900  to  1600  as 
a  function  of  solar  wind  velocity. 

expected  large  amplitude  pulsation  activity  Is 
not  always  observed.  However,  the  use  of  the 
minimum  value  of  egg  during  a  one  hour  Interval 
may  be  misleading, for  8gB  nay  attain  Its  minimum 
value  only  briefly  during  the  hour  li.ierval.  One 
would  not  necessarily  expect  to  see  large  ampli¬ 
tude  activity  for  hours  during  which  9gB  nin  was 
small  for  short  Intervals.  Moreover,  large  am¬ 
plitude  activity  may  not  have  been  observed  If 
the  conditions  for  resonance  in  the  magnetosphere 
were  not  attained  at  these  times,  or  our  obser¬ 
vation  point  may  not  have  been  near  the  reso¬ 
nance  region. 

Other  exceptions  to  the  general  trend  of  de¬ 
pendence  on  solar  wind  velocity  occurred  during 
times  of  small  solar  wind  velocities  when  large 
amplitude  pulsation  activity  was  observed.  For 
example.  In  Figure  2  there  are  11  points  In  this 
class.  If  the  cutoffs  are  chosen  for  amplitudes 
>4y  and  velocities  <500  km/sec.  He  examined  the 
few  large  amplitude  points  that  occurred  at  low 
velocities  to  see  if  they  might  be  occurring  at 
times  of  small  6xb  min.  Small  8gB  min  should  be 
associated  with  large  amplitudes  according  to 
Greenstadt  and  Olson,  but  these  large  amplitude 
values  were  not  well  ordered  by  8jgj  min.  How¬ 
ever,  7  of  11  of  these  hours  of  large  amplitude 
at  low  velocity  were  found  to  occur  on  the  lead¬ 
ing  edge  of  a  high  velocity  stream.  The  leading 
edge  of  a  high  velocity  stream  Is  often  associ¬ 
ated  with  the  highest  level  of  Alfve'hlc  wave 
activity  according  to  Belcher  and  Davis  (1971), 
and  the  solar  wind  wave  amplitude  may  be  related 
to  the  increased  pulsation  activity  amplitude  at 
Calgary  (Vinogradov  and  Parkhomov,  1974).  Geo¬ 
magnetic  activity  was  moderately  high  (2£Kp<4~) 
during  the  seven  hours  associated  with  the  lead¬ 
ing  edge  of  the  high  velocity  stream,  as  expec¬ 
ted  (Sawyer  and  Haurwlti,  1976).  Therefore, 


another  possible  explanation  for  these  exceptions 
is  that  they  are  generated  within  the  magneto¬ 
sphere,  not  without. 

Finally,  we  examined  the  exceptions  to  the 
general  8xb  min  trend  where  large  amplitude  activ¬ 
ity  was  associated  with  large  0XB  min.  The  pur¬ 
pose  was  to  test  whether  large  solar  wind  veloc¬ 
ity  could  account  for  these  large  amplitude 
events.  We  chose  the  class  of  Pc  4  events  with 
amplitudes  >4y  and  8xb  min  >30*.  The  cut-offs 
are  somewhat  arbitrary;  however,  other  values 
were  chosen  without  substantially  changing  the 
results.  The  average  velocity  for  the  IB  events 
In  the  large  amplitude,  large  8xb  *ln  class  Is 
544+17  km/sec.  This  Is  In  contrast  to  the  aver¬ 
age  velocity  of  449+8  km/sec  for  all  134  events 
considered  In  this  study.  In  addition  the  veloc¬ 
ity  distribution  for  the  large  amplitude,  large 
8jtl  min  class  Is  highly  skewed  towards  large 
velocities  as  opposed  to  all  cases  which  have  a 
velocity  distribution  which  is  roughly  a  gaus- 
slan  about  the  mean.  Therefore,  large  solar 
wind  velocities  may  account  for  some  of  the  scat¬ 
ter  In  the  results  of  G  and  0. 

Discussion 

He  have  made  plots  similar  to  Figures  1  and  2 
for  Pc  3  and  4  amplitudes  as  a  function  of  solar 
wind  density  and  dynamic  pressure,  as  a  function 
of  IMF  magnitude  and  as  a  function  of  geomagnetic 
activity  as  expressed  by  the  Kp  Index.  None  of 
these  parameters  correlated  with  pulsation  ampli¬ 
tude  as  clearly  as  did  the  solar  wind  velocity. 

The  correlations  that  were  apparent  were  consis¬ 
tent  with  the  known  correlation  between  the  vari¬ 
ous  parameters  and  the  solar  wind  velocity.  Thus, 
the  amplitude  of  Pc  3  and  4  magnetic  pulsations, 
at  least  as  observed  at  Calgary  during  day  time 
hours,  appears  to  be  controlled  mainly  by  two 
external  boundary  conditions,  the  direction  of 
the  IMF  and  the  solar  wind  velocity. 

The  dual  control  can  be  understood  In  a  frame¬ 
work  in  which  disturbances  in  the  magnetosheath 
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Figure  2.  Hourly  maximum  amplitude  of  Pc  4 
waves  (45-150  sec)  measured  at  Calgary, 
Alberta  for  local  times  from  0900  to  1600  as 
a  function  of  solar  wind  velocity. 
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are  amplified  at  the  magnetopause  by  the  Kelvin- 
Helmholtz  Instability  (Landau  and  Llfahitz,  1960; 
Chandrasekhar,  1961;  Southwood,  1968;  Boiler  and 
Stolov,  1970,  1973,  1974).  The  larger  the  solar 
wind  velocity, the  sore  likely  is  wave  growth  at 
the  magnetopause.  The  source  of  the  magnetosheath 
disturbances  that  are  amplified  by  this  instabil¬ 
ity  may  be  the  quasi-parallel  shock  oscillations 
described  by  Creenstadt  (1972a)  or  any  other 
source  of  waves  xa  the  dayside  magnetosheath. 

Thus,  this  mechanism  is  consistent  with  the  re¬ 
ported  control  of  Pc  3,4  pulsation  periods  by 
the  magnitude  of  the  Interplanetary  magnetic  field 
as  well  as  the  control  of  amplitudes  by  the  di¬ 
rection  of  the  IMF.  The  Kelvln-Helmholtz  insta¬ 
bility  would  simply  act  to  amplify  these  already 
existing  fluctuations.  These  waves  would  then 
couple  to  resonant  regions  inside  the  magneto¬ 
sphere  as  discussed  by  Southwood  (1973,  1974, 

1975)  and  by  Chen  and  Hasegawa  (1974a  and  h). 
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Hourly  distributions  of  the  angle  =  arccos  (X  ■  B),  where  X  is  the  solar  ecliptic  earth-sun  axis  and  B 
is  the  interplanetary  magnetic  field  (IMF),  have  been  compared  with  hourly  Pc  3  and  Pc  4  signal  activity 
at  Calgary  (or  198  hours  in  September,  October,  and  November  1%9.  Hours  whose  distributions  were 
concentrated  above  or  below  tx ,  =  $0°  were  correlated  with  a  clear  separation  of  inactive  from  active 
micropulsation  intervals.  The  correlation  was  stronger  for  the  Pc  3  than  for  the  Pc  4  band,  and  hours  with 
significant  fractions,  say  10-20%  of  their  tx a  distributions  below  30°,  corresponded  particularly  well  with 
high  amplitudes  of  micropulsation  signals.  Signals  in  the  Pc  3  period  range  on  the  ground  can  therefore  be 
used  as  a  crude  monitor  of  certain  extremes  of  IMF  orientation. 


Introduction 

For  several  years  an  international  effort  has  been  underway 
to  define  relationships  between  micropulsation  properties  at 
the  earth's  surface  and  plasma  properties  in  the  magnet¬ 
osphere  and  solar  wind  with  the  ultimate  aim  of  establishing 
reliable  diagnostics  of  exospheric  parameters  obtainable  from 
sensitive  magnetometers  on  the  ground.  A  number  of  promis¬ 
ing  correlations  have  been  discovered,  but  gaps  in  theory  and 
methodology  still  bar  the  way  to  the  desired  objective.  This 
report  describes  an  advance  in  data  analytical  methodology 
applied  to  one  particular  pair  of  correlation  parameters. 

One  relationship  long  under  investigation  has  been  that 
between  interplanetary  magnetic  field  (IMF)  orientation  and 
micropulsation  occurrence  in  the  period  ranges  10-45  s(Pc  3) 
and  45-150  s  (Pc  4).  Very  briefly,  the  results  are  that  micro¬ 
pulsation  events  are  initiated  or  terminated  when  the  IMF 
assumes  or  departs  from  certain  preferred  orientations.  Statis¬ 
tical  scatter  has  been  large,  and  the  exact  details  are  not  wholly 
confirmed  from  one  worker  to  another.  All  seem  to  agree  that 
Pc  4  events  occur  more  or  less  when  the  IMF  aligns  itself  with 
the  sun-earth  line.  Pc  3  events,  on  the  other  hand,  have  been 
found  to  occur  when  the  IMF  either  is  along  the  nominal 
spiral  angle  [ Bol'shakova  and  Troitskaya,  1968],  is  aligned  with 
the  sun-earth  line  [ Greenstadt  and  Olson.  1976;  Nourry.  1976], 
or  possesses  an  enhanced  radial  component  [Webb  and  Orr. 
19761. 

Traditionally,  investigators  have  identified  quasi-sinusoidal 
micropulsation  ‘events’  at  an  observatory  and  then  looked  for 
the  IMF's  that  correspond  to  them.  Statistically  significant 
accumulations  of  such  events  have  demonstrated  the  patterns 
cited  above,  and  a  few  striking  individual  events  have  been 
illustrated  in  the  literature  [Bol'shakova  and  Troitskaya.  1968; 
Nourry,  1976].  IMF  orientation  has  been  defined  in  terms  of 
field  vector  latitude  X,  and  longitude  separately  and  in 
combination  [Troitskaya  el  ai,  1971;  Bol'shakova  and  Troit- 
skaya.  1968;  Gul'elmi,  1974;  Greenstadt  and  Olsonr  1976; 
Nourry,  1976].  Only  rarely,  and  recently,  has  micropulsation 
activity  been  treated,  without  preselection  of  events,  by  exam¬ 
ining  broadband,  mixed  frequency  pulsation  signal  level  and 
by  applying  spectral  analysis  [Greenstadt  and  Olson.  1976; 
Webb  and  Orr,  1976;  Barker  et  ai.  1977;  Webb  el  ai.  1977; 
Arthur  et  ai,  1977]. 
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The  ongoing  study  from  which  this  report  originated  is  an 
attempt  to  lest  a  model  according  to  which  approximate  align¬ 
ment  of  B  with  the  solar  wind  permits  large-amplutide,  quasi¬ 
parallel  bow  shock  waves  to  be  delivered  from  the  subsolar 
region  of  the  shock  to  the  subsolar  region  of  the  magneto¬ 
pause,  providing  periodic  or  nearly  periodic  perturbations  at 
the  dayside  magnetospheric  boundary  as  a  possible  instigator 
of  daytime  geomagnetic  pulsations.  Arguments  leading  to  the 
model  have  been  given  elsewhere  [Greenstadt,  1972;  Greenstadt 
and  Olson.  1976]  and  are  not  central  to  this  communication. 

It  is  worth  noting  here,  however,  that  recent  results  showing 
the  global  nature  of  pulsation  enhancements  have  strongly 
reinforced  the  association  of  Pc  3-4  origin  with  solar  wind 
properties,  particularly  changes  in  IMF  direction  [Barker  et 
ai.  1977;  Webb  et  ai.  1977], 

The  initial  thrust  of  the  study  has  aimed  at  a  simple  demon¬ 
stration  of  whether  or  not  field  alignment,  measured  by  0*», 
the  angle  between  B  and  the  solar  ecliptic  X  axis  (approxi¬ 
mately  the  solar  wind  flow  line),  is  correlated  with  signal  level 
in  the  Pc  3  and  Pc  4  bands  at  Calgary  observatory.  By  signal 
level  we  mean  the  peak-to-peak  value  of  the  horizontal  field 
variation  in  filtered  Pc  3  and  Pc  4  frequency  ranges.  Correla¬ 
tion  was  sought  earlier  between  hourly  maximal  peak-to-peak 
Pc  3-4  signal  level,  bB,  and  hourly  minimal  0XB,  by  using  data 
from  a  particular  3-month  interval,  September-November 
1969.  The  results  were  favorable  but  beset  by  so  much  scatter 
that  convincing  physical  conclusions  were  impossible.  Numer¬ 
ous  potential  explanations  of  the  weak  correlations,  both 
physical  and  methodological,  were  proposed  [Greenstadt  and 
Olson,  1976;  Singer  et  ai,  1977],  The  research  segment  covered 
here  has  attempted  to  reduce  the  scatter  by  attacking  the  most 
severe  methodological  distortion  introduced  earlier,  namely, 
the  misrepresentation  of  each  hour's  IMF  readings  by  the 
hourly  minimal  0**.  We  use  the  same  SB  as  before,  which 
represent  the  hourly  maximum  signal  level  in  either  X  or  Y 
components,  whichever  was  larger,  at  Calgary. 

The  method  chosen  for  improving  the  data  was  to  replace 
each  hour's  0x*m„,  by  an  actual  distribution  of  0* a,  of  which 
there  were  nominally  44  samples,  each  hour.  The  distributions 
were  approximated  by  calculating,  for  each  hour  in  which  at 
least  three  quarters  of  the  data  were  recorded,  the  percent  of 
the  hour's  0*»  that  fell  below  every  10*  increment  in  0X(  from 
0°  to  90°.  The  remainder  of  this  paper  describes  and  discusses 
the  picture  drawn  by  these  cumulative  (i.e.,  integral)  percentile 
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distributions  of  SXB.  Histogram  averages  indicated  by  arrows.  Results  are  shown  for  hours  when  (a)  all  #*«  <  50°  and  all 
>  50°  and  (ft)  all  <  30°  and  all  «„  >  60”. 


distributions  and  the  differential  distributions  derivable  from 
them. 

Results 

Separation  of  extremes.  The  first  and  most  urgent  appli¬ 
cation  of  the  percentile  computation  was  immediately  success¬ 
ful.  Hourly  distributions  were  found  in  which  all  Sx,  were 
clustered  near  one  extreme  or  the  other  of  the  angle  scale,  and 
the  corresponding  hourly  amplitudes  were  tabulated.  The  re¬ 
sults  appear  as  histograms  in  Figure  I.  In  Figure  Is  the  two 
classes  of  hours  all  9X„  <  50”  and  all  9xt  >  50°  are  displayed 
for  Pc  3  and  Pc  4  bands.  The  9X,  <  50”  histograms  have  been 
plotted  upside  down  for  clarity. 

The  dissimilarities  in  the  distributions  of  amplitudes  in  the 
two  classes  are  clear.  Both  the  modes  and  the  averages  of  the 
two  amplitude  distributions  are  displaced,  in  the  predicted 
direction,  for  both  micropulsation  bands.  The  displacement  is 
even  more  striking  when  more  extreme  clustering  is  demanded. 
Figure  1ft  shows  the  results  for  hours  when  all  <  30°  and 
all  9xb  >  60°.  Inthe  case  of  Pc  3,  not  only  are  the  averages  and 
modes  displaced,  but  the  amplitude  distributions  don't  even 
overlap  on  the  amplitude  scale.  The  population  of  the  classes  is 
admittedly  small,  however,  since  the  selected  cases  represent 
relatively  unusual  behavior  for  the  IMF. 

The  separation  of  hours  of  enlarged  SB  from  hours  of  essen¬ 


tially  background  SB  by  extreme  clustering  of  the  directional 
distributions  of  the  IMF's  seems  to  have  defined  necessary 
conditions  under  which  Pc  3-4  were  stimulated  or  suppressed. 
We  wish  also  to  examine  the  more  difficult  problem  of  deter¬ 
mining  sufficient  conditions.  Were  enhanced  signals  produced 
only  when  the  IMF  was  markedly  aligned  with  the  X  axis  or 
was  only  a  temporary  alignment  part  of  each  hour  enough  to 
trigger  activity? 

Grouped  distributions.  Figure  2  presents  an  overview  of  the 
198  hourly  distributions  of  9xt,  divided  into  groups  according 
to  the  amplitudes  of  the  corresponding  hourly  pulsation  max¬ 
ima.  That  is,  all  hourly  distributions  of  $XB  corresponding  to 
hours  in  which  maximal  Pc  3  amplitude  SB  was  0.5  y  or  less 
were  grouped  together;  all  distributions  in  which  SB  was  be¬ 
tween  0.5  and  1.0  y  were  grouped  together,  and  so  on.  The 
range  of  SB  defining  each  group  was  0.5  y  until  SB  became  too 
high  for  0.5  -y  to  define  a  large  enough  sample  of  distributions; 
then  larger  increments  of  SB  were  used  to  sort  the  remaining 
distributions.  The  bottom  value  of  SB  for  each  group  is  printed 
in  the  figure. 

Each  small  panel  of  Figure  2  is  an  overlay  of  all  the  cumula¬ 
tive  hourly  distributions  of  9XS  that  accompanied  the  corre¬ 
sponding  SB.  The  abscissa  is  in  angular  degrees;  the  ordinate 
measures  the  percent  of  9xt  recorded  during  a  given  hour 
below  or  equal  to  a  corresponding  9X >.  Accumulations  were 
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Fig.  2.  C umulalive  percentile  distributions  of  hourly  sets  of  Bx,  arranged  acording  to  groups  of  hours  whose  maximal 
Pc  3  or  Pc  4  signal  amplitudes  iB  fell  in  selected  ranges.  The  bottom  of  each  amplitude  range  is  indicated  by  the  iB  inserted 
in  each  small  graph;  e.g.,  I.Oy  signifies  the  group  of  hours  for  which  1.0  £  iB  S  1.5  >.  'A'  's  indicate  average  of  the  group. 
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Fig  3.  Differential  percentage  distributions  of  the  group  distribution 
averages  for  the  Pc  3  panels  in  Figure  2. 


computed  every  10°.  An  average  cumulative  distribution  made 
up  of  the  average  percentile  at  each  10°  point  was  calculated 
for  each  group  and  indicated  in  the  figure  by  a  string  of  A’s. 

The  variability  of  the  hourly  distributions  for  a  given  SB  is 
readily  apparent,  but  the  tendency  for  most  8xt  to  have  oc¬ 
curred  between  30°  and  90°  when  5 B  was  smallest,  and  be¬ 
tween  0°  and  60°  when  &B  was  largest,  is  easily  discernible. 
The  variability  of  the  distributions  in  a  given  group  made  the 
group's  average  distribution  (the  A's)  a  rather  poor  indication 
of  any  clustering  tendency  among  the  curves.  Nevertheless,  the 
averages  can  be  used  for  describing  some  of  the  statistical 
properties  of  the  distributions  free  from  the  clutter  in  Figure  2. 

Behavior  of  the  averages.  Figure  3  shows  the  histograms  of 
the  average  differential  percent  distribution,  in  10°  intervals, 
for  half-gamma  groupings  of  signal  amplitude.  The  average 
distributions  of  Pc  3  show  very  clearly  the  shift  from  concen¬ 
tration  at  high  9xi.  with  virtually  no  population  at  low  8Xi.  at 
the  bottom,  to  a  wider  spread,  with  substantial  occupation  of 
the  0°-30°  sector,  at  the  lop.  The  shift  occurred  as  the  hourly 
maximal  amplitude  increased.  Comparable  histograms  for  Pc 
4,  not  shown,  reflected  the  more  erratic,  less  self-consistent 
trend  visible  in  Figure  2’s  Pc  4  panels. 

In  another  view  of  the  same  data  we  can  clump  10°  in¬ 
crements  of  9 xe  in  groups  of  three  and  examine  the  bottom-to- 
top  properties  of  the  distributions  in  Figure  3  in  three  vertical 
bands  of  30°  each.  Figure  4  shows  how  the  percentage  of  fx» 
in  each  of  three  sectors,  0°-30°,  30°-60°,  and  60°-90°,  varied 
with  micropulsation  signal  amplitudes  of  corresponding 
hours. 


The  patterns  of  Figure  4  tell  us  the  following: 

1.  Minimal  micropulsation  signal  in  both  Pc  3  and  Pc  4 
channels  was  unmistakably  associated  with  negligible  hourly 
percentages  (<  10%)  of  $xb  less  than  30°. 

2.  Maximal  micropulsation  signal  in  both  Pc  3  and  Pc  4 
channels  was  unmistakably  associated  with  substantial  per¬ 
centages  (20-30%)  of  9xb  less  than  30°. 

3.  A  substantial,  even  predominant,  percentage  of  8Xb  be¬ 
tween  30°  and  60°  was  not  necessarily  associated  with  signifi¬ 
cant  micropulsation  signal  levels.  Restated  from  the  opposite 
point  of  view,  it  was  not  necessary  that  9Xb  be  overwhelmingly 
concentrated  near  90°  in  order  to  cut  off  Pc  3  and  Pc  4  signal. 

4.  The  trend  of  low-to-high  amplitude  with  low-to-high 
percentages  of  9Xb  in  the  0°-30°  range  was  more  pronounced 
and  consistent  in  the  Pc  3  than  in  the  Pc  4  band. 

Figure  4  suggests  that  the  extent  of  IMF  orientation  near  the 
sun-earth  line  was  the  critical  factor  in  influencing  Pc  3-4 
amplitudes.  This  impression  is  reinforced  if  we  isolate  the  low- 
angle  statistics  by  imagining  vertical  bands  of  0°-l0°  and  0°- 
20°  in  Figure  3.  Figure  5 a  shows  the  percentages  of  hourly  8xt 
less  than  20°  versus  amplitude.  We  sec  that  the  more  stringent 
the  alignment  we  require,  the  more  definite  the  trends  of 
percentage  versus  SB  become  for  both  Pc  3  and  Pc  4  period 
ranges.  This  tendency  is  consistent  with  the  highly  restrictive 
angle  requirements  and  pronounced  separations  described  ear¬ 
lier  in  Figure  I . 

Individual  exceptions.  Finally,  we  turn  from  the  over¬ 
simplifications  represchlcd  by  the  averages  and  return  again  to 
the  way  in  which  cases  were  actually  distributed.  Figure  6 
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Fif .  4.  Percentages  of  average  distributions  falling  in  three  ranges  of 
6xa»  versus  IB. 
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Fig.  5.  Percentages  of  average  distributions  falling  in  the  lowest 
ranges  of  versus  bB:  (a)  SXB  =  0°-20°  and  (4)  9,,  -  0°-I0°. 


compares  low-amplitude  cases  with  high-amplitude  cases  for 
both  Pc  3  and  Pc  4.  We  are  asking  here  a  more  rigorous 
question  than  heretofore,  namely,  how  many  individual 
(hourly)  exceptions  were  there  to  the  apparent  correlation 
between  pulsation  amplitude  and  angular  distribution? 

We  separated  all  essentially  uneventful  from  all  eventful 
hours.  For  Pc  J,  hours  when  the  peak  amplitude  was  0.5-1  y 
were  designated  uneventful;  hours  when  the  peak  amplitude 
was  greater  than  2  y  were  designated  eventful  in  the  sense  that 
some  Pc  3  signal  definitely  had  occurred.  The  analogous  sepa¬ 
ration  for  Pc  4  was  for  peak  amplitudes  ofO.5-1.5  y  (unevent¬ 
ful)  and  peak  amplitudes  greater  than  3  y  (eventful).  These 
divisions  correspond  to  separation  of  signals  above  the  aver¬ 
age  peak  amplitudes  from  signals  below  the  average  peak 
amplitudes  as  determined  by  averages  over  the  whole  study 
interval  as  described  by  (ireensiadt  and  Olson  (1977). 

We  want  to  know  how  many  limes  a  small  percentage  o(8XB 
readings  fell  at  low  8X„  even  though  there  was  no  pulsation 
event  and.  conversely,  how  many  times  a  small  percentage  was 
not  accumulated  until  high  8X„  even  though  there  was  pulsa¬ 
tion  activity  during  the  hour.  In  Figure  ba  the  distributions  are 
shown  of  the  number  of  cases  in  which  the  10  and  20  percent¬ 
iles  of  the  hourly  angle  measurements  fell  in  each  10°  in¬ 
crement  of  8XB  for  the  Pc  3  band.  For  example,  in  the  bot- 
. ommost  graph  there  were  5  cases  in  which  no  more  than  10% 
of  the  B  field  directions  0XB  were  below  30°  or,  conversely,  5 
cases  in  which  at  least  90%  of  the  angle  measurements  were 
above  30°. 

The  bottom  pair  of  histograms  declares  that  not  one  case  in 
which  Pc  3  activity  was  undetectable  had  as  much  as  10%  of 
the  corresponding  8XB  below  10°  and  only  5  out  of  40  cases 
had  as  much  as  10%  of  the9*a  below  30°;  further,  every  case  in 
which  Pc  3  activity  was  pronounced  during  the  hour  had  10% 
or  more  of  its  0XB  below  70°,  and  only  4  out  of  35  cases  failed 
to  have  at  least  10%  of  thefl*«  below  50°,  while  over  half,  or  20 
out  of  35  cases,  had  10%  or  more  8XB  below  30°.  Thus  the 
trends  indicated  by  the  averages  in  Figures  4  and  5  were 
strongly  reinforced  by  the  actual  distributions  of  cases  making 
up  the  averages. 

Comparable  statements  can  be  read  from  the  other  histo- 
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grams  of  Figure  6,  which  display  the  10-  and  20-percentile 
distributions  of  uneventful  and  eventful  cases  for  both  Pc  3 
and  Pc  4  bands.  The  total  effect  of  Figure  6  is  that  the  hourly 
distributions  of  8XB  for  the  two  pulsation  conditions  corre¬ 
spond  clearly  to  different  behaviors  of  the  IMF.  The  mode 
shifted  toward  0°  for  all  but  one  of  the  eventful  distributions, 
compared  to  the  uneventful,  while  the  median  8X,  was  45°  or 
more  for  every  uneventful  distribution  but  45°  or  less  for  every 
eventful  distribution.  The  distinctions  were  appreciably  more 
pronounced  for  Pc  3  than  for  Pc  4  periods. 

Of  great  importance  to  this  study  is  what  Figure  6  does  not 
show.  There  is  no  evidence  in  the  figure  that  any  value  in  the 
range  8XB  =  50°-60°  is  essential  to  separating  active  from 
inactive  conditions,  especially  for  Pc  4.  The  uneventful  cases, 
while  deficient  in  low  8XB,  include  numerous  cases  in  which 
20%  or  more  of  the  measurements  were  below  50°,  with  the 
pulsation  amplitude  at  background.  Conversely,  there  were 
active  cases  in  which  all  but  10%  or  less  of  the  8XB  were  above 
50°.  There  is  nothing  in  this  form  of  presentation,  then,  that 
supports  the  shock  origin  model  in  terms  of  the  notion  of  a 
subsolar  point  threshold  at  8XB  *  50°.  This  is  the  nominal 
angle  below  which  the  subsolar  point  is  in  the  quasi-parallel 
portion  of  the  bow  shock  ( Greensiadi ,  1976J. 

There  is  a  visually  attractive  demarkation  in  Figure  6,  how¬ 
ever.  All  the  uneventful  hours  seem  to  have  had  a  sharp  cutoff 
at  30°,  only  a  small  fraction  of  cases  having  their  10  or  20 
percentiles  of  8XB  below  this  angle.  Correspondingly,  all  the 
eventful  hours  appear  to  have  had  a  subdistribution  below  30° 
that  peaked  in  the  10° -20°  interval.  The  striking  separations 
of  Figure  1  will  be  recalled.  In  terms  of  the  question  posed 
earlier,  a  substantial  percentage  ( 10-20%)  of  points  below  30° 
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Hg  6.  Comparison  of  histograms  of  the  distributions  of  10  and  20 
percentiles  of  0**  for  eventful  and  uneventful  hours  of  {«)  Pc  3  activity 
and  (/>)  Pc  4  activity. 
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seemed  to  be  sufficient  to  produce  Pc  3,  4,  but  a  similar 
percentage  of  9XB  below  50°  was  not  sufficient  for  the  same 
effect .  The  statistical  significance  of  such  a  separation  cannot 
be  established  without  further  data.  Physically,  9XB  =  30°  is 
the  value  below  which,  roughly,  the  ‘forward  cap,'  i.e.,  the 
section  of  the  bow  shock  corresponding  to  the  sunward  or 
positive  X  axis,  rather  than  just  the  subsolar  point,  assumes 
quasi-parallel  structure.  This  occurs  because  every  cross  sec¬ 
tion  of  the  shock  parallel  to  the  plane  containing  B  and  X  has 
its  quasi-perpendicular  section  moved  back  behind  the  line 
X  =  0.  Why  such  a  configuration  should  be  significant  for 
micropulsations  has  not  been  considered  before,  but  this  ques¬ 
tion  will  bear  further  investigation  should  the  30°  effect  be 
confirmed.  Certainly,  such  a  widespread  production  of  large- 
amplitude  noise  in  the  shock  would  be  compatible  with  the 
global  scale  of  Pc  3,  4  range  signals  reported  by  Webb  tt  at. 
[1977],  but  the  model  we  are  studying  does  not  require  it 
[Greensladl.  1972). 

Discussion 

The  hourly  distribution  data  supported  the  notion  of  0XB  as 
a  factor  correlating  with  the  amplitude  of  micropulsation  sig¬ 
nal.  The  correlation  was  clearest  for  extreme  conditions.  That 
is,  concentrations  of  SXB  above  60°  were  coupled  with  minimal 
pulsation  signal;  concentrations  below  30°  were  coupled  with 
distinctly  enhanced  pulsation  signal.  The  data  did  not,  how¬ 
ever,  indicate  a  straightforward  connection  between  9XB  and 
bB  in  the  sense  of  any  functional  relationship,  since  there  was  a 
wide  assortment  of  distributions  associated  with  any  given 
range  of  amplitude.  Thus  it  seems  certain  that  at  least  one 
factor  other  than  0XB  played  a  significant  role  in  linking  the 
solar  wind  to  Pc  phenomena.  One  such  factor,  solar  wind 
velocity,  has  already  been  shown  to  be  important  [Singer  el  al.. 
1977). 

The  evidence  presented  in  this  report  does  strongly  support 
the  result  previously  enunciated  that  signals  of  Pc  3  periods 
followed  the  same  pattern  as  signals  of  Pc  4  periods,  tending  to 
be  enhanced  when  0XB  was  low,  although  cases  of  large  ampli¬ 
tudes  also  occurred  when  the  IMP  was  near  the  nominal 
stream  angle.  The  distinction  between  Pc  3  and  Pc  4  reported 
by  some  observers  [ Bol'shakova  and  Troitskaya,  1968;  Kovner. 
1 974 J  was  not  apparent  here. 

The  result  of  this  study  brings  the  advent  of  micropulsalion 
diagnostics  of  the  solar  wind  a  little  closer,  albeit  crudely.  The 
following  rules  were  implied,  with  reasonable  confidence,  by 
these  data: 

1.  If  in  a  given  hour  no  signal  is  recorded  above  1.0  y  in 
the  Pc  3  band  or  1.5  y  in  the  Pc  4  band  peak  to  peak  on  the 
ground  (at  Calgary),  it  is  improbable  (9  cases  in  94)  that  the 
IMP  has  been  oriented  within  30°  of  the  solar  wind  flow  more 
than  12  min  out  of  the  hour. 

2.  If  in  a  given  hour  no  (Pc  3)  signal  above  1.0  y  or  no  (Pc 
4)  signal  above  1.5  y  is  recorded,  it  is  highly  improbable  (1 
case  in  94)  that  the  IMP  has  been  aligned  within  10°  of  the 
solar  wind  flow  more  than  6  min  out  of  the  hour. 

3.  If  in  a  given  hour  Pc  3  amplitude  reaches  above  2  y,  and 
Pc  4  amplitude  above  3  y,  it  is  virtually  assured  (74  cases  in  82) 
that  the  IMP  has  been  within  60°  of  the  solar  wind  Row  for  at 
least  12  min  out  of  the  hour. 


4.  If  in  a  given  ,  ,  Pc  3  amplitude  reached  above  2  y  or 

Pc  4  amplitude  above  3  y,  there  is  an  even  chance  (42  cases  in 
82 )  that  the  I M  F  has  been  aligned  within  30°  of  the  solar  wind 
flow  at  least  6  min  out  of  the  hour  and  better  than  one  chance 
in  three  (31  cases  in  82)  that  the  IMF  has  been  within  30°  of 
the  solar  wind  flow  at  least  12  min  out  of  the  hour. 

These  results  clearly  imply  that  dissection  of  the  underlying 
angular  distributions  provides  a  considerable  improvement 
over  hourly  averages  or  hourly  extremes  of  9XB  in  studies  of 
this  type. 

Rough  as  the  suggested  rules  above  may  be,  further  quan¬ 
tification,  Oner  diagnosis,  generalization  to  other  surface  ob¬ 
servation  points,  and  clues  to  theoretical  description  of  the 
micropulsation  transfer  process  seem  to  be  coming  within 
reach. 
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Abstract*  Separate  studies  using  the  same  mi- 
cropulsation  data  base  in  the  period  range  10- 
150  s  have  shown  earlier  that  signal  levels  re¬ 
corded  during  September,  October,  and  November 
1969  at  Calgary  correlated  positively  with  both 
solar  wind  alignment  of  the  IMF  and  solar  wind 
speed,  but  each  correlation  contained  enough 
scatter  to  allow  for  influence  of  the  other  fac¬ 
tor.  In  this  report,  joint  correlations  of  velo¬ 
city  and  field  direction  with  parameters  repre¬ 
senting  hourly  distributions  rather  than  minima 
of  IMF  orientation  angle  display  the  relative 
effect  of  the  two  agents  on  magnetic  pulsation 
signal  levels.  The  joint  correlations  reduce  the 
overall  scatter  and  show  that  solar  wind  speeds 
above  200-300  km/s  and  angles  between  the  IMF  and 
the  sun-earth  line  of  less  than  50°-60®  are  as¬ 
sociated  with  enlarged  magnetic  pulsation  ampli¬ 
tudes.  These  threshold  effects  tend  to  support 
both  the  bow  shock  origin  and  the  Kelvin-Helmholt z 
amplification  of  daytime  signal  transients  in  the 
Pc  3,  4  period  ranges. 

Introduction 

In  separate  examinations  of  the  same  magnetic 
pulsation  observations  recorded  at  Calgary  Obser¬ 
vatory,  Greenstadt  and  Olson  [1976,  1977]  and 
Singer  et  al.  [1977]  demonstrated  correlations 
between  interplanetary  magnetic  field  (IMF)  direc¬ 
tion  and  the  amplitude  of  Pc  3,  4  signals  (6B)  and 
between  solar  wind  speed  (V8W)  and  <5B.  The  IMF 
correlations  supported  some,  but  not  all,  of  the 
same  relationships  reported  earlier  by  Soviet  re¬ 
searchers  [Bol'shakova  and  Troitskaya,  1968; 
Plyasova-Bakunina ,  1972],  who  had  worked  with  se¬ 
lected  Pc  3  and  Pc  4  wave  train  events  rather  than 
unselected  signal  levels.  More  pointedly  for  the 
motivation  of  this  study,  the  results  of  Greenstadt 
and  Olson  lent  mild  support  to  the  prospect  that 
transient  enhancements  of  signal  in  the  Pc  3,  4 
period  range  were  related  to  the  transfer  of  large- 
amplitude  shock  pulsations  from  the  subsolar  bow 
shock  to  the  midday  magnetopause.  Studies  by 
Webb  and  Orr  [1976]  and  Webb  et  al.  [1977],  con¬ 
ducted  in  parallel  with  our  own,  found  correla¬ 
tions  of  IMF  component  changes  with  global  pul¬ 
sation  events  that  were  compatible,  although  not 
Identical,  with  the  statistical  correlations  we 
reported  earlier. 

Meanwhile,  the  velocity  correlation  of  Singer 
et  al.  supported  the  prospect  that  the  same  Pc  3, 
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4  signals  were  related  to  excitation  and  ampli¬ 
fication  of  perturbations  at  the  magnetopause  via 
the  Kelvin-Helmholt z  instability,  whose  influence 
would  depend  on  Vsv. 

Clearly,  it  is  desirable  to  determine  whether 
the  IMF  or  the  Vsw  correlation  is  the  ’true*  one 
or  whether  both  are  before  placing  confidence  In 
either.  The  angle  ©xb  *  arc  cos  (cos  Ag  cos  4>g), 
where  Ag  and  $g  are  the  solar  ecliptic  latitude 
and  longitude  of  the  IMF,  respectively,  was  used 
in  the  cited  reports  as  the  measure  of  IMF  orien¬ 
tation.  Since  in  the  solar  wind,  tan  a 
we  have  cos  ©xb  *  vsw  cos  ^B^^aw  +  k2)*/2,  where 
k  is  a  proportionality  constant,  and  one  correla¬ 
tion  could  have  been  deterministically  controlled 
by  the  other.  Also,  the  scatter  in  either  cor¬ 
relation  was  more  than  enough  to  have  incorporated 
the  effects  of  the  other  variable  if  two  indepen¬ 
dent  mechanisms  were  indeed  present,  and  it  would 
be  desirable  to  separate  the  Influence  of  both 
while  improving  the  coherence  of  each. 

In  this  report  we  show  that  ©xb  and  Vgw  appear 
to  have  operated  independently  on  6B,  each  in  a 
manner  consistent  with  the  mechanism  though  to  be 
responsible  for  its  correlation.  We  remind  the 
reader  that  the  angle  ©xb  though  to  affect  the 
stimulation  of  daytime  magnetic  pulsations  when 
the  IMF  B  is  less  than  about  50*,  unaberrated,  or 
perhaps  60®  from  the  sun-earth  line,  allowing  for 
solar  wind  aberration,  by  providing  a  quasi¬ 
parallel  bow  shock  structure  in  the  subsolar 
region  delivering  large-amplitude  magnetic  oscil¬ 
lations  directly  to  the  magnetopause.  The  desig¬ 
nation  'pulsations  of  shock  origin*  (PSO)  has  been 
suggested.  The  solar  wind  speed  is  thought  to 
promote  the  generation  of  pulsations  in  the  mag¬ 
netosphere  by  supporting  amplification  of  PSO's 
(or  any  other  perturbations)  at  the  magnetopause 
when  Vsw  is  above  a  threshold  criterion  associated 
with  the  Kelvin-Helmholt 2  instability. 

Data 

Sources.  The  data  set  described  here  was  re¬ 
corded  in  September,  October,  and  November  1969. 
Magnetic  pulsation  data  from  Calgary  were  passed 
through  digital  filters  that  separated  the  Pc  3 
(T  ■  10-45  s)  and  Pc  4  (T  ■  45-150  s)  bands,  and 
the  maximal  amplitude  in- X  or  Y  In  each  band  in 
each  daylight  hour  was  chosen  to  represent  the 
hour's  activity.  IMF  data  were  obtained  from  a 
tape  recording  of  Explorer  35  82-s  averages. 

Solar  wind  velocities  were  represented  by  hourly 
averages  from  the  composite  Solar  Wind  Plasma 
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Fig.  1.  Development  of  pulsation  versus  velocity  correlation  with  nail  relaxation 
of  unfavorable  angular  distributions.  Panels  show  hourly  maximal  amplitude  6B  of 
pulsations  versus  average  solar  wind  velocity  for  (top)  Pc  3  and  (bottom)  Pc  4. 


Data  Set  supplied  by  the  National  Space  Science 
Data  Center. 

Additional  details  of  the  sources,  processing, 
and  characteristics  of  the  data  can  be  found  in 
earlier  reports  cited  in  the  introduction. 

Method.  The  key  to  separation  of  velocity 
from  field  orientation  effects  was  an  improvement 
in  methodology  Introduced  by  Greenstadt  and  Olson 


[1978],  who  examined  the  actual  distributions  of 
0yB  within  each  hour's  worth  of  data.  In  previ¬ 
ous  papers  the  hourly  behavior  of  the  IMF  direc¬ 
tion  had  been  characterized  by  finding  the  lowest 
value  of  6xb  recorded  for  the  hour.  The  Idea  be¬ 
hind  this  selection  of  'minimal  6jb'  was  that 
transient  enhancements,  and  hence  hourly  maxima, 
of  il  could  have  been  produced  by  transient  rota- 
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Fig.  2.  Further  development  of  pulsation  versus  velocity  correlation  with  in¬ 
creasingly  favorable  angular  distributions  (same  format  as  Fig.  1).  Equations  des¬ 
cribe  least  squares  straight  lines. 
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Fig.  4.  Contrasting  correlations  of  5B  versus  cos  Byg  for  hours  of  (left)  low 
V8W  and  (right)  high  Vgw.  (Top)  Pc  3.  (Bottom)  Pc  4.  The  horizontal  scales 
refer  to  the  cosine  of  the  angle  below  which  20Z  of  each  hour's  0xb  w“»  observed. 

those  at  corresponding  Byg  in  the  low-velocity  tion  at  Byg  makes  detection  of  any  single  thresh- 

plots.  This  general  elevation  of  amplitudes  on  old  level  of  angle  rather  more  difficult  than  in 

the  right  affirms  the  importance  of  solar  wind  the  case  of  velocity,  which  tends  to  remain  near- 

velocity  in  the  mechanism  that  brought  Pc  3,  4  ly  constant  during  most  hours  and  hence  repre¬ 
signals  to  the  Calgary  sensors.  Especially  sig-  sentable  by  one  number.  It  is  not  therefore  sur- 

nificant,  however,  was  the  contrast  between  the  prising  that  Fig.  4  fails  to  exhibit  any  clear 

slight  upward  trend  of  6B  with  cos  Byg  for  cases  value  of  cos  Byg  dividing  high  6b  from  low.  There 

in  which  Vw  was  at  the  low  end  of  its  range,  at  could  conceivably  be  a  different  threshold  Byg 

left,  and  the  more  pronounced  upward  trend  when  for  each  percentile  or  for  only  a  single,  'right' 

Vsw  was  near  the  high  end  of  its  range,  at  right.  percentile  or  a  threshold  percentile  for  a  parti- 

It  thus  appears  that  elevated  solar  wind  velocity  cular  Byg,  etc • 

was  a  necessary  factor  in  enabling  the  amplitude-  However,  in  Fig.  5  we  follow  the  guidance  of 

angle  correlation  to  operate,  just  as  the  tempo-  Fig.  1  and  look  again,  in  the  same  format  as  in 

rary  rotations  of  Bsw  for  which  a  few  minutes  of  Fig.  4,  at  the  distributions  of  SB  versus  cos  Byg 

0XB  were  below  50°  or  60°  (Fig.  1)  were  necessary  for  cos  0yg  at  the  10X  Instead  of  the  201  level, 

to  expose  the  amplitude-velocity  correlation.  Presumably,  if  there  were  a  threshold,  it  would 

Were  thresholds  indicated?  The  patterns  of  have  been  most  evident  under  slight  relaxation  of 

points  in  Fig.'s  1-3  and  the  positive  horizontal  the  cutoff  condition.  There  is  Indeed  some  lndl- 

lnteraects  of  the  best  fit  lines  in  Fig.  2  indi-  cation,  in  three  of  the  four  panels  of  Fig.  5  (all 

cate  quite  clearly  that  there  was  a  velocity  below  but  that  at  upper  left),  that  high  SB  occurred 

which  no  significant  SB  could  have  been  expected.  only  to  the  right  of  an  imaginary  vertical  line 

The  most  favorable  angular  selections,  namely,  at  about  the  middle  of  each  horizontal  scale.  We 

those  that  should  have  allowed  the  'real'  velocity  consider  it  particularly  noteworthy  that  with  the 

correlation  to  appear  with  the  least  admixture  of  exception  of  a  single  point,  every  point  in  either 

unfavorable  cases,  were  in  Fig.  3  and  in  the  right  right-hand  panel  that  exceeded  all  points  in  the 

panels  of  Fig.  2.  These  plots  suggest  that  the  corresponding  left-hand  panel  occurred  when  cos 

minimal  Vgw  for  certifiable  Pc  3,  4  signals  at  @XB  0.62,  or  Byg  S  51.°.  The  dashed  lines  out- 

Calgary  was  in  the  neighborhood  of  250-350  km/s.  line  this  behavior.  It  is  thus  suggested  that 

The  patterns  of  Fig.  1  suggest,  albeit  weakly,  truly  heroic  levels  of  5B  were  attained  at  high 

that  some  Byg  below  about  60°  had  to  occur  in  velocities  only  when  some  EBW  happened  to  rotate 

order  that  significant  pulsation  amplitudes  be  within  about  50°  of  the  solar  wind  flow  direction, 

recorded.  Obviously,  the  necessity  of  represent-  Such  behavior  would  be  entirely  consistent  with 

lng  hourly  IMF  orientations  as  a  bivariate  quantl-  the  postulated  requirement  that  quasi-parallel 
ty  involving  both  Byg  and  ®  percentage  dlatrlbu-  shock  pulsations  in  the  subsolar  region  supply  the 
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Pc  3-4  excitation  mechanism.  The  subsolar  point 
of  the  bow  shock  assumes  a  quasi-parallel  struc¬ 
ture  when  0xb  £  50*. 

Summary 

In  our  view,  all  three  issues  raised  in  the 
inquiry  section  above  were  answered  In  the  affir¬ 
mative.  The  most  self-evident  correlations  were 
in  amplitude  versus  velocity,  but  the  correlations 
were  strongest  for  the  most  favorable  distribu¬ 
tions  of  angle  and  appeared  to  require  some  favor¬ 
able  angles  to  emerge  at  all.  The  correlations 
of  amplitude  versus  angle  were  more  apparent  at 
high  velocity  than  at  low,  and  amplitudes  in  gen¬ 
eral  were  elevated  at  the  highest  velocities  over 
those  at  the  lowest  velocities.  Finally,  there 
was  evidence  that  between  50°  and  60°  there  was 
a  threshold  above  which  ©xb  did  cut  off  any  am¬ 
plitude  correlation  with  velocity,  while  it  was 
abundantly  clear  that  velocities  below  about  300 
km/s  would  be  correlated  only  with  insignificant 
amplitudes  regardless  of  angle. 

Fig.  6  summarizes  the  overall  trends  in  a 
pair  of  schematic  three-dimensional  drawings  in 
which  the  various  correlations  are  represented  by 
the  best  fit  straight  lines,  at  the  20Z  level,  of 
Fig.'s  2  and  4.  Other  combinations  of  percentile 
and  velocity  grouping  would  have  produced  dif¬ 
ferent  outlines  with  the  same  general  tendencies, 
some  weaker,  some  stronger,  and  we  would  not  em¬ 


phasize  the  quantitative  content  of  these  draw¬ 
ings  until  it  is  confirmed  by  additional  data. 

Discussion 

The  angle  threshold  of  about  50*-60° ,  the 
tightened  correlation  of  6B  versus  Vw  when  many 
©XB  S  30*,  and  the  6B  versus  cos  ©xb  correlations 
were  consistent  with  the  hypothesis  that  quasi- 
parallel  shock  oscillations  in  the  subsolar  re¬ 
gion  supply  perturbations  to  the  midday  magneto¬ 
pause  that  can  be  transferred  to  the  surface  as 
pulsations  in  the  Pc  3,  4  period  bands.  At  the 
same  time  the  6B  versus  V8W  correlations,  with 
their  cutoff  of  significant  6B  at  about  200-300 
km/s,  were  consistent  with  the  amplification  of 
perturbations  at  the  magnetopause  by  the  Kelv in- 
Helmholtz  mechanism:  Boiler  and  Stolov  [1970, 
1973]  have  shown  that  increased  solar  wind  velo¬ 
cities  should  produce  increased  magnetospheric 
disturbance  according  to  a  Kelvin-Helmholtz  sta¬ 
bility  criterion.  Velocities  above  a  threshold 
determined  by  certain  MHD  parameters  of  the 
plasmas  on  either  side  of  the  magnetopause  should 
produce  wave  growth.  In  their  investigation. 
Boiler  and  Stolov  [1973]  found  positive  support 
for  the  threshold  effect  in  satellite  data,  and 
their  tables  indicated  a  minimal  solar  wind  velo¬ 
city  over  many  cases  of  about  300  km/s  or  a  bit  less. 

The  instability  criterion  they  cited  is  inde¬ 
pendent  of  frequency,  so  it  should  apply  equally 


CASES  FOB  V$w  -  300  TO  >*00  fcn/SEC  CASES  FOR  V$w  -  500  TO  600  Km/ SEC 

5B(v)  Pc  3 


Fig.  5.  Same  as  Fig.  4,  but  the  horizontal  scales  refer  to  the  cosine  of  the 
angle  below  which  10X  of  each  hour's  ©yB  was  observed.  The  horizontal  dashed 
lines  mark  the  highest  points  observed  at  low  speed,  at  left.  The  vertical  dashed 
lines  emphasize  that  at  high  speed,  at  right,  all  values  of  6B  above  the  horizontal 
line  (except  one)  occurred  for  cos  ©xg  *  0.61. 
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Fig.  6.  Diagrammatic  summary  of  joint  6B  versus  cos  SgB  and  V sw  correlations, 
represented  by  coablnlng  the  two-dimensional  straight-line  fits  of  Fig.’s  2  and  4. 


to  Pc  3  and  Pc  4  and,  for  that  gutter,  to  Pc  5 
signals,  if  we  Ignore  possible  dispersion  and 
resonance  affecting  delivery  of  the  signals  to 
the  surface.  The  data  of  the  present  study  there¬ 
fore  support  the  applicability  of  the  Kelvin- 
Helmholtz  criterion  in  still  another  way:  the 
stronger  was  the  <5B-Vhw  correlation,  the  closer 
to  one  another  were  the  two  SB  «  0  cutoffs  of  the 
Pc  3  and  Pc  4  scatter  plots.  The  differences 
between  the  crossings  of  the  SB  -  0  axis  by  the 
Pc  4  and  Pc  3  least  squares  lines  for  201  0gg  < 

60",  50",  and  30"  (Fig.  2)  and  for  50Z  OgB  <  30" 
(Fig.  3)  were  138,  148,  72,  and  27  km/s,  respec¬ 
tively.  Hence  the  data  Imply  that  when  the  angles 
were  sufficiently  favorable  to  expose  the  most  re¬ 
liable  SB-Vsw  correlations,  the  Pc  3  and  Pc  4  sig¬ 
nal  thresholds  occurred  close  to  a  common  minimal 
solar  wind  velocity,  as  the  Kelvin-Helraholtz  cri¬ 
terion  would  suggest.  In  Figure  3,  Incidentally, 
the  two  least  squares  lines  would  cross  at  - 
250  km/a,  where  SB  -  0.44  y,  a  little  below  the 
putative  instrument  sensitivity. 

One  of  the  main  sources  of  weakness  in  the  cor¬ 
relations  shown  here  was  the  presence  of  too  many 
low  SB  at  both  high  velocity  and  small  angle;  i.e., 
enhanced  signals  often  failed  to  appear  under  what 
should  have  been  the  most  favorable  circumstances. 
Although  there  was  much  improvement  over  earlier 
scatter  diagrams,  we  suspect  the  most  probable 
single  cause  of  this  ’non-reBponse’  was  the  loca¬ 
tion  of  Calgary  Observatory  under  the  plasmapause, 
where  relatively  small  changes  in  magnetospheric 
geometry  could  have  removed  the  stations  from  the 
best  signal  path.  An  expanded  data  set  soon  to  be 
studied  should  resolve  this  issue. 

For  the  present  the  data  support  the  notion 
that  PSO’s  were  identifiable  phenomena  in  the  Pc  3 
and  4  period  ranges  during  the  interval  described 
in  this  and  its  preceding  reports. 
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ABSTRACT 

Observations  have  supported  but  not  proved  the  postulate  that  the 
stimulation  of  medium-period  geomagnetic  pulsations  occurs  in  the  interaction 
between  the  solar  wind  and  the  geomagnetic  field.  Both  IMF  orientation  and 
the  solar  wind  speed  have  been  demonstrated  to  correlate  with  the  level  of 
pulsation  activity.  Although  the  correlations  are  not  yet  free  of  substantial 
variability  caused  by  other  factors,  probably  both  physical  and  methodological, 
they  provide  some  insight  into  the  shortcomings  of  existing  theoretical  for¬ 
mulations,  the  appropriate  direction  of  further  observational  analysis,  and 
the  operation  of  magnetospher Ic  processes. 
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INTRODUCTION 


In  the  following  paragraphs,  the  current  state  of  development  and  credi¬ 
bility  of  the  principal  models  connecting  the  solar  wind  to  midperiod  (T  = 
10-300  sec),  daytime,  geomagnetic  pulsation  signals,  i.e.,  Pc  3,^,5,  is  sum¬ 
marized.  The  term  signals  is  deliberately  chosen  to  emphasize  interest  in 
wave  energy  in  the  bands  of  interest,  meant  to  include,  but  not  be  limited  to, 
traditional  monochromatic,  or  quasi-periodic  pulsation  events.  Selected, 
rather  than  exhaustive,  references  are  given,  and  suggestions  for  future 
steps  are  enunciated. 


PULSATION  MODELS 

Models  of  midperiod,  geomagnetic  pulsations  can  be  divided  into  two 
broad  categories:  internal  and  external.  Internal  models  account  for  pul¬ 
sations  in  terms  of  processes  within  the  magnetosphere,  where  opportunities 
exist,  for  example,  at  the  plasmapause  boundary,  for  generating  oscillations 
through  local  instabilities.  External  models  attribute  pulsations  to  forces 
outside  the  magnetosphere  because  circumstantial  evidence  links  their  habits 
of  occurrence  to  properties  of  the  solar  wind. 

Each  model -category  in  turn  can  rest  on  either  of  two  factors,  or  both, 
in  making  its  case:  stimulation  or  control.  Oscillations  may  be  generated  in 
one  regime,  say  at  the  magnetopause,  but  the  ingredients  required  to  activate 
the  generation,  or  to  transfer  it  to  the  earth's  surface,  for  example,  a  reso¬ 
nant  amplification,  may  be  controlled  by  the  other  regime. 

Table  1  displays  the  foregoing  classification  scheme,  listing  sample 
items  in  each  subcategory.  The  crossed  arrows  at  the  center  indicate  the 
interplay  between  stimulation  and  control  that  may  exist  among  various  items. 
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Table  1.  Classes  of  Models  of  Mid-Period 
( 1 0—300  sec,  Pc3~5)  Geomagnetic 
Pul  sat  ions 


The  ensuing  remarks  will  be  devoted  exclusively  to  the  external  models,  with 
the  particular  intent  of  outlining  their  present  status  and  future  prospects. 

STATUS  OF  EXTERNAL  MODELS 

Figure  1  is  an  attempt  to  illustrate  schematically  both  the  conceptual 
logic  of  the  two  we  1 1 -developed ,  externally-oriented  models  and  the  current 
state  of  attempts  to  verify  them  by  observation.  The  vertical  stages  of  the 
drawing,  labeled  from  the  top,  Solar  Wind,  Foreshock,  Magnetosheath,  etc., 
represent  distinct  regimes  forming  the  plasma-physical  chain  linking  the 
solar  wind  with  the  earth's  surface.  There  are  others,  such  as  ionosphere 
and  polar  cusp,  omitted  because  they  do  not  play  a  distinctly  designated 
role  in  either  model.  The  various  stages  are  linked  by  lines  of  different 
types,  as  explained  in  the  figure,  and  below. 
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Surface  Waves  Signals 


—  - - --Quantitative  Model 

. . Qualitative  Model 

—  *  —  *  —  Observations  consistent  with  model  implication 

- Observations  consistent  with  model 

—  - -Observations  consistent  with  model,  but 

incomplete 


5  ■  Perturbed  Quantity 


Status  of  external  model  verification.  Symbolic  con¬ 
nections  indicate,  outside  the  lettered  columns,  the 
physical  steps  making  up  the  respective  models,  and, 
between  the  lettered  columns,  the  qualities  of  obser¬ 
vational  evidence  connecting  the  physical  steps. 
Verification  of  either  model  will  require  solid  lines 
in  the  same  pattern  as  outside  connect  ions . 
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At  far  left  and  right  of  the  diagram,  the  logical  connections  needed  to 
operate  the  respective  models  are  shown,  i.e.,  each  indicated  connection  or 
subconnection  is  a  necessary  part  of  its  model's  physical  operation  and  must 
be  separately  verified  by  observation,  preferably  at  the  same  time  as  every 
other,  to  validate  the  model.  The  central  lines  summarize  the  experimental 
situation  that  presently  attends  either  model.  In  the  following  paragraphs, 
the  two  columns  are  described  one-by-one. 


Surface  waves.  The  surface-wave,  or  Kel vin-Helmhol tz  model,  at  left, 
begins  with  fluid-like  waves  at  the  magnetopause.  The  waves  arise  because 
the  modified  solar  wind  in  the  innermost  magnetosheath  flows  along  the  rela¬ 
tively  stationary  magnetopause,  like  wind  over  water.  Small  surface  waves  are 
amplified  when  the  outside  velocity  Vuc  exceeds  a  critical  value,  and  their 
wave  properties  are  transferred  through  the  magnetosphere  to  the  surface 
subject  to  appropriate  local  and  global  propagation  laws,  including  not  in¬ 
considerable  attenuation.  Along  the  way,  however,  conditions  surrounding  one 
or  another  of  the  finite  geomagnetic  field  lines  may  provoke  resonance, 
selecting  out  a  particular  period,  or  band  of  periods,  thus  creating  an  en¬ 
larged  amplitude,  monochromatic  oscillation  that  can  appear  at  the  foot  of 
the  corresponding  line  as  a  traditional  geomagnetic  pulsation.  The  locus- 
dependent  sheath  velocity  (and  other  parameters,  like  density)  is  tied  to  the 
solar  wind  velocity  by  fluid  theory  (Spreiter  and  Alksne,  1969);  the  Kelvin- 
Helmholtz  conditions,  the  waves'  connection  with  the  magnetosphere,  and  the 
resonances  of  the  latter,  have  been  developed  in  a  large  body  of  theoretical 
papers  of  which  only  a  sampling  is  referenced  here  (Sen ,  1965;  Radoski ,  1966; 
Atkinson  and  Watanabe,  1966;  Sou t hwood ,  1968;  Chen  and  Hasegawa,  197A; 
Lanzerotti  et  al.,  I97M. 
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The  surface  wave  model  has  been  supported  experimentally,  in  parts,  but 
by  no  means  established  in  total:  the  magnetosheath's  plasma  flow  velocities 
have  been  related  to  V^y  in  reasonable  conformity  with  theory  in  the  steady 
state  (Wolfe  and  McKibbin,  1968;  Howe  and  Binsack,  1972);  the  solar  wind  veloc- 
ity  V^w  has  been  related  to  field  variability  at  the  magnetopause,  6BMp,  but 
not  necessarily  to  surface  waves  5R,  by  inference  of  VMS  and  estimation  of 
local  Kelvin-Helmhol tz  critical  velocities  at  the  magnetopause  (Boiler  and 
Stolov,  1973);  V^,  estimated  on  the  basis  of  constant  V^y,  has  been  related 
to  6Bm  just  inside  the  magnetosphere  through  estimated  Kelvin-Helmholtz  criti¬ 
cal  velocities  (Wolfe  and  Kauffman,  1975);  perturbations  of  the  magnetopause 
resembling  surface  waves  propagating  with  VMS  have  been  identified  (Aubry  et 
al . ,  1971);  theoretically-predicted  polarization  characteristics  of  waves  pro¬ 
pagating  through  the  magnetosphere  from  waves  defined  at  the  magnetopause  by 
their  electrical  perturbation  6E^p  (Southwood ,  197*0,  have  been  recorded  at 
synchronous  orbit  (Arthur  et  al . ,  1977c;  Hughes  et  al . ,  1978)  and  on  the  ground 
(Samson  et  a  I . ,  1971;  Lanzerotti  et  al , ,  1972;  Olson  and  Rostoker,  1978);  and 
at  both,  simultaneously  (Arthur  and  McPherron,  1977a);  solar  wind  velocity  has 
been  positively  correlated  with  ground  pulsations  (Sai to,  1964;  Singer  et  al . , 
1977).  It  is  important  to  mention  that  most  of  the  cases  analyzed  by  Lanzerotti 
et  aj_.  (1972)  gave  polarization  behavior  contrary  to  that  predicted  by  Southwood 
(1374),  but  consistent  with  Chen  and  Hasegawa  (1974)  (Lanzerott i  et  al . ,  1974). 

Signals.  The  signal  transfer  model,  at  the  right  of  Figure  1,  is  based 
on  the  idea  that  waves  of  tens  of  seconds  period  already  present  in  the  modi¬ 
fied  solar  wind  of  the  magnetosheath  should,  by  driving  the  magnetopause,  or 
penetrating  it,  make  their  presence  known  in  the  magnetosphere  where  continuing 
waves,  or  similar  waves,  excited  at  similar  periods,  would  appear  at  the  earth 
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as  traditional  midday  geomagnetic  pulsations.  The  notion  arose  from  observe- 

< 

tions  that,  subject  to  local  orientation  of  the  IMF,  the  bow  shock  and  magneto¬ 
sheath  can  contain,  or  consist  of,  large  amplitude  magnetic  pulsations  of  tens 
of  seconds  period  that  may  occasionally,  but  not  always,  reach  the  magneto¬ 
pause  (Greenstadt ,  1972).  Since  the  same  conditions  that  control  the  pulsation, 
now  called  quasi-parallel,  shock  also  produce  upstream  waves  (the  foreshock) 
in  the  solar  wind,  it  has  also  been  suggested  that  these  waves  may  simply  pro¬ 
pagate  to  the  earth,  regime  by  regime,  in  some  way,  left  undefined,  but  in¬ 
cluding  amplification  at  the  bow  shock  (Kovner  et  al . ,  1976).  Theory  has  pro¬ 
vided  for  transfer  of  waves  across  the  magnetopause,  although  not  at  very  high 
efficiency  (Wolfe  and  Kauffman,  1975). 


The  experimental  situation  is  comparable  to  that  of  the  surface  wave 
model,  namely,  supportive  in  parts.  The  IMF  does  control  the  location  and 
extent  of  quas i -paral lei  structure  in  the  midday  bow  shock  (Greenstadt  et  al . , 
1970),  which  becomes  highly  disturbed  (Greenstadt  et  al . ,  1977);  the  disturbed 
magnetosheath  can  exhibit  large  quasi -per iodic  waves  in  the  vicinity  of  the 
magnetopause  (Greenstadt,  1972;  Siscoe  et  al . ,  1967;  Kauffman,  1970);  the  ap¬ 
propriate  IMF  orientation,  namely  0^g  S  50®,  where  0^g  is  the  angle  between 
the  IMF  and  the  X-axis  (i.e.,  the  solar  wind  direction),  has  been  correlated 
with  the  occurrence  of  pulsation  events  and  amplitude  of  signals  at  synchro¬ 
nous  altitude  (Arthur  and  McPherron,  1977b)  and  on  the  earth's  surface 
(Troi tskaya  et  al . ,  1971;  Nour ry ,  1976;  Greenstadt  and  Olson,  1977;  Webb  and 
Orr ,  1976;  Wol fe  et  al . ,  1978);  specific  upstream  wave  events  have  been  cor¬ 
related  with  specific  pulsation  events  on  the  earth  (Bol'shakova  ana  Troitskaya, 
1968;  Nourry ,  1976). 
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Note  that  the  two  external  models  are  by  no  means  mutually  exclusive: 
the  Kel v in-Helmhol tz  model  requires  an  initial  perturbation,  and  the  signal 
model  needs  to  get  its  perturbation  across  the  magnetopause,  preferably  with 
amplification,  or  at  least  not  too  severely  attenuated.  Joint  correlations 
of  and  with  ground  signals  have  suggested  that  both  models  could  be 
operative  simultaneously  (Greenstadt  et  al . ,  1979;  Wol fe  et  al . ,  1978). 

summary,  neither  of  the  two  viable  external  models  of  geomagnetic 
pulsations  has  been  verified.  The  best  that  can  be  said  of  both  is  that  the 
correlations  between  0^,  and  6B^  leave  no  doubt  that  solar  wind  parame¬ 

ters  play  a  role  in  control  of  pulsations.  However,  direct  observations  of 
each  intermediate  process  will  be  necessary  to  demonstrate  external  stimulation 
of  pulsations.  Taken  individually,  the  signal  model  has  the  advantage  that 
no  counterexamples  have  appeared,  except  what  may  be  implicit  in  the  large 
scatter  in  some  of  the  referenced  correlations.  Much  of  the  scatter  was 
methodological  and  not  unexpected.  In  the  case  of  the  surface  wave  model, 
there  have  been  incompatible,  but  not  contradictory,  observations,  depending 
on  whether  appropriate  refinements  of  the  model  were  employed. 

EXTERNAL  CONTROL 

External  control  of  pulsations  actually  generated  by  some  internal  mecha¬ 
nism  can  be  just  as  intricate  and  difficult  to  verify  as  external  stimulation. 
One  example  would  be  the  possible  triggering  of  instability  in  the  magneto¬ 
sphere  by  particle  beams  from  the  sheath.  Quasi-paral  le.l  shock  structure  is 
accompanied  by  diffuse  beams  of  suprathermal  (several  KeV)  ions  seen  both  in 
the  sheath  and  upstream  (West  and  Buck,  1976;  Gosling  et  al.,  1978),  If 
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these  energetic  ions,  under  certain  orientations  of  the  IMF,  have  a  pathway 
into  the  morning  or  midday  magnetosphere,  they  might  provide  the  means  to 
drive  an  instability  that  would  then  produce  waves  and  associated  resonances 
in  the  Pc  3-5  period  range. 

A  second  example  would  involve  a  variable  V^.  The  correlations  between 
the  IMF  orientation  and  ground  pulsations  have  been  studied  at  a  limited  number 
of  stations.  Since  Pc  amplitudes  tend  to  rise  toward  auroral  latitude  and 
also  increase  with  magnetic  activity  and  V^,  and  since  the  auroral  latitude 
drops  with  activity,  amplitudes  at  a  fixed  station  may  indeed  increase  statis¬ 
tically  with  Vsu,  as  observation  has  indicated,  without  implying  anything 
about  the  mechanism  of  Pc  stimulation  at  the  magnetopause. 

It  may  be  assumed  that  further  investigations  of  the  external  models 
will  automatically  separate  the  factors  of  stimulation  and  control. 

WAVES  IN  THE  MAGNETOSPHERE 

Either  of  the  external  models  is  tied  to  questions  of  what  the  proper¬ 
ties  governing  waves  are  in  the  magnetosphere,  at  least  in  the  dayside.  The 
Kelvin-Helmhol tz  model,  for  examples,  has  already  raised  some  unresolved 
issues.  Figure  2a  shows  an  equatorial  cross  section  of  the  magnetopause,  with 
surface  waves  beginning  about  where  theory  predicts  they  should,  and  obser¬ 
vation  has  inferred  they  do  (Bol ter  and  Stolov .  1973;  Wolfe  and  Kauffman, 

1975) ,  i.e.,  about  40°  on  either  side  of  the  sun-earth  line.  The  surface 
waves  presumably  propagate,  and  grow,  in  the  direction  of  If  waves  seen 

to  propagate  in  the  magnetosphere  and  on  the  ground,  east  and  west  from  noon, 
were  indications  of  a  Kelvin-Helmhol tz  process  (Hughes  >?t  al ■ ,  1978;  Olson 
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and  Rostoker,  1978),  how  did  these  waves  get  from  their  locations  on  the  mag¬ 
netopause  flanks  to  the  noon  line  and  then  turn  around  again  to  produce  the 
observed  pattern? 

One  answer  may  be  implied  by  Figure  2b.  The  modified  solar  wind  may 
also  reach  critical  velocities  at  about  1*0°  from  the  subsolar  (stagnation) 
point  in  the  vertical,  noon-midnight  meridian  plane.  If  magnetic  activity 
conditions  are  such  that  the  cusp  is  northward  of  the  *»0°  point,  surface 
waves  might  arise  that  would  grow  and  propagate  both  with  Vuc  and  away  from 
the  noon  meridian  (in  and  out  of  the  plane  of  the  picture)  In  conformity  with 
observation.  A  northward  or  southward  component  of  the  IMF,  or  effects  of 
seasonal  geometry,  might  create  an  asymmetric  response  in  the  northern  and 
southern  hemispheres.  Details  of  such  possibilities  have  yet  to  be  worked 
out  theoretically,  but  it  is  clear  that  elaboration  of  them,  and  observational 
tests,  if  favorable,  will  advance  our  understanding  of  the  way  waves  propa¬ 
gate  between  the  magnetopause  and  the  ground,  in  the  east  and  west  directions, 
and  between  the  northern  and  southern  hemispheres. 

FUTURE  EFFORTS 

Theory.  Two  theoretical  efforts  are  called  for,  one  related  to  each 
external  model.  First,  the  surface  wave  model  needs  to  be  focussed  on  the 
operation  of  the  Kelvin-Helmhol tz  instability  in  the  noon  meridian  plane: 
wave  propagation,  seasonal  effects,  magnetic  activity  effects.  This  aspect 
of  the  problem  was  set  aside  because  the  theoretical  requirement  of  a  thin 
boundary  was  thought  to  apply  best  at  the  equator  (Southwood ,  197*0.  As 
Figure  2b  suggests,  this  requirement  may  not  be  violated  at  higher  latitudes 


at  noon,  at  least  some  times.  Second,  the  signal  model  needs  to  be  developed 
quantitatively  to  encourage  easier  testing  of  its  operation  and  its  credibi¬ 
lity. 

Observation.  The  measurements  required  to  determine  the  validity  of 
the  external  models  can  best  be  read  directly  from  Figure  1.  Any  theoretical 
connection  between  any  two  or  more  stages  needs  to  be  matched  by  a  correspond 
ing  solid  line  where  there  is  none  now,  which  is  almost  everywhere.  In  par¬ 
ticular,  tests  of  the  surface  wave  model  would  benefit  greatly  by  an  effort 
to  find  simultaneous  surface  and  ground  waves  which,  under  the  influence  of 

critical  V  ,  appear  and  disappear  together,  while  tests  of  the  signal  model 
MS 

should  be  designed  around  the  concept  of  detecting  the  progress  of  signals 
from  the  sheath  to  the  ground  when  the  appropriate  IMF  conditions  prevail. 

The  IMS  data  offer  a  bright  prospect  for  advancing  these  efforts. 
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ABSTRACT 
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obtained  in  September  and  October  1970,  confirm  and  extend  similar  results 
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evidence  for  a  velocity-controlled  contribution,  such  as  the  Kelvin-Helmhol tz 
instability,  to  pulsation  activity. 


«J 


Page  E.k 


INTRODUCTION 

Observations  have  linked  occurrence  and  intensity  of  medium-period  geo¬ 
magnetic  pulsations  (T  -  10  to  300  sec)  to  Interplanetary  magnetic  field  (IMF) 
orientation  and  solar  wind  speed  (V^y) .  The  motivation  for  these  observations 
Is  the  belief  that  magnetospher ic  pulsations  are  excited  by  perturbations 
brought  by  the  solar  wind  to  the  dayside  magnetopause  where  they  may  be  ampli¬ 
fied  by  the  Kelvin-Helmholtz  instability.  Data  supporting  these  notions  is  not 
easily  obtained,  so  new  cases,  augmented  statistics,  and  fresh  approaches  are 
being  actively  sought.  The  purpose  of  this  letter  is  to  show  that  earlier 
results  on  the  velocity  correlation  are  confirmed  when  the  same  analysis 
technique  is  extended  to  a  new  data  interval,  a  second  station,  and  a  third 
band  of  pulsation  periods. 

IMF  orientation  has  been  represented  by  latitude  and  longitude  angles 
of  the  IMF,  Xg,  4»B ,  in  the  solar  ecliptic  (SE)  system,  and  by  angle  9Xg 
between  the  IMF  and  V,.y  or  between  the  IMF  and  SE  X-axis  as  an  approximation 
to  the  solar  wind  flow  direction  (cos  0^  ■  cos  Ag  cos  ♦g) .  Pulsations  appear 
when  Ag  =  0*;  $  s  0",  180";  9Xg  S  50®  (Gul 'clmi ,  197*»,  and  references  therein; 
Nourry  and  Watanabe,  1973;  Nourry ,  1976;  Webb  and  Orr,  1976;  Greenstadt  and 
Olson,  1977),  and  correlated  positively  with  V^y  (Singer  et  al . ,  1977).  Most 
recently,  effort  has  focused  on  examining  0Xg  and  V^y  together  (Greenstadt 
et_  al_. ,  1979;  Wolfe  et  al . ,  1979). 

Various  techniques  for  defining  pulsation  activity  and  solar  wind  parame¬ 
ters  have  been  employed,  from  isolating  Individual  events  to  comparing  power 
spectra  and  hourly  averages.  Our  present  study  has  examined  correlations  of 
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Pc  3  and  4  at  Calgary  observatory  with  both  0^  and  Vsy,  separately  and 
jointly,  as  already  cited.  The  study  has  also  developed  the  application  of 
hourly  distributions  of  (Greenstadt  and  Olson ,  1979),  using  satellite  and 
ground  data  from  September,  October,  and  November  1969. 

The  new  data  interval  was  September,  October  1970,  which  satisfied  the 
criterion  that  recordings  were  still  available  from  Calgary,  permitting  a 
cross-check  for  replication  of  previous  results.  The  second  station  was  Leduc 
Observatory,  which  was  recording  simultaneously  slightly  farther  north 
(corrected  geomagnetic  latitude  61.2*;  L  *  **.2)  than  Calgary  (58.7*;  L  *  3-7). 
The  third  band  of  periods  was  150  to  300  seconds,  corresponding  to  Pc  5  pul¬ 
sations.  It  had  originally  been  intended  to  enlarge  the  correlations  of 
0XB  as  well  as  VJW  using  the  new  data,  but  the  natural  behavior  of  the  INF 
didn't  support  the  effort,  as  explained  below. 

DATA 

Measurements  of  horizontal  and  vertical  field  variations  made  at  Calgary 
and  Leduc  stations  were  passed  through  digital  filters  separating  Pc  3,  **, 
and  5  frequency  bands,  and  the  maximal  amplitude  6B,  peak-to-peak,  for  each 
hour  was  used  to  represent  the  level  of  occurrence  of  magnetic  pulsation 
activity  during  the  hour.  The  hourly  SB's  from  the  two  stations  were  inter- 
compared,  axis  by  axis,  for  the  entire  data  interval,  to  test  whether  activity 
at  the  two  stations  was  comparable. 

The  solar  wind  velocity  data  used  in  this  study  consisted  of  three-hour 
averages  from  a  merged  Vela  data  set,  including  Vela's  2-6,  which  Is  described 
by  King  (1 977a, b) . 


For  comparison  of  pulsation  activity  with  solar  wind  parameters,  the 
largest  6B  in  either  X  or  Y  component  each  hour  was  used.  Contributions  of 
vertical  (Z)  components  were  always  smaller  and  were  not  used.  Corresponding 


IMF  data  registered  by  Explorer  35,  in  lunar  orbit,  were  analyzed  to  find  the 
percentile  distributions  of  0XB  between  0*  and  90*  each  hour.  The  distribu¬ 
tions  were  used  to  separate  groups  of  "IMF-favorable"  hours;  i.e.,  most  0 

aD 


near  0*,  from  "IMF-unfavorable"  hours;  I.e., 


most  0XB  near  90*, 


for  various 


percentile  criteria.  Only  two  really  unfavorable  hours  were  found  in  which 


all  0VD  >  60*.  and  both  were  at  low  velocity,  so  favorable  and  unfavorable 

aD 

IMF  orientation  effects  on  the  velocity  dependence  of  6B  could  not  be  con¬ 


trasted.  Relaxation  of  the  criterion  of  unfavorable  conditions  failed  to  im¬ 


prove  these  statistics  in  any  useful  way.  Similarly,  no  dependence  of  6B  on 
0  could  be  developed  for  either  high  or  low  velocities  since  the  unfavorable 

AD 

ends  of  the  scale  had  no  populations.  Readers  interested  in  further  details 
of  instrumentation  or  data  handling  procedures  are  referred  to  Greenstadt  and 
Olson  (1976,  1977,  1979),  Singer  et  al.  (1977),  and  Greenstadt  et  al .  (1979). 


RESULTS 

Station  Correlations.  The  correlations  of  hourly  peak  amplitudes  of 
X-components  In  Pc  3  and  Pc  4  bands  at  Leduc  and  Calgary  are  displayed  in 
Figure  1.  Correlation  coefficients  were  .498  and  .79  respectively.  Similar 
or  better  results  were  obtained  for  the  other  components  and  for  Pc  5.  The 
Y-components  of  Pc  3  and  Pc  4  gave  correlation  coefficients  of  .679  and  .865 
respectively;  the  X  and  Z  components  of  Pc  5  gave  coefficients  of  .942  and 
.912  respectively.  Hourly  averages,  rather  than  peaks,  also  yielded  similar 
correlations. 

Note  that  the  slope  of  the  Pc  4  plot  in  Figure  I  is  >  1.0.  This  was 


true  also  for  the  Y  axis  of  Pc  4  and  for  both  Y  and  Z  axes  of  Pc  5*  Slopes 
of  all  component  correlations  of  Pc  3  were  <  1.0. 
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Velocity  Correlations.  We  present  the  data  in  the  same  'Wiost  favorable" 
format  used  by  Greenstadt  et  al.  (1979,  Figure  3).  Scatter  plots  of  6B  vs 
Vsw  are  displayed  for  hours  In  which  50  percent  or  more  of  the  0XB  were  within 
30*  of  the  X  axis.  Figure  2  shows  the  scatter  plots  and  least-square  lines 
for  hourly  maximal  pulsation  amplitudes  vs  hourly  average  V.  in  Pc  3  and 
Pc  4  period  ranges  at  both  Calgary  and  Leduc  Observatories.  The  dashed  lines 
are  the  linear  least  square  best  fits  found  for  the  earlier  data  from 
September-November  1969  (Greenstadt  et  al . ,  1978,  Figure  3).  The  correlations 
were  weaker  at  Calgary  in  the  1970  Interval  than  In  the  1969  interval,  although 
stronger  for  Pc  4  than  for  Pc  3,  just  as  in  the  earlier  study.  The  correla¬ 
tions  at  Leduc  were  as  strong  (Pc  4)  or  stronger  (Pc  3)  than  at  Calgary.  In 
fact,  the  present  Pc  3  correlation  at  Leduc  was  stronger  than  the  Pc  3  cor¬ 
relation  at  Calgary  for  the  1969  interval. 

Scatter  diagrams  for  Pc  5  are  shown  in  Figure  3-  Correlations  of  ampli¬ 
tude  with  velocity  were  clearly  strong  at  both  stations. 

DISCUSSION 

Both  Calgary  and  Leduc  stations  are  at  magnetic  latitudes  corres¬ 
ponding,  on  average  and  depending  on  the  signal  period, 

to  a  strong  latitude  gradient  in  pulsation  amplitude,  with  the  maximum  some¬ 
where  northward  of  both  observatories  (Samson  et  al . ,  1971;  Wertz  and 
Campbel 1 ,  1976).  Our  best  fit  lines  having  slopes  >  1  for  Pc  4  and  5  were 
consistent  with  this  latitude  effect.  The  original  choice  of  Calgary,  the 
southernmost  of  the  University  of  Alberta  chain  of  stations,  as  the  surface 
instrument  for  our  earlier  analysis,  was  based  on  a  desire  to  depart  no  more 
than  necessary  from  the  latitude  of  Borok  (53*),  the  station  used  for  similar 
analyses  by  Soviet  workers  (Greenstadt  and  01  son,  1976;  1977,  and  references 
therein).  We  also  desired  to  avoid  possible  confusion  introduced  by  highly 
geomagnet i cal ly-act ive  auroral  phenomena.  The  rather  clean  results  obtained 
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here  from  Leduc  plus  the  well-known  increase  in  Pc  3,  k  amplitudes  with  ap¬ 
proach  to  auroral  latitudes,  suggest  that  stronger  correlations  with  solar 
wind  parameters  of  the  type  we  are  studying  might  be  obtained  by  using  data 
from  stations  still  further  north  than  Leduc,  bearing  in  mind  that  the  best 
station  at  any  particular  time  might  depend  on  local  time  (Samson  et  al., 

1971)  and  overall  geomagnetic  activity.  The  data  used  in  this  report  all 
came  from  relatively  Inactive  conditions,  i.e.,  I  <  <.  19.  In  fact,  seven 

out  of  18  hours  in  Figures  2  and  3  were  designated  Q  or  QQ  (NOAA,  Solar  Geo¬ 
physical  Data).  The  position  of  Leduc,  incidentally,  is  comparable  to  the 
northernmost  observatory  utilized  by  Webb  and  Orr  (1976). 

It  is  important  to  note  that  Pc  3  correlations  of  Calgary  and  Leduc  did 
not  show  the  same,  let  alone  a  more  pronounced,  slope  effect,  particularly  in 
view  of  the  report  by  Fukunishi  and  Lanzerotti  (197k)  who  found  Pc  3  larger 
than  Pc  k  when  L  ^  k.  Also,  the  correlation  of  Pc  3  amplitudes  with  velocity 
was  much  poorer  than  those  found  for  the  earlier  data  interval  or  for  the 
other  bands  in  this  data  interval.  The  limited  range  of  Pc  3  and  the  smaller 
population  of  points  may  have  combined  to  obscure  the  velocity  dependence  for 
these  relatively  small  pulsations.  However,  these  deficiencies  and  inconsist¬ 
encies  in  the  Pc  3  results  probably  underscore  an  inherent  methodological 
difficulty  in  pulsation  analysis:  there  can  be  multiple  resonances  at  separate 
latitudes  in  the  piasmapause-plasmatrough  region  for  any  given  Pc  period  range, 
and  any  particular  resonance  will  appear  at  different  latitudes  at  different 
times,  depending  on  magnetic  activity  (Orr  &  Webb,  1975)  and  local  time  (Orr, 
1979).  Figures  1  and  3  of  Orr  (op  ci t . )  illustrate  these  problems  clearly. 

It  may  be  necessary  to  extract  a  Pc  3  "index"  from  recordings  at  several 
stations  to  obtain  a  quantifiable  measure  of  signal  activity  comparable  with 
itself  on  the  same  baseline  from  one  day  to  the  next. 
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Certainly  for  the  Pc  k  band  the  earlier  B  vs  Vgy  results  are  confirmed 
here,  as  they  have  also  been  by  Wolfe  et  al .  (1979)  who  analyzed  data  from  a 
station  at  Pittsburg,  New  Hampshire  (L  *  3.5).  Thus,  the  contribution  of  the 
velocity-dependent,  Kelvin-Helmhol tz  mechanism  to  pulsation  generation  remains 
consistent  with  measurement,  on  application  of  the  previous  analytic  technique 
to  a  new  data  interval.  This  consistency  adds  to  the  persuasiveness  of  the  Kelvin- 
Helmhol tz  model,  at  least  under  favorable  conditions  of  IMF  orientation,  and 
Is  reinforced  by  the  correlations  obtained  here  In  the  Pc  5  band,  which  en¬ 
compassed  the  periods  expected  from  this  model  (Atkinson  and  Matanabe,  1966; 
Southwood ,  1968;  Chen  and  Hasegawa,  197*0.  However,  alternative  models,  less 
developed  quantitatively  than  Kelvin-Helmhol tz ,  are  by  no  means  ruled  out. 

Russell  and  Elphic  (1979)  have  proposed  a  dayside  flux  transfer  process  as  a 
possible  source  of  pulsations  at  the  magnetopause,  and  of  course  the  perverse 
possibility  persists  that  elevated  V$w  is  so  correlated  with  reduced  latitudes 
of  the  auroral  edge  of  the  closed  part  of  the  dayside  magnetosphere  that  at 
any  fixed  station  SB  rises  with  regardless  of  pulsation  origin,  simply 
because  of  increased  sensitivity  of  the  station.  Confirmation,  and  further 
study,  of  the  dependence  of  SB  on  IMF  orientation  will  have  to  await  a  more 
extensive  and  varied  data  population  than  was  provided  by  the  September- 
October  1970  interval. 


ACKNOWLEDGMENTS 

Discussions  with  J.  C.  Samson,  A,  Wolfe,  and  R.  L.  McPherron,  and 
important  comments  by  the  reviewer,  have  been  helpful.  This  study  is  spon¬ 
sored  by  AFOSR  Contract  F  **9620-77-0-0018  (at  TRW),  the  National  Research 
Council  of  Canada  (at  the  University  of  Alberta),  and  NSF  Grant  ATM7**-23*»64 
(at  the  University  of  California,  Los  Angeles). 


Page  E.IO 


REFERENCES 

Atkinson,  6.,  and  Uatanabe,  T.,  Surface  waves  in  the  magnetospher ic  boundary 
as  a  possible  origin  of  iong  period  geomagnetic  micropulsations.  Earth 
Planet.  Sci.  Lett.,  J_,  89,  1966. 

Chen,  L.,  and  Hasegawa,  A.,  A  theory  of  long-period  magnetic  pulsations,  steady 
state  excitation  of  field  line  resonance,  J.  Geophys.  Res.,  79,  102k, 

197k. 

Fukunishi,  H. ,  and  Lanzerotti,  L.  J.,  ULF  pulsation  evidence  of  the  plasmapause. 
1.  Spectral  studies  of  Pc  3  and  Pc  k  pulsations  near  L  •  k,  J.  Geophys. 
Res.,  79,  1 k2 ,  197k. 

Greenstadt,  E.  W. ,  and  Olson,  J.  V.,  Pc  3,k  activity  and  interplanetary  field 
orientation,  J.  Geophys.  Res.,  8?_,  591 1,  1976. 

Greenstadt,  E.  W. ,  and  Olson,  J.  V.,  A  micropulsation  contribution  in  the 

Pc  3-k  range  correlated  with  IMF  radial  orientation,  J.  Geophys.  Res., 

82,  k991.  1977. 

Greenstadt,  E.  W. ,  and  Olson,  J.  V.,  Correlation  of  geomagnetic  pulsation 

signals  in  the  10-  to  15~second  period  range  with  concentration  of  IMF 
orientations  near  the  Sun-Earth  line,  J.  Geophys.  Res.,  8k,  lk93,  1979. 

Greenstadt,  E.  W.,  Singer,  H.  J.,  Russell,  C.  T.,  and  Olson,  J.  V.,  IMF 

orientation,  solar  wind  velocity,  and  Pc  3"k  signals:  A  joint  distri¬ 
bution,  J.  Geophys.  Res.,  8k^  527,  1979. 

Gul'elmi,  A.  V.,  Diagnostics  of  the  magnetosphere  and  interplanetary  medium 
by  means  of  pulsations,  Space  Sci.  Rev.,  16,  331,  197k. 


Page  E.  1 1 


King,  J.  H.,  Interplanetary  Medium  Date  Book,  NSSDC/WDC-A-R&S ,  77-04,  1977a. 

King,  J.  H.,  Interplanetary  Medium  Date  Book-Appendix,  NSSDC/WDC-A-R6S , 

77-04a,  1977b. 

Mourry,  G.  R. ,  Interplanetary  magnetic  field,  solar  wind  and  geomagnetic  micro¬ 
pulsation,  Thesis,  Unlv.  of  British  Columbia,  Department  of  Geophysics 
and  Astronomy,  1976. 

Nourry,  G.  R. ,  and  Watanabe,  T. ,  Geomagnetic  micropulsations  and  interplanetary 
magnetic  field  (Abstract),  EOS,  5^,  1179,  1973. 

Orr,  D. ,  Geomagnetic  pulsations  and  their  relationship  to  structure  within  the 
magnetosphere,  Magnetospher ic  Study  1979,  Proc.  Int'l  Workshop  on  Selected 
Topics  of  Magnetospheric  Physics,  Tokyo,  124,  Japanese  IMS  Committee,  1979. 

Orr,  D.  and  D.  C.  Webb,  Statistical  studies  of  geomagnetic  pulsations  with  periods 
between  10  and  70  sec  and  their  relationship  to  the  plasmapause  region 
Planet.  Space  Sci.,  £3,  1169,  1975. 

Russell,  C.  T.,  and  Elphic,  R.  C.,  ISEE  observations  of  flux  transfer  events 
at  the  dayside  magnetopause,  Geophys.  Res.  Lett.,  £,  33,  1979- 

Samson,  J,  C.,  Jacobs,  J.  A.,  and  Rostoker,  G.,  Latitude-dependent  characteris¬ 
tics  of  long-period  micropulsations,  J.  Geophys.  Res.,  76,  3675,  1971. 

Singer,  H.  J.,  Russell,  C.  T.,  Kivelson,  M.  G.,  Greenstadt,  E.  W. ,  and  Olson, 

J.  V.,  Evidence  for  the  control  of  Pc  3,4  magnetic  pulsations  by  the 
solar  wind  velocity,  Geophys.  Res.  Lett.,  4^,  377,  1977- 

Southwood,  D.  J.,  The  hydromagnet ic  stability  of  the  magnetospheric  boundary. 

Plan.  Space  Sci.,  J_6,  587,  1968. 

Webb,  D.  and  Orr,  D.,  Geomagnetic  pulsations  (5“50  MHz)  and  the  interplanetary 
magnetic  field,  J.  Geophys.  Res.,  8l_,  5941,  1976. 

Wertz,  R.,  and  Campbell,  W.  H.,  Integrated  power  spectra  of  geomagnetic  field 
variations  with  periods  of  .3-300  s,  J.  Geophys.  Res.,  8j_,  5131,  1976. 

Wolfe,  A.,  Lanzerotti,  L.  J.,  and  Maclennan,  C.  G.,  Dependence  of  hydromag- 
netlc  energy  spectra  on  interplanetary  parameters,  J.  Geophys.  Res., 

In  press,  1979- 


Page  E.12 


FIGURE  CAPTIONS 


Figure  1.  Comparisons  of  hourly  maximal  pulsation  amplitudes  SB  (X-axis) 
at  Leduc  and  Calgary  Observatories  for  September,  October  1970; 
Pc  3  at  left.  Pc  It  at  right.  Note  that  the  Pc  A  amplitudes 
tended  to  be  larger  at  Leduc  than  at  Calgary. 


Figure  2.  Correlations  of  hourly  SB  with  V ^  at  Calgary  (left)  and  Leduc 
(right)  for  Pc  3  (upper  panels)  and  Pc  4  (lower  panels).  Each 
point  represents  an  hour  In  which  50  percent  or  more  of  the 
0Xg's  were  less  than  30°. 

Figure  3.  Correlations  of  SB  with  for  the  Pc  5  band,  for  hours  in 
which  >.  50  percent  of  0^^  <  30®. 
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Figure  1.  Comparisons  of  hourly  maximal  pulsation  amplitudes 
6B  (X-axis)  at  Leduc  and  Calgary  Observatories  for 
September,  October  1970;  Pc  3  at  left.  Pc  4  at 
right.  Note  that  the  Pc  4  amplitudes  tended  to  be 
larger  at  Leduc  than  at  Calgary. 
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Figure  2.  Correlations  of  hourly  SB  with  V^y  at 
Calgary  (left)  and  Leduc  (right)  for 
Pc  3  (upper  panels)  and  Pc  1*  (lower 
panels).  Each  point  represents  an 
hour  In  which  50  percent  or  more  of 
the  0's  were  less  than  30°. 
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ASSOCIATION  OF  LOW-FREQUENCY  WAVES  WITH  SUPRATHERMAL  IONS 
IN  THE  UPSTREAM  SOLAR  WIND 

G.  Paschmann1,  N.  Sckopke1,  S.  J.  Bame2,  J.  R.  A  sb  ridge2, 

J.  T.  Gosling2.  C.  T.  Russell5  and  E.  W.  Greenstadt4 


Abstract.  Observations  obtained  upstream  of  the  earth’s  bow  shock 
with  the  LASL/MPI  plasma  instruments  and  the  UCLA  magneto¬ 
meters  on  ISEE-1  and  2  have  revealed  a  sinking  relationship  be¬ 
tween  the  presence  of  low-frequency  fluctuations  in  solar  wind 
density  and  field  strength  and  the  different  types  of  distribution 
functions  of  upstream  ions.  Waves  are  absent  when  the  ions  have  the 
beamlike  distribution  of  the  •‘reflected” ions.  Large-amplitude  waves 
are  present  only  in  conjunction  with  the  'diffuse”  ions,  which  are 
characterized  by  flat  energy  spectra  and  broad  angular  distributions. 
The  waves  are  largely  compressive,  showing  very  good  correlation 
between  oscillations  in  magnetic  field  strength  and  plasma  density. 

Introduction 

Ions,  apparently  of  bow  shock  origin,  with  energies  several  times 
the  solar  wind  flow  energy  are  known  to  exist  in  the  region  upstream 
of  the  earth’s  bow  shock,  generally  on  the  dawn  side  ( Asbndge  et  ai, 
1968).  It  has  been  suggested  that  the  low-frequency  (0.01  to 
0.05  Hz)  magnetic  waves  observed  in  the  same  region  might  be 
caused  by  these  ions  as  they  stream  through  the  solar  wind  plasma 
( Fairfield ,  1969;  Barnes  1970;  Greenstadt  et  ai,  1970).  A  direct  cor¬ 
relation  between  the  presence  of  the  waves  and  the  upstream  ions 
was  reported  by  Scarf  et  al.  (1970).  Ions  of  considerably  higher 
energies  have  also  been  observed  in  the  upstream  region.  Lin  et  al. 
(1974)  reported  the  nearly  continuous  presence  of  protons  with 
energies  above  —  30  ke V  on  field  lines  likely  to  be  connected  to  the 
shock.  Even  at  energies  above  -  100  keV,  bursts  of  protons  are 
quite  common  in  the  upstream  region  (West  and  Buck.  1976;  Sarris 
et  ai,  1976).  Recently  it  was  shown  ( Gosling  et  ai.  1978)  that  there 
actually  are  two  distinctly  different  populations  of  upstream  ions  at 
energies  between  —  1  and  -40  keV,  termed  "reflected”  and  "dif¬ 
fuse”  ions  by  these  authors.  The  reflected  ions  arc  characterized  by 
their  narrow  spectral  and  angular  extent,  i.e.  they  form  a  beam 
directed  outward  from  the  shock  along  interplanetary  field  lines. 
These  ions  are  thus  the  type  originally  described  by  Asbridge  et  al 
(1968).  The  diffuse  ions,  on  the  other  hand,  are  characterized  by 
rather  flat  energy  spectra  and  broad  angular  distributions.  Peak 
intensities  are  much  smaller  for  the  diffuse  ions,  although  the  total 
densities  are  nearly  equal. 

It  is  the  purpose  of  this  paper  to  show  that  the  noted  difference  in 
upstream  ion  distributions  is  associated  with  a  dramatic  difference  in 
the  occurrence  of  the  upstream  low-frequency  waves.  In  addition, 
the  paper  will  demonstrate  the  good  correlation  between  magnetic 
field  and  plasma  density  in  the  wave  events.  The  observations  were 
obtained  with  the  LASL/MPI  plasma  experiments  and  the  UCLA 
magnetometers  on  ISEE-1  and  2.  Detailed  descriptions  of  the 
instruments  and  data  analysis  techniques  may  be  found  elsewhere 
( Bame  et  al..  1 978;  Paschmann  et  ai,  1 978,  Russell,  1 978). 
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Observations 

Figure  1  shows  a  24  minute  interval  of  plasma  and  magnetic  field 
data  obtained  on  8  December  1977.  The  satellites  at  this  time  were 
located  near  18  RE  geocentric  distance  and  300°  ecliptic  longitude, 
and  moved  inward,  encountering  the  bow  shock  after  21:30  UT. 
The  separation  vector  in  GSE  coordinates  is  (51.  ’349,122)  km  with 
1SEE-2  closer  to  the  earth.  The  magnetic  field  is  shown  in  terms  of 
its  total  strength,  B.  The  plasma  data  are  presented  in  terms  of  the 
density.  NE.  of  solar  wind  electrons  and  the  density  Np,  of  supra- 
thermal  ions.  The  electron  density,  NE,  was  computed  from  the  two- 
dimensional  distribution  function  obtained  during  one  spacecraft 
revolution  (—  3  seconds).  To  avoid  contamination  by  photoelec¬ 
trons.  only  electrons  with  energies  above  approximately  10  eV  were 
included.  Since,  in  addition,  spacecraft  charging  was  neglected,  ab¬ 
solute  electron  densities  are  probably  accurate  only  within  a  factor 
of  -1.5. 

The  suprathermal  ion  densities.  NP.  were  obtained  by  integration 
over  only  those  energy  channels  clearly  above  the  solar  wind  protons 
and  alphas.  In  the  case  of  Figure  1,  the  integration  extended 
from  -  3.1  keV  to  -  40  keV  for  ISEE-1  and  from  -  2.6  keV  to 
~  10  keV  for  ISEE-2.  The  difference  in  limits  results  from  the  fact 
that  the  two  instruments  were  operated  in  different  energy  modes  at 
this  time. 

Figure  1  shows  that  initially  no  suprathennaf  ions  are  present  and 
that  the  solar  wind  is  very  quiet,  with  a  density  near  10  cm-*,  a  field 
strength  near  5y.  and  a  bulk  speed  (not  shown)  of  300  km  sec-1. 
Near  20:59  UT  a  low  flux  of  reflected  ions  first  appears,  persisting 
until  -  21:02  UT.  (The  seemingly  higher  densities  seen  on  ISEE-2 
are  an  artefact  caused  by  the  slightly  lower  integration  limit  on 
ISEE-2,  as  was  noted  above).  The  appearance  of  these  reflected  ions 
has  no  noticeable  effect  on  NE  and  B.  At  21 : 02  UT  there  is  a  sudden 
onset  of  magnetic  field  and  plasma  density  oscillations,  coincident 
with  an  equally  sudden  increase  in  the  suprathermal  ion  density.  An 
inspection  of  the  ion  distribution  functions  (not  shown)  indicates 
that  these  ions  are  of  the  diffuse  type.  The  wave  amplitudes  are  very 
high,  with  both  ANE/NE  and  AB/B  being  of  order  one.  The  trans¬ 
verse  oscillations  of  the  field  (not  shown)  are  even  larger  than  the 
compressional.  The  electron  temperature  stayed  nearly  constant  at 
18  x  10*  K  during  the  entire  interval.  Note  the  good  correlation 
between  the  NE  and  B  oscillations  and  between  the  ISEE-1  and  -2 
records.  The  wave  periods,  in  the  spacecraft  frame,  are  seen  to  be 
slightly  less  than  one  minute  but  with  considerable  variability. 
Detailed  examination  of  the  magnetic  records  shows  that  the  phases 
observed  by  ISEE-2  are  delayed  by  0.3 -1.0  sec  relative  to  ISEE-1. 

A  detailed  investigation  of  the  waves  is  beyond  the  scope  of  this 
paper.  It  requires  that  the  propagation  velocity  vector,  frequency, 
and  polarization  be  determined  in  the  plasma  frame  of  reference. 
While  the  plasma  velocity  is  accurately  known,  the  propagation 
direction,  determined  by  minimum  variance  analysis  of  the  magnetic 
data  appears  to  be  highly  variable,  and  sometimes  nearly  at  right 
angles  to  the  spacecraft  separation  vector.  Thus  the  observed  phase 
delay  between  the  two  satellites  cannot  be  converted  to  a  reliable 
wave  propagation  speed  by  use  of  a  simple  model  of  plane  phase 
fronts  sweeping  past  the  spacecraft.  The  wave  polarization  is  also 
found  to  be  complicated  and  variable. 

In  the  data  just  discussed  there  was  some  suggestion  that  the 
presence  of  the  waves  was  correlated  with  the  diffuse  ions,  but  not 
with  the  reflected  ions.  The  brief  duration  and  lower  densities  of  the 
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UT  20.55  se  03  03  It  IS  21.18 


UT  20.  SS  58  03  o?  It  IS  21.19 

Fig.  1.  Upstream  wave  event  on  8  December  1*177.  The  top  three 
panels  show  the  field  strength  B  (in  gamma),  the  solar  wind  electron 
density  Nt  (in  cm  ')  and  the  suprathermal  ion  densities  Np  (in  cm  ") 
as  observed  on  ISEE-1.  The  lower  three  panels  show  B.  Nl:  and  Np 
from  1SEE-2.  Spacing  for  N,  and  NP  is  3  seconds.  B  is  sampled 
every  64  msec.  The  occurence  of  reflected  and  diffuse  ion  popula¬ 
tions  are  identified  by  the  bars  above  the  Np  plot.  Spacecraft  posi¬ 
tion  was  near  18.0  Rj  geocentric  distance.  306°  longitude  and  23.5° 
latitude  in  ecliptic  (GSE)  coordinates. 

reflected  ions  could,  however,  have  accounted  for  the  effect.  The 
next  example  will  demonstrate  that  the  situation  is  indeed  strikingly 
different  for  the  two  types  of  ions  f  igure  2  shows  ISEF-2  electron 
and  suprathermal  ion  densities  and  magnetic  field  strength  for  one 
hour  on  Id  November  77  (An  ion  spectrogram  for  this  interval  is 
shown  in  Figure  2  of  Gosling  et  al .  1078)  At  the  beginning  the 
satellite  was  located  at  17.8  R,.  and  321  longitude,  moving  inward 
and  crossing  the  bowshock  near  20:50  I  IT.  The  ion  density  has  been 
computed  for  energies  between  -  2.5  and  -  10  ke\  Figure  2  shows 
that  the  low-frequency  waves  do  not  appear  until  about  10  01  UF, 


UT  19  3n  40  50  oo  to  20  19.30 


Fig.  2.  Magnetic  field  strength,  electron  and  suprathermal  ion  densi¬ 
ties  from  ISEE-2  for  a  wave  event  on  Id  November  1977.  Spacing  of 
plasma  and  magnetic  field  measurements  is  12  and  0,25  seconds, 
respectively.  Spacecraft  position  initially  was  17.8  Rj.  geocentric 
distance,  320.5°  longitude  and  24  2”  latitude  in  GSE  coordinates. 


I  SEE  l  ORBIT  8  OUTBOceC  8  NOV  1977 


UT  19.45  S5  06  15  25  36  2o.4S 


UT  23.00  10  2o  30  40  50  24.00 


Fig  3.  Electron  and  suprathermal  ion  densities  and  magnetic  field 
strength  from  1SEE-!  for  three  hours  on  8  November  1177.  Note 
the  large  time  gap  between  the  second  and  third  intervals  shown. 
The  spacecraft  crosses  the  bow  shock  near  18:53  UT,  at  16.1  Rj 
geocentric  distance.  299°  longitude  and  22.2°  latitude  in  GSE  co¬ 
ordinates.  At  the  end  of  the  interval  shown  the  spacecraft  has 
proceeded  to  a  distance  of  Id, 4  R,  Again,  occurences  of  reflected 
and  diffuse  ion  populations  are  indicated  by  the  bars  above  the  Np 
plot  In  addition.  Iherc  is  a  period  (between  ld:20  and  20:00  UT) 
where  ion  distributions  are  characterized  by  iniermediate  widths  in 
energy  and  angle,  indicated  by  the  bar  denoted  1M 

although  suprathermal  ions  of  substantial  density  are  observed 
throughout  most  of  this  one-hour  interval  Before  wave  onset, 
however,  the  ions  were  of  the  reflected  type,  whereas  those  observed 
afterwards  were  of  the  diffuse  type.  Note  the  presence  of  high- 
frequency  (£0.5  Hz)  fluctuations  in  field  magnitude  in  conjunction 
with  reflected  ions  for  densities  exceeding  0. 1  cm  ’. 

Further  evidence  for  the  intimate  relationship  between  the 
presence  of  the  waves  and  the  nature  of  the  ion  distribution  function 
is  provided  by  Figure  3.  which  shows  a  total  of  three  hours  of 
ISEE-I  plasma  and  magnetic  field  data  on  8  November  1977.  An 
ion  spectrogram  for  the  same  data  is  shown  in  Figure  1  of  Gosling  el 
a I  (1978).  Near  18:53  UTthe  bowshock  is  crossed,  as  evidenced  by 
the  large  changes  tn  electron  density  and  average  field  strength. 
Note  the  large  time  jump  between  the  second  and  third  sections  of 
the  figure.  The  previously  noted  absence  of  waves  in  conjunction 
with  the  presence  of  upstream  ions  of  the  reflected  type  is  clearly 
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l  SEE  1  ORBIT  15  IKBOLfO  27  NOV  1977 


UT  02:39  41  40  45  47  49  02.51 

Fig.  4.  Field  strength,  electron  and  suprathermal  ion  densities  from 
ISEE-1  on  27  November  1977,  showing  abrupt  termination  of  low- 
frequency  waves  and  diffuse  ions  -90  seconds  before  bow-shock 
crossing  (at  02:49  UT).  At  the  beginning  of  the  interval,  the  space¬ 
craft  was  located  at  -  t4.6  Rt.  geocentric  distance.  320°  longitude 
and  23. 4C  latitude  in  GSE  coordinates. 

evident  during  the  long  interval  between  -  23: 10  UT  and  --  23:57 
UT.  and  during  two  brief  bursts  near  18:55  and  19:02  UT.  In 
contrast,  large  amplitude  density  and  magnetic  field  fluctuations  are 
seen  between  -  20:02  and  -  20:34  UT.  when  ions  of  the  diffuse 
type  are  present.  In  addition  to  these  features  which  already 
emerged  from  the  previous  examples.  Figure  3  contains  a  new 
feature.  Between  19  19  UT  and  20 : 00  UT  the  upstream  ions  have  a 
distribution  which  is  broader  in  energy  and  angle  than  normally 
observed  for  the  reflected  ions,  but  also  considerably  narrower  than 
is  typical  for  the  diffuse  ions.  (See  Figure  1  of  Gosling  et  al.,  1978.) 
During  this  interval  of  "intermediate”  ion  distributions,  significant 
fluctuations  in  Nt  and  B  are  observed,  although  at  much  reduced 
amplitudes.  Note  that  the  density  of  suprathermal  ions  is  nearly  the 
same  (-0.1cm  ’)  under  all  conditions.  To  check  whether  the 
orientation  of  the  IMF  was  the  controlling  factor  on  the  appearance 
of  these  upstream  particles  and  waves  we  calculated  the  angle 
between  the  IMF  and  the  local  expected  shock  normal.  When  there 
were  no  upstream  waves  the  angle  was  58°  or  greater.  When 
reflected  particles  were  present,  it  averaged  about  50 J.  When  inter¬ 
mediate  particles  were  present,  it  averaged  about  40°  and  when 
diffuse  particles  were  present  it  averaged  about  30°. 

Finally  Figure  4  shows  an  example  of  high-resolution  data  where 
the  low-frequency  waves  (in  Nj  and  B)  and  the  suprathermal  ions 
(of  the  diffuse  type)  terminate  just  shortly  (-  90  seconds)  before  the 
bow  shock  is  crossed  ( -  02 : 49  UT).  Again,  as  in  the  first  example, 
the  field  and  density  fluctuations  are  in  phase.  The  dominant  period 
is  -  20  seconds.  At  the  time  when  the  low  frequency  waves  dis¬ 
appear,  the  high-resolution  magnetic  field  data  show  an  onset  of 
high-frequency  waves,  extending  to  the  shock.  The  latter  appear  to 
be  the  whistler  waves  studied  by  Fairfield  (1974).  The  noted  change 
in  character  of  the  waves  is  associated  with  a  change  in  the  direction 
of  the  magnetic  field.  When  the  low  frequency  waves  are  present, 
the  field  makes  an  angle  of  about  30°  to  the  expected  shock  normal. 
When  the  high  frequency  waves  are  present  and  the  diffuse  ions  are 
absent  the  field  makes  an  angle  of  60"  to  the  shock  normal.  In  both 
cases  the  field  direction  is  such  that  connection  to  the  shock  seems 
inevitable  Therefore  it  is  unlikely  that  the  sudden  change  in  particle 
and  field  behavior  could  have  resulted  from  disconnection  from  the 
bow  shock 

Discussion 

In  the  previous  section  we  have  demonstrated  that  the  presence  or 
absence  of  low-frequency  waves  in  the  upstream  solar  wind  is  inti¬ 
mately  related  to  the  type  of  suprathermal  ion  distribution.  Large 


amplitude  waves  were  exclusively  observed  in  conjunction  with  the 
broad  distributions  of  the  “diffuse”  ion  population,  and  not  in  con¬ 
junction  with  the  narrow-beamed  “reflected”  population.  This  fact 
bears  directly  on  the  generation  mechanism  of  the  waves.  Barnes 
( 1 970)  has  argued  that  the  waves  are  magnetoso nic  waves  caused  by 
a  cyclotron  resonance  instability  of  a  fast  ion  beam  travelling  up¬ 
stream  through  the  solar  wind  Fie  predicts  wave  amplitudes  which 
are  proportional  to  beam  density  times  its  speed  (in  the  solar  wind 
frame).  Flowever,  our  observations  indicate  that  fast  ion  beams  of 
substantial  density  can  propagate  through  the  solar  wind  for  ex¬ 
tended  periods  (30  minutes  and  longer)  without  generating  the 
waves.  Quite  obviously,  if  Barnes*  mechanism  is  correct,  there  must 
be  additional  criteria  for  instability  to  occur.  Nevertheless,  once  the 
waves  are  generated,  they  would  strongly  pitch-angle  scatter  the  ions 
(Barnes,  1970),  resulting  in  a  more  isotropic  angular  distribution.  It 
is  therefore  tempting  to  think  of  the  diffuse  ions  as  the  end  result  of 
the  interaction  between  the  original  beam  and  the  waves.  The  fact 
that  the  two  populations  are  almost  never  observed  simultaneously, 
yet  have  approximately  the  same  densities,  fits  this  picture  well.  A 
natural  consequence  would  also  be  a  correlation  between  wave 
amplitude  and  width  of  the  angular  distribution,  in  agreement  with 
the  “intemediate”  case  noted  in  the  discussion  of  Figure  3.  How¬ 
ever,  although  the  broad  angular  distributions  of  the  diffuse  ions 
might  be  explained  this  way.  their  flat  energy  spectra  extending  to 
>40  keV  are  not.  Perhaps  the  electrostatic  fluctuations  reported  by 
Scarf  et  al..  (1970)  to  be  associated  with  the  low-frequency  wave 
events  can  cause  the  energization.  Of  course,  there  might  be  no 
causal  relationship  between  the  two  ion  populations  at  all.  Indeed, 
the  correlation  of  the  appearance  of  diffuse  upstream  ions  and  waves 
with  the  orientation  of  the  field  with  respect  to  the  local  shock 
normal  suggests  that  all  upstream  ions  are  accelerated  at  the  shock, 
and  that  their  different  properties  are  associated  with  the  nature  of 
the  shock  structure  at  their  point  of  origin.  In  this  context  it  should 
also  be  pointed  out  that  ions  similar  to  the  diffuse  upstream  type  are 
observed  downstream  of  the  shock,  showing  a  similar  correlation 
with  fluctuations  in  plasma  properties  ( Asbhdge  etal .  1978). 

No  firm  conclusion  about  the  origin  of  the  different  population  of 
upstream  ions  and  the  mechanism  for  generating  the  low-frequency 
waves  appears  possible  at  this  lime.  Future  work,  for  example  a 
detailed  comparison  of  the  momentum  and  energy  fluxes  of  the  two 
upstream  ion  populations,  should  help  to  clarify  wether  a  direct 
causal  relationship  exists.  Studies  of  the  geometry  of  the  intersection 
of  field  lines  with  the  bow-  shock  and  estimations  of  conditions  at  the 
projected  source  point  of  the  ions  are  also  necessary. 

Acknowledgments.  The  authors  wish  to  thank  M.  Halbig,  J.  Papa- 
mastorakis.  E.  Sckopke.  and  E.  Tech  for  programming  support  and 
data  management.  Helpful  discussions  with  B.U.O.  Sonnerup  are 
gratefully  acknowledged.  Max-Planck-Institut  portions  of  this  work 
were  supported  by  the  Bundesministerium  fur  Forschung  und  Tech¬ 
nology.  Los  Alamos  portions  of  this  work  were  conducted  under  the 
auspices  of  the  U.S.  Department  of  Energy  with  NASA  support 
under  S- 50864  A.  Work  at  UCLA  was  supported  under  NASA  con¬ 
tract  NAS  5-20064  and  at  TRW  under  NASA  contract  NASW- 
3087. 

References 

Asbridge.  J.  R..  S.  J.  Bame  and  I.  B.  Strong,  Outward  flow  of  pro¬ 
tons  from  the  earths  bow  shock,  J.  Geophys.  Res..  73.  5777-5782. 
1968. 

Asbridge,  J.  R.,  S.  J.  Bame.  J.  T.  Gosling,  G.  Paschmann,  and 
N.  Schopke.  Energetic  plasma  ions  within  the  earth's  magneto¬ 
sheath,  Geophys.  Res.  i.ett..  5. 953-955, 1978 . 

Bame.  S.  J..  J.  R.  Asbridge.  H.  E.  Felthauser,  J.  P.  Glore.  G  Pasch¬ 
mann.  P.  Hemmerich.  K.  Lehmann  and  H  Rosenbauer,  ISEE-1 
and  ISEE-2  fast  plasma  experiment  and  the  ISEE  1  solar  wind 


212 


Paschmann  et  al.:  Association  of  Low-Frequency  Waves 


Page  F.J» 


\ 


experiment ,  IEEE  Trans.  Geosci.  Electronics,  GE-16,  216-220, 
1978. 

Barnes,  A.,  Theory  of  generation  of  bow -shock-associated  waves  in 
the  upstream  interplanetary  medium;  Cosmic  Elec.,  1,  90-114, 
1970. 

Fairfield,  D.  H.,  Bow  shock  associated  waves  observed  in  the  far  up¬ 
stream  interplanetary  medium,  J.  Geophys.  Res.,  74,  3541-3553, 
1969. 

Fairfield,  D.  H.,  Whistler  waves  observed  upstream  from  collision¬ 
less  shocks.  J.  Geophys.  Res..  79.  1368-1378,  1974. 

Gosling,  J.  T.,  J.  R.  Asbridge,  S.  J.  Bame,  G.  Paschmann  and 
N.  Sckopkc,  Observations  of  two  distinct  populations  of  bow  shock 
ions  in  the  upstream  solar  wind,  Geophys.  Res.  Lett.,  5,  957-960, 
1978. 

Greenstadt.  E.  W.,  !.  M.  Green,  G.  T.  lnouye,  D.  S.  Colburn,  J.  H. 
Binsack  and  E.  F.  Lyon.  Dual  satellite  observations  of  earths  bow 
shock,  Cosmic  Elec.,  J,  279-296,  1970. 

Lin,  R.  P.,  C.  I.  Meng  and  K.  A.  Anderson,  30-  to  100-keV  protons 
upstream  from  the  earth’s  bow  shock,  J.  Geophys.  Res.,  79,  489  to 
498,  1974. 


Paschmann,  H.,  N.  Sckopke,  G.  Haerendel.  J.  Papamastorakis,  S.  I. 
Bame,  J.  R.  Asbridge,  J.  T.  Gosling,  E.  W.  Hones,  Jr.  and  E.  R. 
Tech,  1SEE  plasma  observations  near  the  subsolar  magnetopause, 
Space  Sci.  Rev.,  in  press,  1978 . 

Russell,  C.  T  ,  The  1SEE  1  and  2  fluxgate  magnetometers,  IEEE 
Trans.  Geosci.  Electronics,  GE-16,  239-242,  1978. 

Sarris,  E.  T.,  S.  M.  Krimigis,  and  T.  P.  Armstrong,  Observations  of 
magnetospheric  bursts  of  high-energy  protons  and  electrons  at 
-  35  Re  with  IMP-7,  J.  Geophys  Res.,  81,  2341-2355,  1976. 

Scarf.  F.  L..  R.  W.  Fredricks,  L.  A.  Franck,  C.  T.  Russell,  P.  J.  Cole¬ 
man  and  M.  Neugebauer,  Direct  correlations  of  large  amplitude 
waves  with  suprathermal  protons  in  the  upstream  solar  wind, 
J.  Geophys.  Res.,  75.  7316-7322,  1970. 

West,  H.  1.  and  R.  M.  Buck,  Observations  of  >  100  keV  protons  in 
the  earth’s  magnetosheath,  J.  Geophys.  Res.,  81,  569-584,  1976 . 


(Received  November  16,  1978; 
revised  January  2,  1979; 
accepted  January  23,  1979.) 


i 


j 


31219-6011 -RU-00 


MAGNETIC  FIELD  ORIENTATION  AND  SUPRATHERMAL  ION 
STREAMS  IN  THE  EARTH'S  FORESHOCK 

by 

E.  W.  Greenstadt\  C.  T.  Russel  1^,  and  M.  Hoppe^ 


Space  Sciences  Department,  TRW  Defense  £  Space 
Systems  Group,  Redondo  Beach,  California  90278 

2 

Institute  of  Geophysics  6  Planetary  Physics, 
University  of  California,  Los  Angeles,  Los  Angeles 
California  90024 


May  1979 


Space  Sciences  Department 

TRW  Defense  6  Space  Systems  Group 
One  Space  Park 

Redondo  Beach,  California  90278 


Page  G.2 


MAGNETIC  FIELD  ORIENTATION  AND  SUPRATHERMAL  ION 
STREAMS  IN  THE  EARTH'S  FORESHOCK 
by 
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M.  Hoppe,  and  G.  Paschmann 

ABSTRACT 

The  different  populations  of  backstreaming  ions  found  earlier  with  ISEE 
instruments  outside  the  earth's  bow  shock  are  organized  by  the  orientation 
of  the  IMF  to  shock  in  the  plane  containing  and  the  solar  ecliptic  X-axis 
(the  B-X  plane).  The  associations  of  ion  classes  with  field  geometry  and 
shock  structure  are  shown  by  means  of  computer -drawn  three-dimensional 
sketches.  Particles  of  several  KeV,  or  about  ten  times  the  solar  wind  stream¬ 
ing  energy,  left  the  shock  as  "reflected"  beams  in  a  guiding  center  direction 
determined  by  Bsw  and  the  interplanetary  drift  velocity.  Their  area  of  de¬ 
parture  was  the  part  of  the  shock  where  the  transition  from  quasi -perpendicular 
to  quasi-parallel  structure  began.  Slightly  farther  toward  the  quasi-parallel 
portion  of  the  shock,  particles  of  lower  energy  escaped  upstream,  widening 
the  ion  spectra,  exciting  small  magnetic  waves,  and  forming  a  class  of  "inter¬ 
mediate"  ions.  Finally,  in  direct  connection  with  the  quas i -paral lei  structure, 
ions  leaving  the  shock  formed  the  "diffuse"  group,  containing  a  wide  distribu¬ 
tion  of  energies  including  high  values  well  beyond  any  recorded  in  the  other 
two  classes.  Sources  of  the  various  ion  groups  are  postulated. 
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INTRODUCTION 

One  of  the  striking  early  successes  of  the  ISEE  1,2  mission  has  been 
the  revelation  that  return  ions,  that  is,  ions  streaming  into  the  solar  wind 
away  from  the  bow  shock,  are  comprised  of  separate  classes  having  distinguish 
able  energy  and  directional  distributions  (Gosling  et^  aj_. ,  1978).  Further, 
each  of  these  classes  of  return  Ions,  which  are  essential ly  al 1  protons,  is 
associated  with  its  own  level  of  upstream  magnetic  wave  activity  (Paschmann 
et  al_. ,  1979). 

The  general  formation  and  variable  location  of  the  wave  and  particle 
foreshocks  have  been  known  to  be  determined  by  the  relationship  of  the  inter¬ 
planetary  magnetic  field  (the  IMF,  or  B<.y)  to  the  bow  shock.  For  an  account 
of  the  pre-ISEE  status  of  upstream  wave  and  particle  observations,  the  reader 
is  referred  to  Greenstadt  (1976)  and  Diodato  et^  aj_.  (1976).  It  was  no  sur¬ 
prise,  therefore,  that  preliminary  calculations  indicated  that  the  classes  of 
ions  were  themselves  associated  with  distinguishable  local  f ield-to-shock 

orientations,  as  represented  by  cone  angle  0_  between  BP1,  and  the  local  nor- 

Bn  bw 

mal  (  Paschmann  et^  aj_. ,  1979;  Gosling  £t  a_k  ,  1979),  and  by  the  more  sophist! 
cated  representations  of  Bonifazi  et^  a_l_.  (1979).  The  purpose  of  this  report 
is  to  supplement  the  rudimentary  calculations  of  Paschmann  et^  aj^  and  Gosling 
£]_  wi th  a  graphic  description  of  the  geometry  of  shock  and  field,  to  show 
that  the  geometry  appears  to  govern  the  particle  and  wave  classes,  and  to 
clarify  the  relationships  between  them  and  shock  structure.  Where  the  com¬ 
panion  paper  by  Bonifazi  et^  a_l_.  assumes  the  geometric  relationship  and 
parametrizes  it  in  a  general  way,  we  attempt  here  to  demonstrate  the  frame¬ 
work  with  specific  particle  pathways. 


Page  G.^ 


This  report  presents  a  series  of  schematic  figures  depicting  the  IMF- 
shock  geometry  for  selected  instants  during  the  first  two  hours  of  ISEE 
particle  data  of  8  November  1977,  published  in  Figure  3  of  Paschmann  et  al . , 
(1979).  Two  examples  were  selected  for  each  of  the  four  identified  types  of 
ion  distribution:  reflected,  diffuse,  intermediate,  and  none,  according  to 
the  nomenclature  of  Paschmann  et^al.  Similar  figures  corresponding  to  the 
various  distributions  in  the  other  figures  of  Paschmann  et  al .  have  been 
examined;  these  corroborated  the  samples  shown  here. 


The  basis  for  the  diagrams  used  herein  is  the  motion  of  charged  par¬ 
ticles  leaving  the  shock  with  guiding  center  velocities  along  B^y,  and  subject 
to  the  electric  field  E<.y  -  -m/qCV-  x  B  ).  The  guiding  center  velocity 


(1) 


of  a  proton  can  be  represented  as 

~gc  "  Po  VSW  ?SW  +  VSW  Sm  °XB 

/\  /\  /\ 

where  unit  drift  vector  D  =  E...  x  Vrll;  p  V„,.  is  the  component  of  V  along 

~SW  -SW  o  SW  gc 

B...  written  as  a  multiple  p  of  the  solar  wind  speed,  and  0vn  is  the  angle 

bW  * M 


XB 


between  X  (i.e.,  -V,.  )  and  B^y.  For  development  of  the  reflected  proton  tra¬ 
jectory,  see  Benson  et^  a_l_.  (1975).  The  result  is  that  ions  returned  to  the 
solar  wind  move  in  a  plane  containing  and  V^y  or,  if  we  neglect  aberration, 
B^y  and  the  solar-ecliptic  X-axis,  i.e.,  the  B-X  plane.  For  observations  at 
a  given  point  on  the  trajectory  of  ISEE,  we  are  therefore  interested  in  the 
particular  B-X  plane  through  the  spacecraft  location  and  in  its  curve  of 
intersection  with  the  bow  shock.  The  shock  was  based  on  the  symmetrized 
model 


Y2  +  Z2  -  .331  [ (X  -  75.25) 2  -  3686] 


(2) 


described  by  Greenstadt  et^  aj_.  (1975),  and  its  location  approximated  in  each 
case  by  scaling  the  subsolar  point  with  the  solar  wind  velocity  Vsw,  density  N, 

and  field  magnitude  measured  by  (SEE  in  the  expression 

R  -  101  (1  +  1.1  — ^NV-J176  (3) 

o  ^M2  SW 

where  Rq  is  the  distance  to  the  subsolar  point  of  the  bow  shock  and  M  is  the 
magnetosonic  Mach  number  (Formisano  et_  ak  ,  1973).  The  constant  101  was  chosen 

to  match  the  nominal  model,  and  the  y  in  the  formula  of  Formisano  et  al  was 

taken  as  5/3-  Values  of  V$w  and  N  were  close  to  275  Km/s  and  13  cm~^  for 
all  the  cases  discussed  here. 

DATA 

The  IMF  Bsw  and  its  three  solar-ecliptic  frame  components  that  accompanied 
the  varieties  of  upstream  ions  defined  by  Gosling  et^  aj^  (1978)  and  Paschmann 
et_  aJL  (1979)  for  8  November  1977  are  plotted  in  Figure  1.  The  shaded  columns 
labeled  R,  IM,  and  D  duplicate  the  designations  given  by  Paschmann  et  al . 

(1979,  Figure  3),  and  signify  reflected,  intermediate,  and  diffuse  ions, 
respectively.  The  qualities  of  these  classes  were  described  in  the  cited 
reports  and  can  be  seen  best  in  Figure  1  of  Gosling  £t  £j_.  (1978).  Some 
characteristics  of  the  classes  will  be  recalled  in  later  paragraphs.  We  as¬ 
sume  that  in  the  interval  1855  to  201*5  UT  the  location  of  the  shock  did  not 
change  radically  after  it  had  been  last  seen  at  1853,  since,  as  already  noted, 
both  ^SW  anc^  ^  remained  essentially  constant  during  the  interval,  as  did  B^w 
also. 

We  see  in  Figure  1  that  each  change  in  the  type  of  particle  observed 
was  accompanied  by  some  shift  or  drift  in  the  average  value  of  one  or  more 
components  of  the  IMF,  although  the  opposite  was  not  necessarily  true;  that  is, 
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not  every  shift  in  B  produced  an  alteration  in  particle  class.  It  is  best 
to  keep  in  mind,  however,  that  while  the  bars  define  valid  generalization  about 
the  classes  of  ions  measured,  there  were  moderate  variations  in  the  actual 
spectra  during  the  tenure  of  each  class,  requiring  somewhat  flexible  defi¬ 
nitions,  and  exhibiting  occasional,  but  rare,  transitory  exchanges  of  category. 

The  lower  case  letters  a  to  h  in  Figure  1  designate  the  eight  times  at 
which  the  field/particle  geometry  are  displayed  and  analyzed  in  this  report. 
They  were  chosen  to  be  typical  of  each  ion  category.  The  number  just  below 
each  letter  gives  the  angle  0gn  between  the  IMF  and  the  local  normal  of  the 
model  shock  of  equations  (2)  and  (3)  at  that  time. 

PRESENTATION 

Figure  2  illustrates  the  field  geometry  in  a  reflected  particle  beam 
and  also  defines  the  format  in  which  the  Information  of  this  report  will  be 
presented.  At  left  is  a  three-dimensional  sketch  of  the  bow  shock  showing 
the  orientation  of  the  B-X  plane  through  the  observation  point  and  its  curve 
of  intersection,  C,  with  the  shock.  The  shock  is  indicated  in  the  figure 
by  two  outlines  intersecting  at  the  subsolar  point,  one  in  the  ecliptic  (X-Y) 
plane,  the  other  in  the  vertical  noon-midnight  (X-Z)  plane.  A  portion  of 
the  X-Y  plane  is  shown  to  contrast  with  the  tilted  B-X  plane  above  it.  The 
point  of  observation  at  the  time  to  which  this  particular  sketch  applies  is 
denoted  by  an  "0"  (actually  a  small  square  computer  symbol  "q")  situated  on 
curve  C,  west  of  the  X-Z  plane. 
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The  sketch  at  right  was  made  by  simply  replotting  the  section  of  C 
containing  the  observation  point  on  an  enlarged  scale,  so  it  is  seen,  as  if 
with  a  "zoom"  lens,  from  the  same  viewpoint  as  at  left.  The  enlarged  figure 
shows  the  directions  of  B  and  the  guiding-center  motion  of  reflected  par- 
tides,  P,  through  0,  in  the  B-X  plane,  according  to  equation  (1).  The 
direction  P  was  obtained  by  assuming  that  return  protons  passed  outward 
through  the  observation  point  with  velocity  components  (anti-)  parallel  to 
B,  equivalent  to  an  observed  energy  of  3  keV.  The  dashed  line  represents  a 
foreshock  boundary  defined  as  the  forwardmost  path  possible  for  3  keV  par¬ 
ticles  leaving  the  shock  along  B. 

DIAGRAMS 

Beam  cases.  The  diagram  shown  in  Figure  2  corresponds  to  case  a,  at 
1855  UT,  in  Figure  1.  The  reflected  beam  appeared  from  the  solar  and  eastern 
quadrants,  principally  the  latter  (Gosling  et  a  1 ■  ,  1978).  The  beam  energy 
ranged  from  about  I  keV  to  over  3  keV  (op  c i t . ) ,  but  an  arrow  describing 
1  keV  particles  would  actually  pass  inside  the  nominal  shock,  so  protons  of 
that  energy  should  not  have  escaped  upstream.  It  seems  likely  that  the  ex¬ 
tremes  of  the  recorded  energy  band  corresponded  to  protons  of  nonzero  pitch 
angle  whose  parallel  components  occupied  a  narrower  energy  range.  That  is, 
we  imagine  escaping  protons  having  had  a  range  of  energies  divided  between 
parallel  and  perpendicular  motion  with  respect  to  B.  Some  particles  at 
either  extreme  of  the  range  of  parallel  energies  would  have  had  nonzero  pitch 
angles  and  have  been  detected  with  total  energies  a  little  above  the  maximum 
or  below  the  minimum  that  applied  to  the  guiding  centers.  The  value  3  keV 
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was  therefore  chosen  as  an  approximate  upper  limit  for  guiding-center  contri¬ 
butions  to  the  measured  beam.  Protons  of  less  energy  would  have  followed 
paths  closer  to  tangency  and  would  have  originated  on  the  shock,  sunward 
from  the  foot  of  the  3  KeV  arrow. 

A  second  case  of  a  particle  beam,  case  b  at  1902  in  Figure  1,  is  shown 
in  Figure  3-  At  that  time  the  B-X  plane,  left  of  Figure  3,  had  a  greater 
tilt  to  the  ecliptic  than  in  case  a,  but  the  relationship  of  B,  P,  and  C 
were  substantially  the  same,  as  seen  in  the  detail  at  right.  In  this  sketch 
the  path  of  2  KeV  protons  has  also  been  depicted,  together  with  its  correspond¬ 
ing  foreshock  boundary.  The  peak  energy  of  the  beam  was  probably  close  to 
this  value. 

Null  cases.  Two  cases  from  intervals  in  which  no  return  particles  were 
recorded  are  drawn  in  Figure  4.  The  cases  are  c  and  h ,  at  I9H  and  2040 
(Figure  l).  Note  that  the  B-X  planes  in  these  cases,  and  the  distances  of 
0  from  C  and  from  the  line  of  symmetry  of  C,  were  quite  different  even  though 
the  null  observation  was  the  same  at  both  times.  The  details  in  Figure  *4 
show  that  in  both  cases  the  3  KeV  paths  did  not  Intersect  the  shock  so  that 
protons  of  even  the  highest  parallel  energies  recorded  in  the  beams  earlier 
could  not  have  reached  the  point  of  measurement.  None  did.  The  dashed  tangent 
lines  representing  the  foreshock  boundary  for  3  KeV  particles  demonstrate  that 
the  observations  were  outside  that  boundary  at  those  times. 

Diffuse  cases.  Two  cases,  f  and  g,  of  diffuse  protons  at  2010  and  2030 
are  conceptualized  in  Figure  5.  Once  again,  the  B-X  planes  and  the  global 
observation  geometry  associated  with  them  were  vastly  different  from  one 
another.  The  local  field  geometry,  however,  which  depended  also  on  the  X- 
component  of  B,  produced  similar  particle  situations  in  the  two  cases.  In 
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case  f,  paths  of  1,  3.  and  10  keV  particles  are  shown,  together  with  the 
1  keV  foreshock  boundary;  the  diffuse  spectrum  ranged  from  about  .6  keV  to 
beyond  the  10  keV  limit  of  the  Los  Alamos  particle  record  (Gosling  et  al . , 
1978).  We  picked  a  1  keV  lower  limit  as  a  round  number,  allowing  for  an 
unknown  contribution  from  non-zero  pitch  angles.  Higher  energy  particles 
would  simply  have  streamed  along  lines  closer  to  8  because  of  proportionately 
smaller  contributions  from  the  drift  velocity.  The  1  keV  foreshock  boundary 
would  presumably  have  been  near  the  closest  tangent  pathway;  higher  energy 
protons  would  have  had  larger  foreshocks. 

In  the  detail  of  case  g,  the  viewing  angle  made  It  impractical  to  dis¬ 
play  more  than  the  1  keV  particle  path,  which  is  already  close  to  B.  Vir¬ 
tually  the  entire  upstream  region  constituted  the  1  keV  foreshock. 

Intermediate  cases.  Figure  6  exhibits  two  cases,  d  and  e,  at  1923 
and  1950,  from  the  section  designated  intermediate  by  Paschman  et^  aj_.  (1979). 
The  IMF  changed  little  during  the  intermediate  section  and  the  two  images  are 
substantially  the  same:  the  particles  ranged  from  a  few  hundred  eV  to  about 
3  keV,  and  the  details  of  Figure  6  show  that  the  observation  points  were  on 
or  just  inside  the  1  keV  foreshock  boundaries. 

DISCUSSION 

Reflected  beams  and  mixed  beams,  During  the  interval  described  here, 
return  particle  beams  contained  energies  from  about  1  KeV  to  3_i<  KeV  and 
entered  the  detector  from  the  eastern  and,  to  a  lesser  extent,  the  solar  quad 
rants.  "Intermediate"  particles  extended  f rom  a  few  hundred  eV  to  about  the 
same  upper  limit  as  the  beams,  3-'*  KeV,  and  entered  the  detector  from  the 
eastern,  and  in  lesser  intensities,  the  antisolar  and  solar  quadrants. 
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The  diagrams  suggest  an  explanation:  the  shock  stopped,  trapped,  and 

circulated  (through  its  "foot")  a  subset  of  the  solar  wind  protons  which  were 

accelerated  by  the  interplanetary  electric  field  to  energies  about  10  times 

the  solar  wind  streaming  energy.  These  protons  slid  westward  along  the  shock 

surface  until,  in  the  region  where  the  angle  between  the  IMF  and  the  surface 

expanded  to  about  40°,  they  were  released,  escaping  upstream  along  B  with 

drift  D.  (In  terms  of  the  conventional  complementary  angle  of  the  shock 

normal,  0  =  50°. )  The  escaping  particles  were  not  released  exactly  along 

B,  but  with  a  moderate  distribution  of  nonzero  pitch  angles.  Some  of  these 

particles,  or  others  trapped  about  where  the  first  were  released,  escaped  a 

little  farther  to  the  west,  up  to  where  the  IMF  made  an  additional  10°  or  so 

with  the  surface  (0„  =  ^O0) .  In  this  extended  region,  protons  of  lower 
Bn 

energies,  i.e.,  lower  guiding  center  velocities,  and  greater  pitch  angles 
could  also  escape.  Thus  the  intermediate  group  appear  to  have  been  a  beam 
that  had  spread  to  lower  energies  and  wider  angles  of  entry  at  the  detector 
by  virtue  of  relaxed  containment  geometry  at  the  shock,  or,  more  accurately, 
a  mixture  of  beams  from  a  larger  area  of  the  shock  surface. 

Diffuse  ions.  The  diffuse  ions  ranged  from  a  few  hundred  eV  to  energies 
beyond  10  keV.  They  entered  the  detector  from  all  quadrants  but  with  greatest 
intensity  from  the  east.  The  diffuse  ions  therefore  were  a  separable,  if 
not  independent  phenomenon  from  the  collimated,  reflected  particles  at  the 
forward  edge  of  the  proton  foreshock.  Strong  geometrical  support  of  this 
view  is  that  the  diffuse  ion  population  included  particles  of  high  energies 
well  above  those  found  in  the  beams.  Such  particles  could  easily  have  escaped 
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the  shock  and  reached  the  satellite  with  the  beams,  if  they  originated  in  the 
same  place,  but  none  was  observed. 

The  angles,  or  range  of  angles,  **0°  £  0^  $  50°  that  appeared  to  sepa¬ 
rate  the  source  of  beam  ions  from  that  of  diffuse  ions  are  recognized  as  the 
same  angles  that  separate  quas i -perpend icular  from  quasi-parallel  shock  struc¬ 
ture.  Although  the  magnetic  profile  of  the  November  8  shock  crossing  at 
1852:45,  which  has  been  displayed  elsewhere  (Russell  and  Greenstadt,  1978, 
Figures  15,  20),  did  not  show  the  typical  quasi-perpendicular  step,  6  a  19 
was  very  high  in  this  case  and  such  a  value  of  8  has  typically  been  associated 
with  rapid,  large  amplitude  fluctuations  in  the  shock,  even  at  large,  nearly 
perpendicular  0^  (Formi  sano  et^  aK  ,  1975).  We  therefore  describe  the  cros¬ 
sing,  with  its  connected  reflection  beam,  as  high-8,  quas i -perpend icular  at 

eBn  *  50°- 

True  quasi-parallel  structure  has  been  tied  closely  to  long  period  up¬ 
stream  waves  (Greenstadt,  1 976) ,  which  were  not  observed  after  the  November  8 

crossing,  but  were  observed  later  in  connection  with  smaller  0_  and  diffuse 

Bn 

ions  (Paschmann  et^  a]_.  ,  1979).  We  therefore  associate  the  diffuse  ions  with 
quasi-parallel  structure,  although  the  latter  structure  could  not  be  directly 
observed  in  these  cases,  the  shock  having  already  been  left  behind  the  space¬ 
craft  . 

Ions  and  waves.  The  linking  of  diffuse  ions,  quasi-parallel  structure, 
and  upstream  waves  is  supported  by  a  simple  calculation  which  at  the  same 
time  relates  the  waves  to  the  low,  rather  than  the  high,  end  of  the  back- 
streaming  ion  spectrum.  The  parameter  p  used  in  former  studies  to  character¬ 
ize  the  foreshock  boundary  (Greenstadt,  1972;  Diodato  et  alt,  1976)  can  be 
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found  in  terms  of  the  guiding  center  energy  of  an  observed  return  proton: 

1  /O 

p  ■  (£r/£’Sy  -  sin  2  0Xg)  '  +  cos  0^.  At  1923  (case  d,  Figure  6),  p  had 

the  typical  values  listed  in  Table  1.  In  that  case,  the  1  keV  foreshock 
boundary  was  tangent  close  to  the  subsolar  point  of  curve  C,  where  p  s  1.6 
according  to  Diodato  et_  a1_.  (1976).  This  implies  that  the  paths  of  particles 
a  little  below  1  keV  were  coincident  with  the  boundary  of  the  wave  foreshock. 
Diodato  et^  £]_.  found  in  general  that  p  averaged  about  2  for  the  shock  overall, 
a  little  less  near  the  subsolar  point,  a  little  more  toward  the  flanks. 
Calculations  similar  to  the  above  were  made  for  the  other  cases,  with  similar 
results.  If  we  note  that  the  intermediate  cases,  such  as  d,  corresponded  to 
the  appearance  of  small  wave  activity  (Gosling  et^  aj_. ,  1979,  Figure  10),  we 
may  infer  that  it  was  the  1  keV,  or  lower,  contributions  to  the  backstreaming 
ion  populations  that  were  closely  associated  with  the  familiar  upstream  wave 
activity  outside  the  quas i -paral lei  shock.  These  contributions  were  marginally 
released  upstream  in  cases  of  intermediate  geometry  and  were  a  major  part  of 
the  diffuse  spectrum,  where  the  waves  were  more  pronounced. 

We  envision  three  models  as  potential  explanations  of  the  diffuse  ions: 

f.  Origin  in  quasi-parallel  structure.  Since  the  q-parallel  phenomena, 
including  the  diffuse  ions,  reign  somewhere  behind  the  foreshock  boundary,  the 
solar  wind  has  already  passed  through  beams  of  reflected  ions  and  electrons 
before  reaching  the  q-parallel  region.  These  beams,  made  up  of  electrons  and 
protons  of  diverse  energies,  are  not  necessarily  co-located,  implying  the 
possible  existence  of  charge  separation  potentials.  Transit  through  the  sun¬ 
ward  foreshock  beams  may  therefore  in  itself  contribute  to  some  scatter  and 
acceleration  of  solar  wind  ions  even  before  the  large  q-parallel  waves  are 
encountered.  Then,  typical  large  amplitude  pulsations  and  upstream  waves  that 
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characterize  the  q-parallel  shock  could,  through  Fermi -mechan i sm  or  multiple 
acceleration  by  plasma  wave  bursts  in  wave  gradients,  energize  a  portion  of 
the  solar  wind  ions  to  tens  of  KeV;  the  pulsations  themselves,  and  their 
associated  upstream  waves,  would  scatter  the  particles  into  a  broad  pitch 
angle  distribution  as  they  drift  both  upstream  and  downstream.  Certainly, 
there  is  ample  evidence  that  solar  wind  ion  distributions  <_re  modified  in  the 
presence  of  the  diffuse  particles  and  their  associated  magnetic  waves  (Bonifazi 
et  aK  ,  1979). 

2.  Origins  in  quasi -parai lei  action  on  q-perpendicular ,  suprathermal 
ions.  Some  of  the  population  of  suprathermal  ions  up  to  10  E ^  created  in 
the  q-perpendicular  shock  and  sojourning  in  its  foot,  of  which  the  remainder 
escape  and  appear  upstream  as  reflected  beams,  could  slip  further  along  the 
shock,  or  be  convected  downstream  into  the  q-parallel  structure  where  they 
could  acquire  additional  energy  through  multiple  accelerations  and  then  be 
scattered  upstream  and  downstream,  essentially  from  the  disturbed  magneto¬ 
sheath  behind  the  parallel  shock. 

3.  Origin  in  the  magnetosphere.  Through  some  magnetic  pathway,  ions 
of  tens  of  keV  in  the  magnetosphere,  or  at  the  magnetopause,  could  find  their 
way  into  the  magnetosheath  where,  in  the  region  subtended  by  the  q-parallel 
shock,  their  guiding  centers  could  then  be  conveyed  along  the  average  field 

direction  in  the  sheath  (together  with  a  drift)  and  into  the  upstream  solar 
wind.  As  a  result  of  strong  scattering  by  q-parallel  sheath  pulsations  they 
would  appear  as  a  population  of  diffuse  suprathermal  ions  both  in  the  sheath 
and  in  the  foreshock. 

Additional  cases  involving  carefully  selected  spacecraft  locations  and 
analysis  combining  particle  distributions  with  global  and  local  field  geometry 
will  test  these  hypotheses. 
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CONCLUSIONS 

1.  The  quasi -perpendicular  section  of  the  bow  shock  accelerated  solar 
wind  ions  up  to  about  10  E ^  and  released  them  into  the  upstream  region  about 
where  8gn  *  50°,  thus  allowing  ions  with  highest  guiding  center  velocity  to 
escape  the  shock. 

2.  Shock-accelerated  ions  appeared  as  a  broader  beam  of  lower  energies 

and  wider  pitch  angles  from  the  section  of  the  shock  where  40  S  0O  <  50°. 

Bn 

3.  Diffuse  ions  were  associated  with  the  quasi -paral lei  section  of  the 
shock  and  included  particles  of  energies  >  10  E ^  not  accelerated  in,  or  more 
certainly  not  released  from,  the  quasi-perpendicular  section. 

4.  If  we  assume  the  concurrent  existence  of  all  specified  upstream 
components  observed  only  one-by-one  in  the  data,  the  results  are  summarized 
by  Figure  7  where,  in  a  typical  B-X  plane  the  reflected  beam  (R) ,  the  mixed 
beam,  or  intermediate  ions  ( I H) ,  and  the  diffuse  ions  ( D )  are  indicated  by 
different  shadings.  Question  marks  designate  areas  where  particle  character¬ 
istics  have  not  yet  been  recorded,  For  example,  the  amount  and  location,  if 
any,  of  penetration  by  the  boundary  beam  R  of  the  foreshock  into  the  diffuse 

ions  is  unknown.  The  field  lines  are  shown  as  wavy  in  the  foreshock  outside 
the  quasi-parallel  region.  The  diffuse  ions  are  known  to  penetrate  the 
shock  (Gosling  et_  aj_. ,  1978),  as  depicted  by  their  shading.  It  must  be 
remembered  that  every  plane  parallel  to  the  one  illustrated  would  presumably 
exhibit  the  same  conf igurat ions ,  forming  together  a  three-dimensional  form 
of  which  the  sketch  would  be  a  single  cross  section.  Finally,  note  that 
this  empirical  sketch  is  valid  strictly  where  all  the  actual  points  of  mea- 
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surements  occurred  near  the  nominal  shock  outline.  The  conceptual  extrapo¬ 
lation  to  greater  distances  and  other  locations  doesn't  take  into  account 
the  more  subtle  distinctions  that  might  emerge  from  careful  consideration 
of  the  differences  between  the  three  proposed  explanations  of  the  preceding 
DISCUSSION  section  or  of  Bonifazi  et^  aJL  (1979).  We  believe  new  observations 
are  too  close  to  achievement  to  justify  more  detailed  speculations  at  this 
time,  and  we  hope  to  determine  before  long  whether  return  ion  scattering  is 
proportional  to  distance  from  the  nominal  shock,  from  the  foreshock  boundary, 
or  from  the  beam  tangent  point. 
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Figure  1 


Figure  2 


Figure  3 


Figure  k 


Figure  5 


Figure  C 


FIGURE  CAPTIONS 

The  magnitude  and  components  of  the  interplanetary  magnetic  field 
(IMF)  during  the  upstream  ion  observations  of  Paschmann  et  al. 

( 1 979 *  Figure  3)-  Shaded  sections  cover  Intervals  during  which 
the  three  classes  of  return  particle  spectra,  reflected,  inter¬ 
mediate,  and  diffuse,  were  measured.  The  lower  case  letters  mark 
times  selected  for  the  examples  in  the  ensuing  figures;  angles 
9„  denote  the  orientations  of  B  to  the  local  nominal  shock  normal 
at  those  times. 

The  geometry  of  case  a,  corresponding  to  a  beam  of  reflected  ions. 
The  three-dimensional  diagram  of  the  bow  shock  at  left  shows  the 
B-X  plane,  its  intersecting  curve  C,  and  the  observation  point  0. 

Tic  marks  are  at  5  intervals.  The  close-up  at  right  of  cross 

section  C  surrounding  the  observation  point  shows  that  the  point 
was  near  the  forwardmost  origin  of  return-protons  of  3  KeV  guid¬ 
ing  center  energy. 

The  geometry  of  case  b,  for  a  reflected  beam,  showing  both  2  keV 
and  3  keV  boundaries  relative  to  the  observation  point. 

The  geometry  of  two  cases,  c  and  h,  in  which  no  return  ions  were 
detected.  The  observation  points  were  outside  the  3  keV  boundaries. 

The  geometry  of  two  cases,  f  and  g,  when  diffuse  return  ion  dis¬ 
tributions  were  recorded. 

The  geometry  of  two  cases,  d  and  e,  when  intermediate  return  ion 
distributions  were  recorded.  The  observat ion  points  were  within 
the  respective  1  keV  boundaries. 


Page  G. 


Figure  7 


A  representation  of  an  inferred  foreshock  structure  in  the  B-X 
plane.  The  cases  (lettered  circles)  are  indicated  as  if  the 
structure,  i.e.,  the  IMF  had  been  fixed  and  the  observation  point 
had  moved  around  the  shock. 
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Figure  ),  The  magnitude  and  components  of  the  interp Janetary  magnetic 
field  (IMF)  during  the  upstream  ion  observations  of 
Paschmann  et_  a_l_.  (1979,  Figure  3).  Shaded  sections  cover 
Intervals  during  which  the  three  classes  of  return  particle 
spectra,  reflected.  Intermediate,  and  diffuse,  were  mea¬ 
sured,  The  lower  case  letters  mark  times  selected  for  the 
examples  in  the  ensuing  figures;  angles  0gn  denote  the 
orientations  of  B  to  the  local  nominal  shock  normal  at  those 
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Figure  2.  The  geometry  of  case  a,  corresponding  to  a  beam  of 
reflected  ions.  The  three-dimensional  diagram  of 
the  bow  shock  at  left  shows  the  B-X  plane,  its 
Intersecting  curve  C,  and  the  observation  point  0. 
Tic  marks  are  at  5  Rg  intervals.  The  close-up  at 
right  of  cross  section  C  surrounding  the  observa¬ 
tion  point  shows  that  the  point  was  near  the  for- 
wardmost  origin  of  return-protons  of  3  keV  guiding 
center  energy. 


Figure  k.  The  geometry  of  two  cases,  c  and  h.  In  which  no 
return  Ions  were  detected.  The  observation 
points  were  outside  the  3  keV  boundaries. 


The  geometry  of  two  cases,  d  and  e,  when  Intermediate 
return  Ion  distributions  were  recorded.  The  obser¬ 
vation  points  were  within  the  respective  1  keV 
boundaries. 
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TRW  Defense  6  Space  Systems  Group 
One  Space  Park 

Redondo  Beach,  California  90278 


Medium-period  geomagnetic  pulsations  are  generated  In  the  interaction 
between  the  solar  wind  and  the  magnetosphere,  in  the  magnetosphere,  or  in  both, 
and  their  characteristics  are  controlled  by  conditions  in  the  solar  wind,  in 
the  magnetosphere,  or  in  both.  Determination  of  what  the  correct  relationships 
are  between  pulsations  and  the  geomagnetic  environment  is  currently  under  study 
by  several  groups,  following  older,  pioneering  work  by  Soviet  investigators 
(Troitskaya  et  al.,  197a;  Gul'el'mi  et_  a]_. ,  1973). 

There  is  little  doubt  that  some  connection  exists  between  the  solar  wind 
and  pulsation  activity  at  earth's  surface.  An  example  is  shown  in  Figure  I.  These 
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Figure  I 

results  by  Saito  et  al .  (1979)  with  hourly  solar  wind  and  pulsation  data  (left 
panels)  virtually  dupTicate  those  by  Greenstadt  et_  al_.  (1979)  and  Wolfe  et  al . , 
(1978).  Pc  3  amplitude  rises  with  Vjy,  top  left,  and  with  decreasing  cone 
angle  bottom  left.  The  inverse  correlation  with  0w»  clearly  improves 
when  Bgg  <  SO*,  but  strong  scatter  has  affected  all  studies.  Cone  angle  0xg 
Is  the  angle  between  the  IMF  and  the  solar  ecliptic  X-axis.  Part  of  the  scatter 
may  arise  from  the  time  scale  of  IMF  orientation  changes,  which  is  generally 
much  less  than  an  hour.  Hence  the  daily  data  in  the  right  panels  show  little 
of  Interest. 
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evidence  that  the  noon-midnight  boundary  is  an  important  source  of  magneto- 
spheric  waves. 

2.  Quantitative  development  of  signal  model.  The  signal  model  needs  to 
be  developed  further  to  produce  quantitative  estimates  that  can  be  observation- 
ally  tested. 

3.  Definition  and  refinement  of  surface-based  "input  level."  Some  means 
of  conditioning  and  combining  pulsation  records  from  ground  stations  is  needed 
that  will  characterize  the  global  pulsation  activity  level  as  ap  does  the  global 
magnetic  storm  activity  level. 

4.  Correlation  of  input  level  with  solar  wind  parameters.  The  global 
pulsation  measure,  once  defined,  should  be  used  to  refine  and  improve  the  one- 
station  correlations,  such  as  those  in  Figure  1. 

5.  Tracing  of  signals  through  plasma  regimes.  Spacecraft  observations 
should  be  used  in  an  attempt  to  actually  document  the  passage  of  specific  sig¬ 
nals  from  the  solar  wind  or  the  magnetopause  to  the  ground. 

A  statistical  effort  to  link  waves  in  the  magnetosheath  to  pulsations 
via  the  same  relationships  already  studied  has  been  made  recently.  In  Figure  3 
we  see  that  the  correlation  of  signal  amplitude  with  cone  angle  has  been  ex¬ 
tended  by  Saito  et^  a]_.  (1979)  to  the  magnetosheath,  central  panel,  where  sig¬ 
nal  power  is  highest  for  lowest  0j(B»  Signal  power  was  also  found  to  correlate 
inversely  with  angular  distance  a  from  the  subsoiar  line  in  the  magnetosheath, 
right  panel.  Such  statistical  correlations  need  to  be  enlarged,  and  particu¬ 
lar  events  in  the  sheath  should  be  sought  and  studied. 


9»-  aa-OEPENO.  OF  P«.  IN  THE  MAGNETOSHEATH 
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The  Air  Force’s  magnetometer  system  may  be  expected  to  make  significant 
contributions  to  items  3,  4,  and  5,  above.  One  necessity  for  improved  data 
on  pulsation  origin,  both  for  statistical  and  individual -case  approaches,  is 
a  set  of  ground  stations  insulated  from  the  strongest  sources  of  geomagnetic 
variability  but  responsive  to  wave  signals.  Some  latitude  may  be  much  better 
than  others,  as  illustrated  in  Figure  4.  Pulsation  period,  which  varies  with 
geomagnetic  latitude,  at  left,  is  also  subject  to  dynamic  variation  caused  by 
changing  plasma  density  in  the  plasmatrough.  This  dynamic  variability  affects 
pulsation  amplitudes  through  changing  resonance  periods  and  is  probably  re¬ 
sponsible  for  much  of  the  scatter  in  correlation  diagrms.  Diurnal  contri¬ 
butions  to  the  plasmatrough  variation  are  seen  if  Figure  4  at  right,  where 
characteristic  density  enhancements  are  visible  at  certain  local  times.  Both 
figures  are  from  Orr  (1979).  The  observatory  symbols  under  the  left  figure 


Figure  4 


show  the  locations  of  the  Air  Force  chain  (AF)  and  most  of  the  stations  re¬ 
sponsible  for  the  solar  wind  parameter  correlations.  The  stations,  or  station 
chains,  are:  Borok;  Pittsburg,  N.  H.;  Great  Britain;  Calgary;  Leduc;  Ft. 
Reliance.  We  may  infer  that  the  Air  Force  chain  is  advantageously  situated 
for  monitoring  Pc  4  comparatively  free  of  plasma  trough  variability.  Ft. 
Reliance,  or  any  station  at  about  70*,  may  also  be  a  good  site  for  monitoring 
Pc  5  in  the  afternoon. 

The  AFGL  system  should  be  equipped  with  the  computing  and  communication 
hardware,  software,  and  staff  necessary  to  take  advantage,  through  internal 
or  external  usage,  of  its  unique  potential. 
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APPENDIX  I 


LETTERS  TO  DR.  SAGALYN 


TRW 


79-4351 .6-59 
Bldg  R-l ,  Rm  1176 
19  April  1979 


Dr.  Rita  Sagalyn,  Chief 

Plasmas,  Particles  &  Fields  Branch 

AFGL 

Hanscom  AFB,  MA  01371 
Dear  Dr.  Sagalyn: 

Thank  you  for  the  opportunity  to  attend  AFGL's  geomagnetic  workshop,  6-7 
April.  I  thought  the  workshop  informative,  productive,  and  worthy  of 
repetition.  What  more  can  be  said? 

Here  are  my  answers  to  the  items  on  your  list  of  questions.  Each  reply  Is 
cast  in  the  frame  of  my  own  project  rather  than  in  the  full  generality  of 
the  questions,  which  I  presume  are  best  answered  in  total  by  a  panel  or  a 
composite  of  such  individual  responses. 

To  define  my  own  project  briefly,  I  am  attempting  to  determine  wheter  geo¬ 
magnetic  pulsations  can  be  related  to  solar  wind  parameters  with  sufficient 
precision,  i.e.,  reproducivility,  to  allow  either  phenomenon  to  serve  as  mea¬ 
sure,  or  even  predictor,  of  the  other.  The  attempt  involves  study  of  both 
the  physical  processes  that  may  be  Involved  in  connecting  the  solar  wind  to 
the  ground  and  of  "rule  of  thumb"  methods  of  correlation  that  may  permit 
utilization  of  some  sort  of  index  while  understanding  of  physical  processes 
is  still  incomplete. 

There  is  a  body  of  experimental  and  model  investigation,  many  years  old,  that 
supports  the  existence  of  a  connection  between  solar  wind  properties  and  Pc 
3,4,5  properties,  but  the  connection  Is,  with  rare  exception,  tenuous  and 
unreliable  because  of  the  multivariate  nature  of  the  problem  and  associated 
methodological  difficulties. 

The  approach  I  see  Is  to  assume  that  pulsation  amplitude  «B  (at  a  point  on 
the  surface)  Is  a  function  p  of  several  variables: 

6B  =  0 (T ,  0^g  Vjjy,  t|_»  aM’  Bz>» 

where  T,  eXB,  Vsw,  t^,  AM,  and  Bz  are  the  pulsation  period,  IMF  cone  angle, 
solar  wind  velocity,  local  time,  magnetic  latitude,  and  north-south  IMF 
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component,  respectively.  The  period  in  turn  is  itself  a  function  t  of  some 
of  the  same  variables  for  possibly  independent  reasons: 

T  =  t(t|^,  A Jy|  ,  B2»  Bjy), 

where  Bsw  is  the  IMF  magnitude.  I  use  Bz  merely  as  a  shorthand  external  indi¬ 
cator  of  geomagnetic  activity  level.  One  task,  for  the  present,  is  to  find 
the  means,  methodologically,  to  hold  as  many  variables  constant  as  possible 
while  searching  for  improved  correlations  among  the  remainder.  The  AFGL 
chain  holds  promise  of  yielding  data  stabilizing  t^,  am,  and  t  (t|_,  A^,  , 

considerably  simplifying  the  investigation. 

Now  the  answers  to  your  questionnaire: 

l.a  Two  types,  more  precisely,  two  scales  of  information  are  needed. 
First,  relatively  brief,  say,  two-hour  to  two-day  Intervals  of  data  correspond 
ing  to  identified  events  or  sequences  of  events  are  needed  for  the  study  of 
the  physical  transmission  of  waves  through  the  magnetosphere  to  the  ground. 
Second,  relatively  lengthy  Intervals,  say  several  hours  of  data  each  day  for 
several  months,  are  needed  to  develop  and  test  pulsation-based  indices  of  mag- 
netospheric  and  solar  wind  parameters.  Within  either  scale,  the  Information 
content  consists  of  the  time  and  amplitude,  or  power,  of  signals  In  selected 
frequency  bands. 

l.b  The  information  can  be  made  available  from  fluxgate  or  induction 
coll  recordings  as  plots,  films,  tapes  (or  listings),  depending  on  the  appli¬ 
cation.  Plots  and  tapes  can  be  made  available  by  mall;  taped  data  can  also 
be  made  available  by  telephone  line  between  computers. 

I.c.i.  Mid-latitude  networks,  i.e.,  stations  equatorward  of  the  plasma- 
pause  like  AFGL's,  are  advantageously  situated  to  escape  some  of  the  extremes 
in  variability  of  pulsation  characteristics  that  affect  stations  connected  to 
the  plasmatrough.  The  most  outstanding  example  is  periodicity:  a  station  at, 
say  60°  magnetic  latitude  may  pick  up  a  resonant  response  In  the  low  end  of 
the  Pc  3  band,  at  a  Pc  3  amplitude,  one  day  and  the  low  end  of  the  Pc  5  band, 
at  a  Pc  5  amplitude,  the  next  day  because  of  changes  In  density  In  the  plasma- 
trough.  How  can  such  measurements  be  compared?  Badly. 
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To  first  order,  the  plasmasphere  above  the  AFGL  chain  tends  to  be  more 
stable  than  the  trough  further  north  so  that  the  stations  ought  to  pick  up 
signals  or  resonances,  if  any,  concentrated  in  the  Pc  4  band.  That  should 
help,  because  the  amplitude  in  the  one  band,  from  day  to  day,  may  be  more 
directly  relatable  to  differences  in  the  level  of  input  perturbation,  or 
signal,  at  the  magnetopause,  and  less  to  local  variations  in  the  magneto¬ 
sphere.  We  do  not  belittle  the  prospect  of  getting  one  variable  out  of  the 
way  in  this  business.  Further,  if  this  approach  works,  it  may  become  pos¬ 
sible  to  calibrate  some  of  the  more  refractory,  high  latitude  data  so  as  to 
make  them  more  usable. 

l.c.i.i.  Additional  reduction  in  variability  can  be  achieved  at  certain 
local  times,  so  an  E-W  chain  increases  the  probability  that  some  station  will 
be  in  the  optimal  location  often  enough  to  build  up  an  acceptable  population 
of  cases.  This  probability  becomes  even  higher  if  Eskdalemuir  (U.K.),  for 
example,  is  considered  as  an  eastern  extension  of  the  AFGL  network. 

l.c.i.i.i.  For  the  explicit  purposes  I  am  enunciating,  a  sampling  in¬ 
terval  of  20  seconds  and  a  processing  interval  of  five  to  ten  minutes  to  ob¬ 
tain,  for  example,  maximal  amplitudes  or  spectral  estimates,  would  be  adequate. 
This  is  relatively  undemanding  and  is  well  within  the  resolution  of  the  AFGL 
system. 

1. d  Ground-satellite  coordinated  measurements  are  necessary  both  for 
studying  cases  and  for  statistical  analysis.  If  the  AFGL  chain  yields  a 
decent  result,  then  network-network  coordination,  particularly  between  the 
AFGL  and  other  stations  at  about  70°  would  be  desirable. 

2.  Pulsations  can  be  monitored  by  visual  Inspection  of  filtered  or  un¬ 
filtered  plots,  by  representations  of  filtered  signals  In  terms  of  amplitude 
levels,  power  spectral  densities,  or  dynamic  spectra  (in  order  of  Increasing 
cost,  complexity,  and  Information  content).  Such  representations  can  be 
displayed  as  plots  or  recorded  on  tape  for  further  processing  and  analysis. 
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Some  types  of  pulsations,  PI  2  for  example,  may  be  predictive  of  substorms 
or  even,  we  may  speculate,  of  geosynchronous  satellite  anomalies.  Predictions 
of  this  sort  are  In  their  Infancy  and  require  Investigation  and  development. 

3. a  As  far  as  pulsations  are  concerned,  the  Inadequacy  of  existing  In¬ 
dices  lies  In  the  current  lack  of  such  Indices.  Years  ago,  Salto  devised  a 
pulsation  Index  he  called  Kc  3,  based  on  Pc  3  measurements  at  one  observatory, 
which  correlated  well  with  solar  wind  velocity.  This  development  was  never 
followed  up  and  never  expanded  Into  a  worldwide  Index.  The  IMS  network  offers 
an  opportunity  and,  I  believe,  an  obligation  to  attempt  such  a  project. 

3. b  New  Indices  should  be  devised,  but  It  will  not  be  an  easy  task. 

It  must  begin  as  simply  as  possible  with  Improved  correlations  that  will  lead 
to  selection  or  establishment  of  stations  most  suitable  for  a  worldwide 
monitor  system.  It's  quite  possible  that  the  Air  Force  chain  will  serve  as 
the  core  bases  for  recording  data  that  can  be  fashioned  Into  a  worldwide  Index, 
or  at  least  a  mid-latitude  Index. 

4. a.&b.  For  the  particular  Investigation  and  application  I  have  In  mind, 
satellite-ground  calibration,  and  the  reverse,  are  Intrinsic  elements. 

5.  Yes  and  No.  Existing  methods  are  adequate,  but  the  AFGL  system 
does  not  appear  at  present  to  be  set  up  to  use  them.  First-order  processing, 
l.e.,  spectral  selection  of  data  has  not  yet  been  developed  for  either  Inter¬ 
nal  or  external  use,  and  phone-line  callup  capability  has  not  been  established 
for  external  use.  There  did  not  appear  to  be  a  research  plan  at  AFGL,  whose 
unfolding  would  automatically  result  In  the  above  facilities.  For  further 
detail,  see  commentary  below. 

6.  a  Current  ground  based  plans  seem  to  be  mixed.  Some  recordings  have 
already  been  terminated;  some  will  be  discontinued  at  the  end  of  the  IMS; 
some  after  MAGSAT;  some  as  budgets  dictate.  Some  stations  are  just  now  de¬ 
signing  Improvements  for  the  future,  signifying  an  Intention  to  continue 
recording  Indefinitely. 
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6.b  What  Is  required  Is  largely  unknown  for  the  simple  reason  that 
there  Is  a  sizeable  phase  lag  between  recording  and  analysis  of  IMS  data. 
Serious  Inadequacies  of  IMS  data  collection.  If  they  exist,  may  not  become 
evident  for  a  year  or  two.  One  of  the  major  accomplishments  of  the  IMS,  In 
which  the  AFGL  chain  shares.  Is  a  great  Improvement  In  acquisition  and  re¬ 
cording  technique,  which  now  Is  ahead  of  the  capacity  of  the  analysis  com¬ 
munity  to  absorb  the  data.  Although  several  significant  results  have  al¬ 
ready  come  out  of  the  IMS,  the  bulk  of  advances  In  both  quantity  and  quality 
will  undoubtedly  arise  In  the  early  80 's. 

My  view  Is  that  recordings  should  continue  toward  solar  maximum  until  analysis 
shows  which  deletions  or  changes  will  be  profitable.  Consider  the  Air  Force 
chain  and  pulsations  as  an  example:  recent  work  has  shown  the  global  nature 
of  some  selected  pulsation  events;  that  Is,  several  stations,  both  merldlanally 
spaced  and  conjugate  to  each  other,  observed  similar  manifestations  at  the 
same  time.  Some  of  the  data  were  displayed  during  the  workshop.  The  distance 
limits  on  such  similarity  have  not  been  established.  It  may  be  (It  may  not) 
that  the  eight  AFGL  stations  could  be  replaced  by  four  with  no  serious  loss 
of  Information  content,  but  I  don't  think  we'll  know  until  well  after  the 
IMS  Is  over  (see  below). 


COtMENTARY 

AFGL  has  Installed  an  advanced  system  of  geomagnetic  data  collection  In 
which  the  acquisition  points  are  advantageously  located  for  at  least  some  1m 
portant  geomagnetic  Investigations.  Further,  the  data  gathering  and  storage 
scheme  creates  a  file  of  data  of  superior  organization  and  accessibility. 

The  file  so  created  contains  gargantuan  amounts  of  Information  which  can  be 
extracted  only  by  advanced  techniques  of  data  processing.  Fortunately  such 
techniques  are  available  today,  and  at  far  less  expense  than  they  were  even 
a  very  few  years  ago. 
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I  was  surprised,  and  disappointed,  to  find  that  certain  kinds  of  pro¬ 
cessing  had  not  been  planned  as  an  Integral  part  of  the  AF6L  facility  from 
the  outset.  One  purpose  of  the  system  Is  to  record  and  apply  field  oscil¬ 
lations,  but  the  basic  content  of  the  recordings  from  which  any  analysis 
begins  Is  a  summary  of  what  oscillations  are  present  at  a  given  time.  Even 
If  the  spectral  content  Is  not  stored  separately,  and  I  think  It  definitely 
should  be,  the  capacity  to  comb  any  given  set  of  tapes  for  spectral  content 
should  be  a  routine  library  function  fundamental  to  virtually  any  exploit¬ 
ation  of  the  pulsation  records,  whether  Internally  or  externally  originated. 

In  a  sense,  the  entire  AFGL  net  may  be  regarded  as  a  single  geomagnetic 
Instrument  whose  application  requires  an  associated  calibration.  I  offer 
two  Illustrations:  suppose  someone  Is  Interested  In  studying  a  selected 
storm  or  substorm  or  a  set  of  them  —  some  IMS  Intervals,  say.  The  dis¬ 
turbance  fields  would  be  obtained  by  subtracting  Sq.  But  what  Is  for  each 
of  the  AFGL  observatories?  Any  Individual  wishing  to  examine  a  particular 
interval  can  hardly  be  expected  to  make  a  bulk  analysis  of  the  entire  library 
to  find  Sq  before  he  can  begin  work  on  a  two-day  storm.  A  set  of  Sq  curves 
or  tables  should  be  regarded  as  a  calibration  that's  part  of  the  system  des¬ 
cription.  Similarly,  there  is  a  diurnal  variation  in  baseline  pulsation  level 
that  can  be  used  to  Indicate  the  best  hours  for  selecting  passband  data  and 
the  level  above  which  pulsation  events  can  be  said  to  have  occurred.  Curves 
and  tables  of  such  variations  should  be  part  of  the  system  description. 

RECOMMENDATIONS 

It  may  be  my  misunderstanding,  but  I  though  I  received  conflicting 
evaluations  of  the  resources  available  for  data  processing  at  AFGL.  Personnel 
at  the  Saturday  afternoon  pulsations  workshop  Indicated  that  the  kind  of  data 
conditioning  I  have  mentioned  would  be  beyond  those  resources.  I  understood 
from  conversation  with  you,  however,  that  It  might  be  possible  to  support 
some  data  manipulation.  Naturally,  I  prefer  your  Interpretation.  I  recom¬ 
mend  an  Investment  In  data  conditioning,  processing,  description,  and 
comnunlcatlon,  as  follows: 
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1.  Library  routines,  or  microprocessor  hardware,  for  digital  filter¬ 
ing  of  the  data  with  selectable  passbands,  start,  and  stop  times  should  be 
developed  and  debugged,  or  bought,  with  the  utmost  dispatch. 

2.  A  capability  of  recording  filtered  output  on  tape,  printed  list¬ 
ings,  and  plotted  graphs  should  be  acquired.  I  really  see  no  reason  AFGL 
shouldn't  be  accumulating,  storing,  and  perhaps  even  mailing  out,  the  same 
sort  of  f-t  diagrams  that  Onagawa  and  Kakloka  are  providing.  I  am  attaching 
a  copy  of  a  sample  page  from  Onagawa. 

3.  Average,  quiet-time  characteristics  of  each  station  at  all  acces¬ 
sible  frequencies  should  be  discovered,  recorded,  listed,  and  printed  as  part 
of  the  system  description. 

4.  A  capability  of  telephone  line  coupling  that  would  permit  external 
computer  access  to  any  stored  data  set  should  be  established. 

In  sum,  there  Is  hope  and,  I  believe,  expectation  that  the  Air  Force  mag¬ 
netometer  net  can  be  a  productive  element  In  geomagnetic  research  and  appli¬ 
cation,  but  some  definite  steps  need  to  be  made  now  to  turn  promise  Into 
reality. 

Sincerely, 


Eugene  W.  Greens tadt 

Space  Sciences  Department 

TRW  Defense  &  Space  Systems  Group 

EWG:mjl 

enc:  as  Indicated 
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Dr.  Rita  Sagalyn 
Air  Force  Geophysics  Laboratory 
Air  Force  Hanscom  Field 
Sudbury,  Massachusetts  01776 

Dear  Rita: 

Thank  you  very  much  for  the  opportunity  to  participate  In 
the  Geomagnetic  Field  Workshop  held  at  the  Air  Force  Geophysics 
Laboratory.  I  think  the  meeting  was  quite  worthwhile  and  I  hope 
our  suggestions  will  be  of  use  to  you. 

My  Impression  of  the  magnetometer  network  was  very  favorable. 
The  network  seems  to  be  working  well  and  the  quality  of  the  data 
returned  Is  high.  The  mini-computer  system  supporting  the  network 
ts  quite  satisfactory.  The  data  presently  being  acquired,  as  well 
as  that  previously  archived,  should  be  of  great  value  In  a  variety 
of  research  projects. 

I  must  say  that  I  am  disappointed  that  the  laboratory  Is  not 
making  greater  use  of  the  network  data  than  It  is.  This  is  sur¬ 
prising  In  as  much  as  the  Air  Force  has  need  for  the  kind  of  In¬ 
formation  this  network  can  provide.  In  fact,  the  Air  Force  Is 
supporting  the  Space  Environment  Laboratory  In  Boulder  to  satisfy 
such  needs. 

I  am  not  sufficiently  familiar  with  the  organization  or  goals 
of  the  AFGL  to  know  the  responsibilities  of  your  group.  However, 

I  suspect  you  are  obligated  to  study  geophysical  phenomena  with 
the  eventual  goal  of  quantl tatl vely  predicting  environmental 
parameters  on  a  basis  of  a  few,  simply-monitored  variables.  An 
example  of  what  I  mean  would  be  predicting  the  size  of  the  auroral 
oval  from  polarcap  magnetic  variations.  It  does  not  appear  that 
the  laboratory  Is  very  deeply  Involved  In  such  activities  at  the 
present  time,  since  no  one  at  AFGL  seems  to  have  published  any 
papers  on  substorms  In  the  last  four  years. 

I  choose  this  criteria  for  evaluation  because  It  has  recently 
become  clear  that  the  fundamental  element  of  geomagnetic  activity 
If  the  magnetospherl c  substorm.  I  believe  most  transient  changes 
In  the  Ionospheric  or  magnetospherlc  environment  of  Interest  to  the 
Air  Force  can  be  related  to  substorms  In  some  way.  While  I  see 
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great  Interest  at  AFGL  in  magnetic  indices  and  their  relationship 
to  Ionospheric  parameters,  no  person  there  appears  to  be  specifi¬ 
cally  studying  the  physical  processes  involved  in  substorms.  As 
a  consequence,  no  one  has  been  able  to  define  uses  for  the  data 
from  the  magnetometer  network. 

To  define  possible  uses  for  the  network  and  its  data  I  must 
try  to  answer  the  first  question  in  your  list  of  discussloh  topics, 
i.e.,  what  are  some  of  the  important  unsolved  problems  of  magneto- 
spheric  physics.  In  my  opinion,  these  Include: 

1)  What  causes  the  onset  of  a  substorm  expansion  phase? 

2)  What  controls  where  the  substorm  will  occur,  how  big 

tt  will  be,  how  long  It  will  last? 

3)  What  determines  the  properties  of  particles  energized 

and  Injected  by  the  substorm? 

4)  What  are  the  mechanisms  for  redistributing  Injected 

particles? 

.5)  What  are  the  loss  mechanisms  for  these  particles? 

6)  How  is  the  energy  released  by  the  substorm  dissipated 

fn  the  ionosphere? 

7)  What  effects  does  this  energy  deposition  have  on  the 

Ionosphere? 

As  you  know,  many  of  us  working  on  the  subject  of  substorms 
have  speculative  answers  to  most  of  these  questions,  but  based  on 
Inadequate  data.  You  also  know  these  models  are  quite  controversial 
and  are  certainly  far  from  being  quantitative. 

With  regard  to  your  second  question  of  what  information  Is 
needed  to  solve  these  problems,  I  suggest  the  following: 

1)  Better  monitoring  instruments  with  higher  time  and  space 

resolution. 

2)  Better  distribution  of  instruments  such  as  ground  networks 

and  multiple  satellites. 

3)  Better  data  acquisition  techniques. 

4)  Better  archiving  and  retrieval  systems. 

5)  More  support  of  researchers  studying  these  problems. 

As  a  specific  example,  I  believe  great  progress  could  be  made 
on  the  question  of  substorm  onset  mechanisms  if  we  had  an  auroral 
Imaging  device  producing  high  resolution  pictures  of  the  aurora  at 
a  rate  of  one  per  minute.  Using  solar  wind,  magnetospheric  and 
ground  network  data  obtained  simultaneously,  the  auroral  pictures 
would  help  Immensely  in  deciding  which  models  to  use  in  interpreting 
our  ground  data. 

To  answer  your  questions  of  what  contributions  the  AFGL  chain 
of  magnetometers  can  make  to  substorm  studies,  we  must  first  ask  In 
what  ways  is  the  chain  unique.  These  include: 
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1)  One  of  the  few  constant  magnetic  latitude  chains  in  the 

world. 

2)  Identical  instrumentation  at  all  stations. 

31  Instruments  have  high  sensitivity  and  low  noise. 

4)  Data  is  acquired  at  very  high  time  resolution. 

5)  Exceedingly  accurate  relative  timing  at  all  stations. 

6)  Data  are  acquired  and  stored  In  real  time. 

7)  There  is  an  extensive  historical  data  file. 

Geomagnetic  phenomena  which  can  be  studied  with  your  data  are  those 
which  produce  strong  magnetic  perturbations  at  sub-auroral  latitudes. 
These  include  substorm  field  aligned  currents  and  pulsation  phenomena 
associated  with  the  plasmapause  such  as  Pi  2  pulsations. 

Both  of  these  subjects  are  currently  of  great  interest  to 
magnetospheric  physicists  and  appear  to  be  closely  related  topics. 

The  onset  Oi*  the  expansion  phase  of  a  magnetospheric  substorm  is 
characterized  by  a  number  of  events  including  a  short  burst  of  Pi  2 
pulsations  and  the  formation  of  a  wedge  of  field  aligned  current 
diverting  a  portion  of  the  near  earth  tail  current  through  the 
auroral  oval.  There  is  some  evidence  that  the  PI  2  burst  is  the 
initial  transient  associated  with  the  resonant  properties  of  the 
circuit  carrying  this  field  aligned  current. 

Since  the  cause  of  the  substorm  expansion  onset  is  not  presently 
known,  it  is  important  to  Identify  the  processes  which  occur  near 
the  region  of  onset,  at  the  time  of  onset.  Pi  2  pulsations  can  be 
used  in  such  studies  because  they  define  expansion  onset  times  pre¬ 
cisely.  The  AFGL  network  with  its  high  sensitivity,  low  noise,  high 
time  resolution,  accurate  timing,  and  location  just  equatorward  of 
the  plasmapause  is  ideally  suited  for  observations  of  Pi  2  pulsations. 
One  Important  application  of  the  network  data  would  be  to  create  a 
list  of  substorin  expansion  onsets  for  use  in  substorm  research  projects. 

It  is  equally  Important  to  know  where  an  onset  occurs  and  how  it 
subsequently  develops.  Information  regarding  this  question  can  be 
obtained  from  the  east-west  magnetic  variations  (D  component)  measured 
along  a  constant  magnetic  latitude  chain.  The  center  of  the  sheet 
of  field  aligned  current  entering  the  ionosphere  in  the  morning  sector 
is  near  the  longitude  of  maximum  negative  deviation  in  D,  the  center 
of  the  outward  current  near  the  longitude  of  maximum  positive  deviation 
In  D,  and  the  center  of  the  current  systems  where  the  D  perturbation 
is  zero.  The  AFGL  network  Is  Ideally  located  for  determining  such 
parameters  when  North  America  passes  through  the  midnight  sector. 

The  preceding  comments  partially  answer  one  of  your  major 
questions,  how  can  magnetic  field  data  be  used  to  monitor  magnetic 
activity,  pulsations,  etc. 

Hy  preceding  comments  are  also  relevant  to  your  question  regard¬ 
ing  possible  new  indices.  Since  the  AFGL  is  of  limited  extent  in 
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local  time,  it  is  not  possible  to  generate  planetary  indices. 
However,  the  data  can  be  used  to  generate  indices  during  fixed 
intervals  of  universal  time.  For  example,  an  important  contribu¬ 
tion  the  network  data  could  make  would  be  to  produce  pulsation 
indices  for  various  types  of  magnetospheri c  wave  activity. 
Specifically,  a  Pi  2  index  of  substorm  activity  could  be  gener¬ 
ated  by  band  pass  filtering  out  the  60  second  period  band.  An 
indirect  monitor  of  solar  wind  characteristics  could  be  gener¬ 
ated  by  band  pass  filtering  the  Pc  3  band. 

The  final  question  I  would  like  to  touch  on  is  that  of  data 
dispersal.  I  must  say  that  it  is  here  that  the  AFGL  seems  to  be 
having  the  most  difficulty.  Until  very  recently,  few  persons 
outside  AFGL  had  ever  seen  the  data.  Further,  to  my  knowledge, 
very  little  had  been  done  with  It  Inside  AFGL.  The  reason  for 
these  problems  seems  to  be  psychological.  The  individuals  who 
have  set  up  the  system  consider  themselves  scientists  not  support 
staff.  Also  it  appears  that  the  pressure  on  these  individuals  to 
utilize  the  data  has  not  been  sufficient,  or  that  they  have  been 
otherwise  occupied.  Thus,  while  It  is  possible  to  acquire  small 
amounts  of  data  in  graphical  form,  a  request  to  study  a  major 
portion  of  the  digital  data  would  not  be  easy  to  satisfy. 

The  cheapest  and  easiest  way  for  AFGL  to  disperse  the  data 
would  be  to  duplicate  the  entire  data  file  as  it  is  and  send  it 
to  the  World  Data  Center  in  Boulder.  Alternatively,  AFGL  can 
act  as  a  data  center  filling  requests  as  they  come  in;  however, 
this  requires  staff  specifically  dedicated  to  this  function. 

Still  another  alternative  would  be  to  provide  a  dial-up  computer 
link  so  outside  users  could  access  the  on-line  data  files.  This 
method  is  limited  by  the  amount  of  data  kept  on  line,  by  the 
data  transmission  rates,  and  by  the  costs  of  dial-up  link.  It 
would  also  require  software  development  to  make  possible  the 
access  by  other  computers. 

Any  of  the  above  methods  is  unlikely  to  result  in  extensive 
use  of  the  data  file  because  of  the  costs  involved  in  using  the 
data.  To  obtain  such  use,  AFGL  will  have  to  support  the  cost  of 
research  on  the  data  base.  Possible  methods  include: 

1)  Reassignment  of  current  AFGL  staff. 

2)  Addition  of  new  AFGL  research  staff. 

3}  Visiting  scientist  positions. 

4)  Summer  research  positions. 

5)  Consulting  arrangements. 

6)  Grants  to  universities  or  research  organizations. 

If  the  data  are  to  be  analyzed  at  AFGL,  additional  support  must  be 
put  into  the  ml ni -computer  facility  In  the  form  of  dedicated  pro¬ 
grammers  and  additional  hardware  and  software.  Without  this,  AFGL' 
staff  or  visitors  will  not  be  able  to  access  the  data.  In  addition 
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Institutional  arrangements  must  be  made  whereby  those  Interested 
in  the  data  can  influence  program  development.  For  example,  in 
my  opinion,  a  major  problem  with  the  existing  data  processing 
scheme  is  that  suitable  averages  and  indices  are  not  generated 
by  the  data  acquisition  and  archiving  program  at  the  time  the 
data  are  acquired.  Since  I  pointed  out  the  need  for  this  four 
years  ago  while  we  were  building  the  instruments  at  UCLA  and  It 
was  not  implemented  by  those  responsible,  I  doubt  they  would 
accept  this  suggestion  today. 

An  attractive  alternative  is  for  the  AFGL  to  sponsor  research 
and  development  grants  In  outside  organizations.  For  example,  the 
algorithms,  programs,  etc.  required  to  analyze  AFGL  data  elsewhere 
could  be  implemented  at  AFGL  by  existing  staff  members  provided 
proper  documentation  of  these  was  generated  by  the  grantees. 

It  seems  likely  that  current  research  areas  in  geomagnetism 
will  eventually  become  matters  of  continuous  importance  to  both 
military  and  civilian  agencies.  When  this  happens,  continual 
monitoring  of  magnetic  activity  with  magnetometer  networks  like 
AFGL,  and  also  real  time  processing  and  display  of  data  from  these 
networks,  will  be  essential.  I  believe  one  important  goal  of  AFGL 
research  activities  should  be  to  incorporate  various  research  tools 
into  the  real  time  data  acquisition  system  so  that  their  value  in 
monitoring  can  be  demonstrated. 

I  hope  these  comments  and  suggestions  will  be  of  some  use  to 
you.  I  would  be  glad  to  discuss  any  of  these  topics  with  you  or 
your  laboratory  in  greater  detail.  Let  me  know  if  I  can  be  of 
further  help. 


Sincerely  yours,  ) 

Robert  L.  McPherron 
Professor  of  Space  Science 
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